Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,1110CT39.2% 95.0 / 33.9 74G219S (GGT→AGT) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base C (34/11);  new base T (0/29);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 7.56e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537968‑4538258
                                                                                                                                                   |                                                                                                                                                     
aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                              >  1:104935/1‑151 (MQ=255)
 aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                              >  1:170267/1‑150 (MQ=255)
       ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccg                                                                                                                                       >  1:530760/1‑151 (MQ=255)
            tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGCTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACtgccgc                                                                                                                                   >  1:154171/1‑150 (MQ=255)
                ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTACATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                                >  2:397076/1‑149 (MQ=255)
                ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                                >  2:332419/1‑149 (MQ=255)
                 cATCTGTTAATGAGACGGAAGCAACCTTATCCCTTTTAAAATTATTCCCGGGAAAATCAATTTAAATTAGGGCACCCGTTCGCAATCCGGGTTGGCCTCCCCCATTAAAAAGCCAGTTTTAATCCCTTTTATTCACAACCGCCGCACCAGc                                                                                                                             <  2:273976/151‑1 (MQ=255)
                 cATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCCCACCAGc                                                                                                                             >  2:79584/1‑151 (MQ=255)
                       gTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTcc                                                                                                                       >  1:253223/1‑151 (MQ=255)
                            gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAc                                                                                                                   >  2:199588/1‑150 (MQ=255)
                            gagaCGGGAGAAACCTTTTTCATTTTAAAATTAATACCGGTAAATTAAAGTAAAAATACGGTATCTGTTCACAATCCTGTTTTGCTTCCGCCATAAAAAAGCCGGTTTTTTCCTCTTTTTTTAACACCCGCCGCACCAGCCGTTCCGTa                                                                                                                    <  2:477427/149‑1 (MQ=255)
                                 gggAGCAAGTTTTTTCCATTTAAAATTATTCCCGGTAAAATAAATTTAAAATCCGGCTCCCTTTCGCAACCCTGTTTTGGTTTCCCCTTAAATAAGCCAGGTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAggg                                                                                                              <  1:731851/150‑1 (MQ=255)
                                  ggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  2:643006/1‑151 (MQ=255)
                                  ggAGAAAGTTTATTCCATTTAAATTTATTACCGGAAAAATAAATTTAAAATACGGCACCCGTTCGCAATCCTGGTTTGGCTTCCCCATTAATAAGCCGGTTTTTATCCCTTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            <  1:830727/151‑1 (MQ=255)
                                   gAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  1:395354/1‑150 (MQ=255)
                                        aGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      >  1:879893/1‑151 (MQ=255)
                                         gttttttCCATTTTAAATTTTTCCCGATAAAATAAATTAAAACACCGACACCCGTTCGCATTCGTGTTTTGGTTGACCCATTAATCGCCCGTTTTTTGCCCCTTTAATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       <  2:284335/144‑1 (MQ=255)
                                         gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  2:427224/1‑151 (MQ=255)
                                         gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  1:585335/1‑151 (MQ=255)
                                           ttATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  1:1003299/1‑149 (MQ=255)
                                             aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACCTCATCCGTTAGCAATTCTGGCTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGATCCGTACCAGGGCGATTCGCCAGCg                                                                                                 >  2:146664/1‑151 (MQ=255)
                                               tCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                                >  2:926221/1‑150 (MQ=255)
                                                 cATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                             >  2:291703/1‑151 (MQ=255)
                                                    tttAAAATTTTACCGGTTAAAATAAGTTAAAAACCGGCCTGCCTTCCGAAATCTGGGTTGGTCTTCGCCATAAATAGGCAGGTTTTGCCCCCTTTTTTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                            <  2:277875/150‑1 (MQ=255)
                                                       aaaTTATTCCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCCCCATAAATAAGCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       <  1:992296/151‑1 (MQ=255)
                                                       aaaTTATTCCCGGTGAAATAAATTTAAAATACGCCACCCTTCCGCAATCCTGTTTTGTCTTCCCCATTAATAAGCAAGTTTTTATCCCTTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                          <  1:610123/148‑1 (MQ=255)
                                                       aaaTTATTCCCGGTAAAATCAATTTAAAATACGGCTCCCGTCCCAAACCCTGTTTGGGCTTCCCCTTAAAAAACCCAGTTTTTACCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa                                                                                         <  1:964841/149‑1 (MQ=255)
                                                               aCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                               >  1:505174/1‑151 (MQ=255)
                                                                cctggaaaaTCAATTTAAAATACGCCACCCTTCCCAATCCTTGTTTTGGTTCCGCAATAAAAAGCCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                               <  2:168560/144‑1 (MQ=255)
                                                                ccGGTAAAATAAAGTTCAAATCCGGCCCCCGTCCGCAACCCTGTTTGGGCTCCCCCATTAATAACCCGGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              <  1:144441/151‑1 (MQ=255)
                                                                ccGGTAAAACAAAGTTAAAATCCGCCACCCGTCCGAATTCCGGTTTTGGCTCCCCATTTAAAAAGCCGGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                                <  1:336688/149‑1 (MQ=255)
                                                                   gaaaaacaaattacAAATCCGCCTTCCGTCCCCATCCCGGTTTGGCCTCCCCCATAAAAAACCCGGTTTTTCTCCCTTTTACCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                                <  2:70658/133‑1 (MQ=255)
                                                                      aaataaatTTCAAATACGGCACCCTTTCGCAACCCTGTTTTGGCTTCCCCATAAATAAGCCAGGTTTTATCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                          <  1:808844/149‑1 (MQ=255)
                                                                      aaaTCAATTTAAATTCGGGCACCCTTTCGAAATCTTGTTTTGTTTTCGCCTTAAAAAAGCCAGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                           <  1:947050/148‑1 (MQ=255)
                                                                      aaaTCAAGTTAAAATAGGCCACCCGTCCCCAACCTTGTTTGGGCTTCCCTATAAAAAAGCCAGTTTTTCTCCCTTTAATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                             <  2:757566/146‑1 (MQ=255)
                                                                                 aaaaACGGCACCCGTTCGCAATCCTGTTTTGTCTCCCCCATTAATAAGCCGGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                                <  1:682053/148‑1 (MQ=255)
                                                                                  aaTAGGGCATCTGTGCGCATTGCGTGTTGGTCTGCCCTATAAATAAGCCTGTTTTTACCCCTTTTCTTAAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                            <  2:782446/151‑1 (MQ=255)
                                                                                     aCGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                          >  2:743077/1‑150 (MQ=255)
                                                                                        ggCATCCGTTCGCAACCCTGTTTTGGCTTCGCCTTAAAAAACCCAGTTTTTATCCCTTTTATTCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                       <  1:66274/151‑1 (MQ=255)
                                                                                           cacccGTTCACACTCGTTTTTTGGTTTCCCCATTAATAAGCCAGTTTTTATCCCTTTTATTCCAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                        <  1:717083/145‑1 (MQ=255)
                                                                                              cctttcccatcccTGGTTTGGCTTCCCCATTAATAACCCGGTTTTTATCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                     <  2:416194/137‑1 (MQ=255)
                                                                                                 ttCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                >  1:681749/1‑150 (MQ=255)
                                                                                                   cGCCATCCTTGTTTGTCTTCCCCTTAAAAAACAGAGTTTTTACCTCTTTTCTTCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                              <  2:509154/150‑1 (MQ=255)
                                                                                                   cGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                                >  1:451184/1‑148 (MQ=255)
                                                                                                   cGCAATCCGGTTTTGGCTTCCCCATTAAAAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                             <  1:157828/151‑1 (MQ=255)
                                                                                                      aaTCCTGTTTTGGTTTCCCCATTAAAAAGCCGGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          <  1:789457/151‑1 (MQ=255)
                                                                                                       aTCCTGTTTTGCCCTCCCTATTAATACCCCAGTTTTTATCCCTTTAATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                           <  2:569578/149‑1 (MQ=255)
                                                                                                         ccTGGTTTGGTTTCGCCATAAAAAAGCCGGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          <  1:963886/148‑1 (MQ=255)
                                                                                                               ttggcctccGCCATTAAAAGGCCAGTTTTTACCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                   <  2:759409/142‑1 (MQ=255)
                                                                                                               ttggccTTCCCCTTAAATAACCCGGTTTTTACCCCTTTTATACAAATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                        <  2:254765/140‑1 (MQ=255)
                                                                                                               ttGGCTTCGCCATAAATAAGCCGGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  <  1:351009/150‑1 (MQ=255)
                                                                                                               ttGGCTCCCCCATTAATAACCCAGTTTTTATCCCTTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  <  2:64677/150‑1 (MQ=255)
                                                                                                               ttGCCTTCCCCATAAATAAGCCATGTTTTACCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                  <  2:928270/150‑1 (MQ=255)
                                                                                                                tGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                                >  1:526658/1‑151 (MQ=255)
                                                                                                                       cGCCATAAAAAACCCGGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           <  1:643276/150‑1 (MQ=255)
                                                                                                                         ccATTAATAGGCAGGGTTTAATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            <  2:253054/147‑1 (MQ=255)
                                                                                                                            ttAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      >  1:271354/1‑150 (MQ=255)
                                                                                                                            ttAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      >  1:748641/1‑150 (MQ=255)
                                                                                                                             ttaaaagcccAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      <  1:376267/142‑1 (MQ=255)
                                                                                                                             ttaaaaacccgttttttATCCCTTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      <  2:530508/137‑1 (MQ=255)
                                                                                                                              aataaCCCAGTTTTTATCACTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                    <  1:500195/150‑1 (MQ=255)
                                                                                                                                taaGCAAGTTTTTATCCCTTTTATTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  <  1:290076/150‑1 (MQ=255)
                                                                                                                                 aaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGAGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  >  1:384331/1‑149 (MQ=255)
                                                                                                                                  aGCCAGGTTTTATCACCTCTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt               >  1:605172/1‑151 (MQ=255)
                                                                                                                                     ccAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              >  2:701927/1‑150 (MQ=255)
                                                                                                                                      cAGTTTTTATCCCTTTTATCCACATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg              <  2:151824/149‑1 (MQ=255)
                                                                                                                                       aGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc           >  2:178537/1‑151 (MQ=255)
                                                                                                                                        gtttttACCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa          <  1:177851/149‑1 (MQ=255)
                                                                                                                                         tttttACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa          <  1:408867/149‑1 (MQ=255)
                                                                                                                                           tttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:549900/1‑151 (MQ=255)
                                                                                                                                           tttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  2:257702/1‑151 (MQ=255)
                                                                                                                                               tCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAttt   >  2:978633/1‑151 (MQ=255)
                                                                                                                                                cACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:438264/1‑151 (MQ=255)
                                                                                                                                                cACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:100202/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537968‑4538258

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: