Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,015,912 T→A 5.5% intergenic (+30/+40) rmf → / ← fabA ribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,015,9120TA 5.5% 1057.7 / 32.5 311intergenic (+30/+40)rmf/fabAribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase
Reads supporting (aligned to +/- strand):  ref base T (142/147);  new base A (8/9);  total (150/161)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 315 total aligned reads displayed.

ATCAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGA  >  NC_000913/1015764‑1016060
                                                                                                                                                    |                                                                                                                                                    
aTCAGGCCGGCATCGCCGGACGCTCACAAGAAATGTGGCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGCGAAGCCATGGCGGACAGGGTAGTACTGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAtt                                                                                                                                                     >  1:2422786/1‑150 (MQ=255)
     gCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGa                                                                                                                                                <  2:3626335/150‑1 (MQ=255)
     gCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGa                                                                                                                                                >  1:41945/1‑150 (MQ=255)
      ccGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACATTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACTGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAg                                                                                                                                               <  1:1433473/150‑1 (MQ=255)
       cGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAgg                                                                                                                                              <  2:2310177/150‑1 (MQ=255)
        ggCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGt                                                                                                                                             <  2:3910444/150‑1 (MQ=255)
        ggCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGt                                                                                                                                             <  2:4836588/150‑1 (MQ=255)
        ggCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGt                                                                                                                                             <  1:5226816/150‑1 (MQ=255)
        ggCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGt                                                                                                                                             <  2:142992/150‑1 (MQ=255)
        ggCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGt                                                                                                                                             <  2:1365273/150‑1 (MQ=255)
          cATCGCCTGACGCTCAAAATAAATGTGTCCGTATCTGACGCTGAATAAATGGTCACATTGGCTGGGCGGGTGGCTAGATGCCATGGCGGAAAGGGTAGTTATGGACTGATTGGGTCTATTTAAATAGAAACCTCCGCATTGCGGAGGttt                                                                                                                                           <  1:3667985/150‑1 (MQ=255)
            tCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCg                                                                                                                                         >  1:3076013/1‑150 (MQ=255)
            tCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTATTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCg                                                                                                                                         <  1:870547/150‑1 (MQ=255)
                 ggACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCttt                                                                                                                                    <  2:4135534/150‑1 (MQ=255)
                  gACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCtttt                                                                                                                                   >  2:2935802/1‑150 (MQ=255)
                      cTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGCCAGGGTAGTACTGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATa                                                                                                                               <  2:2431948/150‑1 (MQ=255)
                       tCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATAc                                                                                                                              >  1:4784037/1‑150 (MQ=255)
                       tCAAAAGAAATGTGTCCCTATCAGACGCGGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATAc                                                                                                                              >  2:2562940/1‑150 (MQ=255)
                         aaaaGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATActc                                                                                                                            >  2:3861664/1‑150 (MQ=255)
                            aGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgt                                                                                                                         >  1:4766123/1‑150 (MQ=255)
                            aGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgt                                                                                                                         <  2:4743857/150‑1 (MQ=255)
                              aaaTGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtct                                                                                                                       >  2:5052116/1‑150 (MQ=255)
                               aaTGTGTCCCTATCAGACGCTGAATCAAAGGTCACACTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtctg                                                                                                                      <  1:1382800/150‑1 (MQ=255)
                                aTGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAACAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGa                                                                                                                     <  1:402188/150‑1 (MQ=255)
                                  gtgtCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATCCTCTGTCTGAtt                                                                                                                   <  2:4264407/150‑1 (MQ=255)
                                     tCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATtata                                                                                                                >  1:976964/1‑150 (MQ=255)
                                        cTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATc                                                                                                             <  2:698413/150‑1 (MQ=255)
                                          ctCAGACGCTGAATCAAAGGTCACAAGGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAg                                                                                                           <  1:4911111/149‑1 (MQ=255)
                                          aTCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAg                                                                                                           <  1:3565297/150‑1 (MQ=255)
                                          aTCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAg                                                                                                           <  1:3909250/150‑1 (MQ=255)
                                             aGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAg                                                                                                        <  1:469958/150‑1 (MQ=255)
                                             aGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAg                                                                                                        <  1:5185233/150‑1 (MQ=255)
                                             aGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAg                                                                                                        >  1:739298/1‑150 (MQ=255)
                                               aCGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGc                                                                                                      >  1:668899/1‑150 (MQ=255)
                                                 gcggAATCAACGGTCACACTGGCTGGGAGGCTGGCGAGCAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAg                                                                                                    <  2:3141674/147‑1 (MQ=255)
                                                   tGAATCAACGGTCACAATGGCTGGGAGGCGGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGAc                                                                                                  <  1:3261424/150‑1 (MQ=255)
                                                   tGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGAc                                                                                                  >  1:1871657/1‑150 (MQ=255)
                                                   tGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGAc                                                                                                  >  1:1822940/1‑150 (MQ=255)
                                                      aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTa                                                                                               >  2:2844431/1‑150 (MQ=255)
                                                       tCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTAt                                                                                              >  1:2158557/1‑150 (MQ=255)
                                                          aaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAGTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCt                                                                                           <  2:3948063/150‑1 (MQ=255)
                                                           aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTg                                                                                          <  1:557418/150‑1 (MQ=255)
                                                             gTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGa                                                                                        >  1:2021296/1‑150 (MQ=255)
                                                              tCACACTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGaa                                                                                       <  1:2218260/150‑1 (MQ=255)
                                                              tCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGaa                                                                                       <  2:1496408/150‑1 (MQ=255)
                                                              tCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGaa                                                                                       >  1:3434766/1‑150 (MQ=255)
                                                              tCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGaa                                                                                       >  1:4454772/1‑150 (MQ=255)
                                                              tCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGaa                                                                                       <  2:313022/150‑1 (MQ=255)
                                                                  aagggCGGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGACTCTATAAAAAGAAACCGCCGCAGTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGa                                                                                   <  2:4263859/147‑1 (MQ=255)
                                                                  aaTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGa                                                                                   <  1:1037427/150‑1 (MQ=255)
                                                                  aaTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGa                                                                                   >  1:3797262/1‑150 (MQ=255)
                                                                      gCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGa                                                                               >  2:5288388/1‑150 (MQ=255)
                                                                        tGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGGTTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACt                                                                             <  2:2567159/150‑1 (MQ=255)
                                                                        tGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACt                                                                             >  1:435749/1‑150 (MQ=255)
                                                                         gggCGGCTGGCGAGAAGCCATGGCGGACTGGGTAGTAATGTCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGTTTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACtt                                                                            <  1:1228264/150‑1 (MQ=255)
                                                                          ggAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACttt                                                                           >  1:4572409/1‑150 (MQ=255)
                                                                            aGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCa                                                                         <  2:4662447/150‑1 (MQ=255)
                                                                            aGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCa                                                                         <  2:541205/150‑1 (MQ=255)
                                                                            aGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCa                                                                         <  2:5358903/150‑1 (MQ=255)
                                                                              gCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAgg                                                                       <  2:3347183/150‑1 (MQ=255)
                                                                                tGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTc                                                                     >  2:3277152/1‑150 (MQ=255)
                                                                                tGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTc                                                                     >  2:3488319/1‑150 (MQ=255)
                                                                                tGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGAGTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTATCAGGTc                                                                     >  2:4392946/1‑150 (MQ=255)
                                                                                 ggCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCg                                                                    >  1:1710203/1‑150 (MQ=255)
                                                                                    gagaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTg                                                                 >  2:118949/1‑150 (MQ=255)
                                                                                     agaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTgg                                                                <  1:3843932/150‑1 (MQ=255)
                                                                                      gaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGc                                                               <  1:4811891/150‑1 (MQ=255)
                                                                                      gaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGc                                                               <  2:446224/150‑1 (MQ=255)
                                                                                          ccATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGta                                                           >  2:4013240/1‑150 (MQ=255)
                                                                                            aTGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGtata                                                         >  1:5301245/1‑150 (MQ=255)
                                                                                               gCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGAt                                                      <  1:5052116/150‑1 (MQ=255)
                                                                                               gCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGAt                                                      <  1:2471545/150‑1 (MQ=255)
                                                                                                 ggACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCa                                                    <  1:2809354/150‑1 (MQ=255)
                                                                                                           gtaATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATc                                          >  1:3315767/1‑150 (MQ=255)
                                                                                                              aTGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAAc                                       <  1:288789/150‑1 (MQ=255)
                                                                                                               tGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAAcc                                      <  2:4948364/150‑1 (MQ=255)
                                                                                                                    tGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCAc                                 >  2:84530/1‑150 (MQ=255)
                                                                                                                    tGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATGATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGccc                                 >  2:110382/1‑148 (MQ=255)
                                                                                                                       ttCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTc                              <  2:304247/150‑1 (MQ=255)
                                                                                                                       ttCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTc                              >  1:3494184/1‑150 (MQ=255)
                                                                                                                        tCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCg                             <  1:2951066/150‑1 (MQ=255)
                                                                                                                              tctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATcc                       <  1:2343785/150‑1 (MQ=255)
                                                                                                                               ctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGTCGTATCGTGGAACTGACCGTCTTTCAGGTCGCTGGCGGTATAGTTCAGTCGACCATCAACCTGCACTTCGCCATCCg                      >  2:582212/1‑150 (MQ=255)
                                                                                                                               ctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCg                      >  2:161705/1‑150 (MQ=255)
                                                                                                                                tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGc                     <  1:867784/150‑1 (MQ=255)
                                                                                                                                tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGc                     >  2:647150/1‑150 (MQ=255)
                                                                                                                                tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGc                     >  2:5236085/1‑150 (MQ=255)
                                                                                                                                tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGc                     >  2:4847043/1‑150 (MQ=255)
                                                                                                                                tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGc                     >  1:1683146/1‑150 (MQ=255)
                                                                                                                                  tAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCa                   <  2:3797262/150‑1 (MQ=255)
                                                                                                                                   aaaaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAg                  >  2:1910610/1‑150 (MQ=255)
                                                                                                                                   aaaaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAg                  <  1:1759238/150‑1 (MQ=255)
                                                                                                                                     aaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGc                <  2:4454772/150‑1 (MQ=255)
                                                                                                                                         aaaCCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCAt            <  2:735695/150‑1 (MQ=255)
                                                                                                                                           aCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCCGGCCCATaa          >  1:1887259/1‑150 (MQ=255)
                                                                                                                                           aCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATaa          >  1:3175071/1‑150 (MQ=255)
                                                                                                                                              tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCa       >  1:1134088/1‑150 (MQ=255)
                                                                                                                                               ccGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCCGGCCCATAATCAg      >  1:1539393/1‑150 (MQ=255)
                                                                                                                                               ccGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAg      <  1:1005700/150‑1 (MQ=255)
                                                                                                                                                  cATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCCGGCCCATAATCAgacg   >  1:2547305/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCAGGCCGGCATCGCCGGACGCTCAAAAGAAATGTGTCCCTATCAGACGCTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAGACGACCATCAACCAGCACTTCGCCATCCGCCAGGCCCATAATCAGACGA  >  NC_000913/1015764‑1016060

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: