Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,647,670 T→C 5.4% intergenic (+19/+35) gor → / ← dinQ glutathione oxidoreductase/UV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,647,6700TC 5.4% 276.6 / 4.8 93intergenic (+19/+35)gor/dinQglutathione oxidoreductase/UV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (45/43);  new base C (2/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATG  >  NC_000913/3647522‑3647818
                                                                                                                                                    |                                                                                                                                                    
tGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCTACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGagta                                                                                                                                                     >  2:50685/1‑150 (MQ=255)
 ggCTTTGGTATGGACGAAATGTTGCAGTGCTTCGCGGTGGCGCTGAAGATGGGGGCTACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGagtag                                                                                                                                                    <  2:1636677/150‑1 (MQ=255)
 ggCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGagtag                                                                                                                                                    >  2:1017196/1‑150 (MQ=255)
  gCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGagtagt                                                                                                                                                   >  1:438927/1‑150 (MQ=255)
       ggTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGct                                                                                                                                              <  1:56616/150‑1 (MQ=255)
        gTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGctc                                                                                                                                             >  2:1233570/1‑150 (MQ=255)
        gTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGctc                                                                                                                                             <  1:189099/150‑1 (MQ=255)
             gACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGc                                                                                                                                        <  1:1230883/150‑1 (MQ=255)
             gACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGc                                                                                                                                        >  1:534000/1‑150 (MQ=255)
              aCGAATTGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCATAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGCACACTACTCTTAGcc                                                                                                                                       <  1:305721/150‑1 (MQ=255)
              aCGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGcc                                                                                                                                       <  1:1639974/150‑1 (MQ=255)
               cGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGccc                                                                                                                                      >  2:194479/1‑150 (MQ=255)
               cGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGccc                                                                                                                                      <  1:227938/150‑1 (MQ=255)
               cGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGCACACTACTCTTAGccc                                                                                                                                      <  1:1451311/150‑1 (MQ=255)
                gAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGCACACTACTCTTAGCCCt                                                                                                                                     <  2:1185346/150‑1 (MQ=255)
                  aaTGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTa                                                                                                                                   <  1:1441076/150‑1 (MQ=255)
                  aaTGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTa                                                                                                                                   <  1:894030/150‑1 (MQ=255)
                        gCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACg                                                                                                                             <  2:61377/150‑1 (MQ=255)
                        gCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGAATTAGACAATACCGTCGCCATTCAGGCAACGGCGGCAGAAGAGTACGTGACAATGGGGTAAATGTTAAAGGGGAAAGAGTACTGTGCTCAAAGCCATTGAATAct                                                                                                                             >  2:1102831/1‑149 (MQ=255)
                            ggCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTc                                                                                                                         <  2:188571/150‑1 (MQ=255)
                               ttCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGc                                                                                                                      <  1:524298/150‑1 (MQ=255)
                               ttCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGc                                                                                                                      <  2:1580238/150‑1 (MQ=255)
                                 cgcgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACGATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTa                                                                                                                    >  1:1704521/1‑150 (MQ=255)
                                  gcgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTAt                                                                                                                   <  1:471658/150‑1 (MQ=255)
                                   cgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACATTGGGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATc                                                                                                                  >  1:1256975/1‑150 (MQ=255)
                                   cgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATc                                                                                                                  <  1:1494066/150‑1 (MQ=255)
                                   cgGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATc                                                                                                                  <  1:961240/150‑1 (MQ=255)
                                      tGGCGCTGAAGATGGGGGCAACCAAAAAATACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTTACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGt                                                                                                               <  2:776181/150‑1 (MQ=255)
                                        gcgcTGAAGATGGGTGAAACCAAAAAAGTCTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGAGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTc                                                                                                             <  1:326372/150‑1 (MQ=255)
                                         cgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCa                                                                                                            >  2:1167398/1‑150 (MQ=255)
                                         cgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCa                                                                                                            <  1:375247/150‑1 (MQ=255)
                                         cgcTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCa                                                                                                            >  2:1196084/1‑150 (MQ=255)
                                           cTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGa                                                                                                          >  2:1486433/1‑150 (MQ=255)
                                           cTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTAGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTACGCTATAAGTTCAGa                                                                                                          >  1:877325/1‑150 (MQ=255)
                                             gAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAg                                                                                                        >  2:1340263/1‑150 (MQ=255)
                                               aGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCt                                                                                                      <  1:1017196/150‑1 (MQ=255)
                                               aGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCt                                                                                                      >  2:434982/1‑150 (MQ=255)
                                                gATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTg                                                                                                     >  2:275497/1‑150 (MQ=255)
                                                 aTGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGa                                                                                                    <  1:793115/150‑1 (MQ=255)
                                                    ggggCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACCATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGc                                                                                                 >  1:240220/1‑150 (MQ=255)
                                                      ggCAACCAAAAAAGAATTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGATAAATGTTAACGGGCTAAGAGTAGTGTGCTCTTAGCCCGTAATTACGTTTCCGCTATCAGTTAAGAAGCTgcagcac                                                                                               >  1:1273734/1‑149 (MQ=255)
                                                          aCCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAg                                                                                           >  2:156514/1‑150 (MQ=255)
                                                          aCCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAg                                                                                           >  1:751584/1‑150 (MQ=255)
                                                                aaaGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATc                                                                                     >  2:161878/1‑150 (MQ=255)
                                                                  aGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGGTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAg                                                                                   >  2:1443814/1‑150 (MQ=255)
                                                                  aGACATCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAg                                                                                   <  2:1116310/150‑1 (MQ=255)
                                                                     cTTCGACAATACCTTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTc                                                                                <  1:9932/150‑1 (MQ=255)
                                                                     cTTCGACAATACCGTCGCCTTTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTc                                                                                >  2:1667076/1‑150 (MQ=255)
                                                                     cTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTc                                                                                <  1:631646/150‑1 (MQ=255)
                                                                       tCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCca                                                                              <  1:275497/150‑1 (MQ=255)
                                                                          aaaaTACCGTCGCCATTCACCCAACGGCGGCATAAGAGTTCGTTACAATGCGTTAAATGTTAACGGGCTAAGCGTCGTGGGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCCGAAGCTGAAGCAGAAAGCGTTCAGTTCcagcag                                                                          <  2:1234385/148‑1 (MQ=255)
                                                                           cAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcagcag                                                                          <  2:248445/150‑1 (MQ=255)
                                                                           aaaTACCGTCGCCATTCACCCAACGGCGTCAGAAGAGTTCGTGACAATGCGTTAACTGTTAAAGGGCTAAGAGTAGTGGGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCcagcag                                                                          <  1:871381/149‑1 (MQ=255)
                                                                               aCCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCa                                                                      <  1:615212/150‑1 (MQ=255)
                                                                                ccGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCaa                                                                     <  1:1671160/150‑1 (MQ=255)
                                                                                ccGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCaa                                                                     >  2:1529550/1‑150 (MQ=255)
                                                                                ccGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCaa                                                                     >  1:124432/1‑150 (MQ=255)
                                                                                  gTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTcc                                                                               <  2:522366/138‑1 (MQ=255)
                                                                                  gTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTcc                                                                               >  1:522366/1‑138 (MQ=255)
                                                                                     gCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAACCAGAAAGCGGAACAGTTCCAGCAGCGCAATTAAc                                                                >  1:1251022/1‑150 (MQ=255)
                                                                                       cATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGc                                                              >  1:1689467/1‑150 (MQ=255)
                                                                                       cATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGc                                                              <  1:194479/150‑1 (MQ=255)
                                                                                         ttCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGccc                                                            >  1:1615853/1‑150 (MQ=255)
                                                                                            aCCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTa                                                         >  1:39872/1‑150 (MQ=255)
                                                                                              ccAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGa                                                       <  1:397621/150‑1 (MQ=255)
                                                                                              ccAACGGCGCCAGAAGAGTTCGTGACAATGCGTTAAATGTTAACGGGCTAAGCCTAGTCAGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAAGTAACGCCCCTAGa                                                       >  1:78480/1‑150 (MQ=255)
                                                                                                aaCGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAAc                                                     >  2:928550/1‑150 (MQ=255)
                                                                                                    gcttcAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGCGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTGCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACgatg                                                 <  2:72233/146‑1 (MQ=255)
                                                                                                      ggcAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACgatgat                                               >  2:1305703/1‑150 (MQ=255)
                                                                                                       gcAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGAtt                                              >  1:884354/1‑150 (MQ=255)
                                                                                                          gaagaGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCt                                           >  1:1182446/1‑150 (MQ=255)
                                                                                                          gaagaGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCt                                           >  1:1438615/1‑150 (MQ=255)
                                                                                                            agagTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTACCGCCCCTAGAACGATGATTGCttt                                         >  2:1119362/1‑150 (MQ=255)
                                                                                                               tttCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatc                                      <  1:968093/149‑1 (MQ=255)
                                                                                                                  cGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTatcaat                                   >  2:1188629/1‑150 (MQ=255)
                                                                                                                     gACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCTCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCAc                                <  2:438927/150‑1 (MQ=255)
                                                                                                                     gACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCAc                                >  2:1315543/1‑150 (MQ=255)
                                                                                                                     gACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCAc                                >  1:138144/1‑150 (MQ=255)
                                                                                                                        aaTGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCg                             >  2:475784/1‑150 (MQ=255)
                                                                                                                        aaTGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCg                             <  2:216294/150‑1 (MQ=255)
                                                                                                                              ttAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTtct                       <  2:1467501/150‑1 (MQ=255)
                                                                                                                               tAAATGTTAAAGGGCTAAGAGTAGAGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCATAGAACGATGATTGCTTAATCAATCACCCGTTTtctc                      >  2:953223/1‑150 (MQ=255)
                                                                                                                                aaaTGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTcc                     <  1:434982/150‑1 (MQ=255)
                                                                                                                                      tAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGa               >  1:574301/1‑150 (MQ=255)
                                                                                                                                       aaaGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAATGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGAt              <  1:156514/150‑1 (MQ=255)
                                                                                                                                       aaaGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGAt              <  1:1486433/150‑1 (MQ=255)
                                                                                                                                         aGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg            <  2:1438615/150‑1 (MQ=255)
                                                                                                                                         aGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg            >  2:1483564/1‑150 (MQ=255)
                                                                                                                                         aGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg            <  2:1182446/150‑1 (MQ=255)
                                                                                                                                         aGGGCTAAGAGCACACTACTCTTAGCCCTTAATTACGTTTCGGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg            >  2:661961/1‑150 (MQ=255)
                                                                                                                                         aGGGCTAAGAGCACACTACTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATgg            >  1:201363/1‑150 (MQ=255)
                                                                                                                                           ggCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAg          <  1:1305703/150‑1 (MQ=255)
                                                                                                                                           ggCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAg          <  2:124432/150‑1 (MQ=255)
                                                                                                                                           ggCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAg          >  2:1699055/1‑150 (MQ=255)
                                                                                                                                            gCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGt         >  1:1013900/1‑150 (MQ=255)
                                                                                                                                             cTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTATCCATGCGATGGAGTg        >  1:474931/1‑150 (MQ=255)
                                                                                                                                               aaGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATAAATCACCCGTTTTCTCCATGCGATGGAGTgag      >  2:1661114/1‑150 (MQ=255)
                                                                                                                                                  agtagtGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCNNAAGCTGANGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAAt   >  1:1702994/1‑150 (MQ=255)
                                                                                                                                                  agtagtGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAAt   >  1:1709507/1‑150 (MQ=255)
                                                                                                                                                   gtagtGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCATCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATg  <  2:1394548/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGAGTAGTGTGCTCTTAGCCCTTAATTACGTTTCCGCTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATG  >  NC_000913/3647522‑3647818

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