Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,020 T→A 6.6% R204* (AGA→TGA) ‡ dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,0200TA 6.6% 247.7 / 7.1 91R204* (AGA→TGA) ‡dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (21/62);  new base A (4/2);  total (27/64)
Fisher's exact test for biased strand distribution p-value = 5.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

ATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGT  >  NC_000913/2476872‑2477166
                                                                                                                                                    |                                                                                                                                                  
aTCGCCACCGCATTATGATTCATCGCGGCAATATCAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCCGCAAATTAACCGCATAATGTTGCGCCGATCTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtctg                                                                                                                                                   <  1:1925612/150‑1 (MQ=255)
  cGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTATCCGCATAATGTTGCGCTCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCtgtctgtc                                                                                                                                                 <  2:1867545/150‑1 (MQ=255)
   gCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTACAGCAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCa                                                                                                                                                >  2:1542292/1‑150 (MQ=255)
   accACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGCTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTACATGCcggtctggaa                                                                                                                                                >  1:26199/2‑147 (MQ=255)
      aCCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGCACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATgg                                                                                                                                             <  2:1118773/150‑1 (MQ=255)
      aCCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATgg                                                                                                                                             <  1:1057338/150‑1 (MQ=255)
       ccccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAg                                                                                                                                           <  2:717670/150‑1 (MQ=255)
        cacAATATGAATAATCGCGTCAATAACAGCTAAATCAGAACGATCAAAGCCATTTACATAATATGTCGGCAAATTAACCGCATAATGTTTCGCACAACTATGCCATTCATCCGTACCCTTGTCTTTGCTCCATGCCTGTCTGTCATGGAg                                                                                                                                           <  2:302918/148‑1 (MQ=255)
         ccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCATTGAGTCATGGAGc                                                                                                                                          >  2:390932/1‑150 (MQ=255)
         ccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGc                                                                                                                                          <  1:1550434/150‑1 (MQ=255)
            atatGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCAAc                                                                                                                                       <  2:1612994/150‑1 (MQ=255)
                gATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTAc                                                                                                                                   <  2:589614/150‑1 (MQ=255)
                gATTCATCGCCGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTAc                                                                                                                                   <  1:70162/150‑1 (MQ=255)
                 aTTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTACATGGCTGTGTGTAATGGAGCAAcattca                                                                                                                                  >  1:1480874/1‑145 (MQ=255)
                   tCATCGCGGCAATAACAGCTAAATCAGACCGATCAAAGCCAATTCCAGCAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGt                                                                                                                                <  1:1183495/150‑1 (MQ=255)
                   tCATCGCGGCAATAACAGCTAAATCAGAACGATCAAATCCACTTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGt                                                                                                                                <  1:1191301/150‑1 (MQ=255)
                   tCATCGAGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGt                                                                                                                                <  2:534926/150‑1 (MQ=255)
                     aTCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTGCCAGGCGGTGTGTGATGGAGCCACGGACACTGa                                                                                                                              >  2:182018/1‑150 (MQ=255)
                      tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGac                                                                                                                             >  2:360543/1‑150 (MQ=255)
                      tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCAGCATGCCTGTCTGTGATGGAGCAACGTACAggggc                                                                                                                             >  1:1286332/1‑146 (MQ=255)
                        gcgGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCAACAAGGCAGGGTGTCATGGGGCCAACTACAGTGacaa                                                                                                                           >  1:1177897/1‑149 (MQ=255)
                          ggCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTCTCCAAGCCTGTCTGTCATTGAGCAACGTACAGTGACACAg                                                                                                                         >  2:246159/1‑150 (MQ=255)
                            cAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGt                                                                                                                       <  2:2032894/150‑1 (MQ=255)
                                    gCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTAcc                                                                                                               <  1:400220/150‑1 (MQ=255)
                                           cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTgg                                                                                                        <  2:765022/150‑1 (MQ=255)
                                            aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTTATCCATGCCTGTATGTCAAGGAGCCACGGTCAGTGACACCGGGTTATTACCGTGTTTGt                                                                                                       >  1:28610/1‑150 (MQ=255)
                                            aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGt                                                                                                       <  2:1908569/150‑1 (MQ=255)
                                            aGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTACATGCCAGTGGGTGATGGAGCCACGGTCAGGGCCCCAGGGTAATTAACGTGTTGGt                                                                                                       >  1:1256778/1‑150 (MQ=255)
                                                cGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCAAGCCAGTGGGTCATCGAGCAACCTACCAGGAGACAGGGTAATTACCGGGTTGGGTaaa                                                                                                   >  2:398030/1‑150 (MQ=255)
                                                  aTCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGc                                                                                                 <  1:774168/150‑1 (MQ=255)
                                                   tCAAAGCCCATTCCAGAAGATGTCGGCAAATTAACCGCATACTGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCa                                                                                                <  2:1194536/150‑1 (MQ=255)
                                                   tCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCa                                                                                                <  2:771771/150‑1 (MQ=255)
                                                   tCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCa                                                                                                <  1:242369/150‑1 (MQ=255)
                                                     aaaGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCGGTCTGTCCTGGTGCCACGTACAGTGACACCGGTTTATTACCGTATTGGTTATAGCATg                                                                                              >  2:331289/1‑150 (MQ=255)
                                                            aTTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGa                                                                                       <  2:1476798/150‑1 (MQ=255)
                                                               ccAGAAGATTTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                    <  2:1004722/150‑1 (MQ=255)
                                                               ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                    <  1:705837/150‑1 (MQ=255)
                                                               ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                    <  2:185375/150‑1 (MQ=255)
                                                               acAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCg                                                                                    <  2:1499997/149‑1 (MQ=255)
                                                                  gaagaTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTAt                                                                                 <  2:1215380/150‑1 (MQ=255)
                                                                     gaTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCCTGGCTGTCTGTCATGGGGCAACGTACAATGACCCAGGTTTATTTACGTGTTGGTTTATTCATGTGTTTTAGCGTAttgc                                                                              >  1:24631/1‑148 (MQ=255)
                                                                           ggCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTa                                                                        <  1:840662/150‑1 (MQ=255)
                                                                            gCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCACCATGCCTGTCTGTCATGGAGCAACGTACACTGACACCGGTTAATTAACGTGTTGGTTAAAGCATTTGTTTGAGCGTATTCCGGGCTaa                                                                       >  2:1489369/1‑149 (MQ=255)
                                                                             cAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACa                                                                      <  2:1246568/150‑1 (MQ=255)
                                                                                 ttAACCGCATAATGTTGCGCCCATCTATGCCATTCATCCGTACCGGAGTCGTTGCTACATGCCTGTCTGTCATGGAGCAACGTACAGTGACACATGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACt                                                                  <  2:1267953/150‑1 (MQ=255)
                                                                                 ttAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACACACt                                                                  <  1:1348888/150‑1 (MQ=255)
                                                                                  tAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCGATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTg                                                                 <  1:1473562/150‑1 (MQ=255)
                                                                                   aaCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTgg                                                                <  2:1443508/150‑1 (MQ=255)
                                                                                    aCCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGc                                                               <  2:824499/150‑1 (MQ=255)
                                                                                       gCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCAGTCTGTCATGGAGCAAAGGACACTGACACAGGTTTATTACCGTGGTTGTTAAATCATGTCTTTGAGCCTATTCCCGGCTTAATACTGGAAgt                                                            >  1:363720/1‑149 (MQ=255)
                                                                                        cATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGa                                                           <  1:1005297/150‑1 (MQ=255)
                                                                                         aTAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCAGGCCTGTCTGTCATGGAGCAACGTACAGTGACACATGTTAATTACCGTGTTGGTTAAAGCATGTATTTGAGAGTATTACGGGGAACATACTGTCAGGAt                                                          >  2:344834/1‑150 (MQ=255)
                                                                                         aTAATGTTGCGACCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGGCAAGGAGCACCGAACAGTGACACAAGGTAATTACCGTGTTGGTTAAAGCATGTCTTTGGCCGAATTGCCGGCTACATACATGCAGGAt                                                          >  2:109044/1‑150 (MQ=255)
                                                                                          tAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAtt                                                         <  2:2124730/150‑1 (MQ=255)
                                                                                          tAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAtt                                                         <  1:1654697/150‑1 (MQ=255)
                                                                                           aaTGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAttt                                                        <  1:798679/150‑1 (MQ=255)
                                                                                           aaTGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGGCTGTCATGGAGCAACGTACAGTGACACAGGTTAAATACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAttt                                                        <  2:185845/150‑1 (MQ=255)
                                                                                             tGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttct                                                      <  1:2067779/150‑1 (MQ=255)
                                                                                               ttGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCAAGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttcttc                                                    <  2:1901897/150‑1 (MQ=255)
                                                                                                tGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTca                                                   <  1:656871/150‑1 (MQ=255)
                                                                                                   gcCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatcc                                                <  1:953252/150‑1 (MQ=255)
                                                                                                    cccAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatcca                                               <  2:599461/150‑1 (MQ=255)
                                                                                                      cAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTTACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatccatc                                             <  1:1461870/150‑1 (MQ=255)
                                                                                                      cAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatccatc                                             <  1:782428/150‑1 (MQ=255)
                                                                                                       aaCTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCa                                            <  1:360543/150‑1 (MQ=255)
                                                                                                              ccATTCATCCGTACCTGAGTAGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGGTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTTCATACTGGCAGGATTTCTTCATCCATCAGTAAGtg                                     <  2:653706/150‑1 (MQ=255)
                                                                                                               cATTCATCCGTACCGGAGTCGTGGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACTCAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtga                                    <  2:643991/150‑1 (MQ=255)
                                                                                                                    aTCCGTCCCGGCGTCGGTGCTCCAAGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGGTAATTACCGTGTTGGTTAAAGCATCTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagt                               <  2:173562/150‑1 (MQ=255)
                                                                                                                     tCCGTACCGGAGTCGTTTCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                              <  2:1534717/150‑1 (MQ=255)
                                                                                                                     tCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtg                              <  2:363720/150‑1 (MQ=255)
                                                                                                                      ccGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTTCAGTGACACAGGTTAATTTCCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtga                             <  1:621453/150‑1 (MQ=255)
                                                                                                                      ccGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtga                             <  1:811487/150‑1 (MQ=255)
                                                                                                                          aCCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTg                         <  1:729593/150‑1 (MQ=255)
                                                                                                                          aCCGGAGTCGTAGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTg                         <  1:390932/150‑1 (MQ=255)
                                                                                                                            cGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGAAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgcg                       <  2:1579202/150‑1 (MQ=255)
                                                                                                                              gAGTCGTTGCTCCATGCCTGGCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAGTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCt                     <  2:1488382/150‑1 (MQ=255)
                                                                                                                               aGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTg                    <  1:1264069/150‑1 (MQ=255)
                                                                                                                                gTCGTTGCTCCATGCCTGTCTTGCATGGAGCAACGTACAGAGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGa                   <  1:873356/150‑1 (MQ=255)
                                                                                                                                   gTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg                >  2:590660/1‑150 (MQ=255)
                                                                                                                                    ttGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgc               <  2:158214/150‑1 (MQ=255)
                                                                                                                                    ttGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgc               <  1:1221937/150‑1 (MQ=255)
                                                                                                                                      gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc             >  1:1606470/1‑150 (MQ=255)
                                                                                                                                      gCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc             >  1:1034145/1‑150 (MQ=255)
                                                                                                                                       cTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca            >  2:1287508/1‑150 (MQ=255)
                                                                                                                                       cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca            >  1:1571954/1‑150 (MQ=255)
                                                                                                                                       cTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGATGGTTAAAGCATGTCTATGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGca            >  1:1466640/1‑150 (MQ=255)
                                                                                                                                         ccATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc          <  1:761599/150‑1 (MQ=255)
                                                                                                                                          cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaact         <  1:1757985/150‑4 (MQ=255)
                                                                                                                                           aTGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCCGTCAGCGcaactc        <  1:901383/150‑5 (MQ=255)
                                                                                                                                            tGCCTGTCTGTCCTGGAGCAACGTACAGTGACACAGGTTAATTACCGCGGTGGTTAAAGCATGTCTTTGAGCGGATTCCGGGCTACATACTGGCCGGATTTCTTCATCCAGCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc       <  1:1211691/150‑1 (MQ=255)
                                                                                                                                            tGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc       >  1:671356/1‑150 (MQ=255)
                                                                                                                                             gCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCAtcatca      >  2:1922748/1‑150 (MQ=255)
                                                                                                                                              ccTGTCTGTCATGGAGCAACGTACAGTTACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCATTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa     >  2:1686105/1‑150 (MQ=255)
                                                                                                                                              ccTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaa     >  2:614135/1‑150 (MQ=255)
                                                                                                                                               ctgtctgtCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaaa    <  2:855531/150‑1 (MQ=255)
                                                                                                                                               ctcactgtcATGGAGAAACGTACAGTGACAAATGTTAATTACCCTGTTGGTTAAAGCATATCTTTGAGCGTATTCCTAGATACATACTGCCACGATTTCATAATACAACAGAAAGTGATGAGTGAGTTGCGCTGACCACGCAACATaaaa    >  2:944158/5‑150 (MQ=255)
                                                                                                                                                tgtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAg   <  1:2085553/150‑1 (MQ=255)
                                                                                                                                                tgtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAg   <  2:1177897/150‑1 (MQ=255)
                                                                                                                                                 gtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGt  <  1:1933558/150‑1 (MQ=255)
                                                                                                                                                 gtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGt  >  2:1520774/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGT  >  NC_000913/2476872‑2477166

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: