Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,448,276 T→A 6.3% intergenic (‑127/+38) rplN ← / ← rpsQ 50S ribosomal subunit protein L14/30S ribosomal subunit protein S17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,448,2760TA 6.3% 352.4 / 7.6 121intergenic (‑127/+38)rplN/rpsQ50S ribosomal subunit protein L14/30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base T (62/51);  new base A (5/3);  total (67/54)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

GCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTC  >  NC_000913/3448128‑3448422
                                                                                                                                                    |                                                                                                                                                  
gCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAttt                                                                                                                                                   >  2:1044278/1‑150 (MQ=255)
   tAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGTGACCCTAAAATCAGTTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAATGATATGGGTAGAAAAAATAAACGGCTCATTTCTg                                                                                                                                                <  2:1546091/150‑1 (MQ=255)
   tAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTg                                                                                                                                                <  2:1787104/150‑1 (MQ=255)
   tAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTg                                                                                                                                                <  1:48180/150‑1 (MQ=255)
   tAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTg                                                                                                                                                <  2:1569341/150‑1 (MQ=255)
     gTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAg                                                                                                                                              <  2:440821/150‑1 (MQ=255)
           tCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGt                                                                                                                                          >  2:368262/1‑148 (MQ=255)
              tgtaTCATTTTAGTGCTCAGATAATGTAAACTATTACTTAGACCCGAAAATCAGTTGTTTAAATATCCCCTTATCGTGGGCGCGGCATTATAACACCGCTTCATGTATATGGGTAGAAAAAATAATCGGCTCATTTCTGAGCCGTTTAtt                                                                                                                                     <  2:1611921/147‑1 (MQ=255)
               ggATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTAtta                                                                                                                                    >  2:707429/1‑149 (MQ=255)
               ggATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTc                                                                                                                                    >  2:811829/1‑150 (MQ=255)
                  tCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCATGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTa                                                                                                                                 <  1:889883/150‑1 (MQ=255)
                  tCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTa                                                                                                                                 <  2:1701913/150‑1 (MQ=255)
                   cATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAAACAGGTCGTTAAAAATCCCCATATCGAGGGCGCGCCATTATAACACCGCTTCAAGGATATGGGGAGAAAAAATAACCGGCTCAGAAATGAGCCGTTTAttttttc                                                                                                                                <  1:1815578/150‑6 (MQ=255)
                     ttttAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGTCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                              <  2:22797/150‑1 (MQ=255)
                     ttttAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTg                                                                                                                              <  1:535394/150‑1 (MQ=255)
                             cTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtt                                                                                                                      >  1:188273/1‑149 (MQ=255)
                             cTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAgtg                                                                                                                      <  1:532651/150‑1 (MQ=255)
                               ccGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTa                                                                                                                    >  1:203630/1‑150 (MQ=255)
                               ccGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTa                                                                                                                    >  2:149671/1‑150 (MQ=255)
                               ccGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTa                                                                                                                    >  2:309866/1‑150 (MQ=255)
                                        tCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTAc                                                                                                           >  1:780339/1‑150 (MQ=255)
                                               ctactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCg                                                                                                    <  2:654920/150‑1 (MQ=255)
                                               ctactGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGGTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCg                                                                                                    >  1:220543/1‑150 (MQ=255)
                                               ctactCAGACCCGAAAATCACGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCCCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCAGTTCTGAGCCGTTTATTAGTATTGAGAGCGTGTAATCTATTACAGAACCg                                                                                                    >  2:1470289/1‑150 (MQ=255)
                                                 actGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCt                                                                                                  <  1:1795314/150‑1 (MQ=255)
                                                 actGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCATATTGAGAGAGTGTACTGTATTACAGAACCGCt                                                                                                  <  1:654431/150‑1 (MQ=255)
                                                  ctGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCtt                                                                                                 <  2:218227/150‑1 (MQ=255)
                                                    gagaCCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTtc                                                                                               <  1:1442482/150‑1 (MQ=255)
                                                      gaCCCGAAAATCAGGTCGTTAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTAAAGAACCGCTTTctct                                                                                            >  2:765040/1‑150 (MQ=255)
                                                        cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                           >  1:735284/1‑150 (MQ=255)
                                                        cccGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTa                                                                                           >  2:1135507/1‑150 (MQ=255)
                                                           gAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACaa                                                                                        >  2:1567982/1‑150 (MQ=255)
                                                                tCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGa                                                                                   <  2:1356223/150‑1 (MQ=255)
                                                                  aGGTGGTTAAAAATCCCAATATCGAGGGGGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAACCGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAAc                                                                                 <  1:559641/150‑1 (MQ=255)
                                                                       gatAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCg                                                                            >  1:439028/3‑150 (MQ=255)
                                                                       gTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCg                                                                            >  1:1709936/1‑150 (MQ=255)
                                                                        ttAAAAATCCCCATATCGAGGGCGCGGCATTATAACCCCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACCACGCGAACCAGCGt                                                                           >  2:1587116/1‑150 (MQ=255)
                                                                           aaaaTCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTcga                                                                        >  1:1122513/1‑148 (MQ=255)
                                                                            aaaTCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAg                                                                       <  1:1025751/150‑1 (MQ=255)
                                                                                ccccATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGAtt                                                                   >  2:1389557/1‑150 (MQ=255)
                                                                                 cccATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGAttt                                                                  <  2:1429316/150‑1 (MQ=255)
                                                                                    atatCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGt                                                               <  1:1754352/150‑1 (MQ=255)
                                                                                     tatCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTc                                                              <  2:801477/150‑1 (MQ=255)
                                                                                      atCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCt                                                             >  2:843066/1‑150 (MQ=255)
                                                                                        cGAGGGCGCGGCTTTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTg                                                           <  1:309866/150‑1 (MQ=255)
                                                                                        cGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTg                                                           >  1:550741/1‑150 (MQ=255)
                                                                                         gAGGGCGCGGCATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTgg                                                          <  1:838174/150‑1 (MQ=255)
                                                                                         gAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTgg                                                          <  1:984002/150‑1 (MQ=255)
                                                                                         gAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGAAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTgg                                                          >  2:1243816/1‑150 (MQ=255)
                                                                                          aGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGa                                                         <  2:1815138/150‑1 (MQ=255)
                                                                                           gggCGCGGCATTATAACACCTCTTCAATGATATGGGTAGAAAAAAAAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGTGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAAACAGCGTCCAAGACTAAATCCTGGAc                                                        <  1:895451/150‑1 (MQ=255)
                                                                                              cgcgGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGc                                                     >  1:1494971/1‑150 (MQ=255)
                                                                                                cgGCATTATAACACCGCTTCAAGTATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAAGACTAAATCCTGGAcgctgg                                                   <  2:138918/150‑6 (MQ=255)
                                                                                                cgGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGctg                                                   >  1:913826/1‑148 (MQ=255)
                                                                                                cgGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCgg                                                   <  2:118630/150‑1 (MQ=255)
                                                                                                   cATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACg                                                >  2:75608/1‑150 (MQ=255)
                                                                                                      tataACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCa                                             <  2:10978/150‑1 (MQ=255)
                                                                                                        taACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCAtt                                           >  2:1474610/1‑150 (MQ=255)
                                                                                                         aaCACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTc                                          >  2:3775/1‑150 (MQ=255)
                                                                                                           cacCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCGTTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACCACGCGCACCCGCTTCCAGGATTAAGACTTGGTCAGCGGCCGGCATTcgc                                        >  2:262483/1‑150 (MQ=255)
                                                                                                           cacCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGt                                                               <  1:1396897/127‑1 (MQ=255)
                                                                                                           cacCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGt                                                               >  2:1396897/1‑127 (MQ=255)
                                                                                                           cacCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcgc                                        >  2:886971/1‑150 (MQ=255)
                                                                                                           cacCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACAAGCGTCCAGGATTTAGTATTGGACAGCGGACGGCATTcgc                                        >  1:1281237/1‑150 (MQ=255)
                                                                                                            acCGCTTCTAGCTAATGGGTAGAAAAAATAAACGCCTCATATCTGAGCCGTTTATTCGTATTGAGTGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcgcg                                       <  2:269849/150‑1 (MQ=255)
                                                                                                            acCGCTTCAAGGATATGGGTAGAAAAAATAAATGGCTCTTTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcgcg                                       <  2:1613845/150‑1 (MQ=255)
                                                                                                            acCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcgcg                                       >  2:742258/1‑150 (MQ=255)
                                                                                                            acCGATTCAAGTATATGGGTAGAAAAAATACCACGCGCATTTCTGAGCCTTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTCCAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTcgcg                                       <  1:1197027/150‑1 (MQ=255)
                                                                                                             ccGCTTCAATGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCgg                                      <  1:1290213/150‑1 (MQ=255)
                                                                                                             ccGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCgg                                      >  1:1530011/1‑150 (MQ=255)
                                                                                                              cGCTTCAAGGATATGGGAAGAAAAAATAAACGGCTCATTTCTGCGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGa                                     <  1:1157457/150‑1 (MQ=255)
                                                                                                                cTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGAtt                                   >  2:1376309/1‑150 (MQ=255)
                                                                                                                 ttCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACTCGAACCAGCGTCCATGATTTAGTCTTGGACAGCGGACTTCATTCGCGGAttt                                  >  1:632423/1‑150 (MQ=255)
                                                                                                                 ttCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGAttt                                  >  2:396170/1‑150 (MQ=255)
                                                                                                                   cAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCa                                >  2:1247177/1‑150 (MQ=255)
                                                                                                                       gATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCa                            >  2:793097/1‑150 (MQ=255)
                                                                                                                         tatGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACg                          >  2:1640163/1‑150 (MQ=255)
                                                                                                                          atGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGt                         <  1:75608/150‑1 (MQ=255)
                                                                                                                          atGGGTAGAAAAAATAAACGGCTCATATCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGt                         <  1:225939/150‑1 (MQ=255)
                                                                                                                            gggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTcc                       >  1:504040/1‑149 (MQ=255)
                                                                                                                             ggTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAc                      <  2:694810/150‑1 (MQ=255)
                                                                                                                              gTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCAcc                     <  1:401155/150‑1 (MQ=255)
                                                                                                                              gTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACAAGCGTCCAGGATTTAGTCTTGGACAGAGGACGGCATTCGCGGATTTCAACCACGTCAcc                     >  2:894639/1‑150 (MQ=255)
                                                                                                                                 gAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGAt                  >  2:1719665/1‑150 (MQ=255)
                                                                                                                                 gAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGAAAGCGGACGGAATTCGCGGCTTTCAACCACGTCACCGAt                  >  1:296823/1‑150 (MQ=255)
                                                                                                                                   acgaaTAAACGGCTCAGAAATGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGCTAc                >  2:60240/4‑150 (MQ=255)
                                                                                                                                    aaaaTAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATAcc               <  2:1206489/150‑1 (MQ=255)
                                                                                                                                    aaaaTAAACGGCTCATTTCTGAGCCGTTTATTAGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGTACAGCGGACGGCATTCTCGGATTTCAACCACGTCACAGAAAcc               >  1:365370/1‑150 (MQ=255)
                                                                                                                                      aaTAAACGGCTCAGAAATGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGc             >  2:1482603/1‑150 (MQ=255)
                                                                                                                                        tAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCAt           >  1:462509/1‑150 (MQ=255)
                                                                                                                                        tAAACGGCTCAGAAATGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCAt           >  2:1624046/1‑150 (MQ=255)
                                                                                                                                         aaaCGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCAtt          >  2:582838/1‑150 (MQ=255)
                                                                                                                                         aaaCGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCAtt          <  2:51304/150‑1 (MQ=255)
                                                                                                                                          aaCGGCTCAGAAATGAGCGGTTTATTAGTATTGAGAGAGTGTACTGAATTACAGAACAGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCCCCGATACCGcttaa         >  2:929692/1‑146 (MQ=255)
                                                                                                                                          aaCGGCTCAGAAATGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGAGTTCAACCACGTCACCGATACCGCATTc         >  2:1094821/1‑150 (MQ=255)
                                                                                                                                              gCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCgttg     <  2:735284/150‑1 (MQ=255)
                                                                                                                                               cTCATTTCTGAGCCTTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACAGCTTTCTCTACAACGCGAACCCGCGTCCAGGAGGTAGTCTTGGACAGCTCACGGCATTCGCGGATTTATACCACCTCAACTATACCGCATTCgtggt    >  2:358456/1‑150 (MQ=255)
                                                                                                                                                 cATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCCCGTCCCCGATACCGCATTCGTTGTtc  >  2:292601/1‑150 (MQ=255)
                                                                                                                                                 cATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtc  <  1:1398909/150‑1 (MQ=255)
                                                                                                                                                 cATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTtc  >  1:1061485/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCATAGTCTGTTCTTGGATCATTTTAGTGCTCCGCTAATGTCAACTACTACTGAGACCCGAAAATCAGGTCGTTAAAAATCCCCATATCGAGGGCGCGGCATTATAACACCGCTTCAAGGATATGGGTAGAAAAAATAAACGGCTCATTTCTGAGCCGTTTATTCGTATTGAGAGAGTGTACTGTATTACAGAACCGCTTTCTCTACAACGCGAACCAGCGTCCAGGATTTAGTCTTGGACAGCGGACGGCATTCGCGGATTTCAACCACGTCACCGATACCGCATTCGTTGTTC  >  NC_000913/3448128‑3448422

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