Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,470,499 T→A 4.2% R204* (AGA→TGA)  dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,470,4990TA 4.2% 271.0 / 2.6 88R204* (AGA→TGA) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (40/43);  new base A (2/2);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

ATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTA  >  CP009273/2470351‑2470646
                                                                                                                                                    |                                                                                                                                                   
aTCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCctgtct                                                                                                                                                     <  2:593166/149‑1 (MQ=255)
    ccACCCCAATATGATTCATCGCGGGAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGAAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCa                                                                                                                                                 >  2:649734/1‑149 (MQ=255)
    ccACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCa                                                                                                                                                 >  2:268487/1‑149 (MQ=255)
    ccACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCa                                                                                                                                                 <  2:156802/149‑1 (MQ=255)
     cACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCAt                                                                                                                                                <  2:27209/149‑1 (MQ=255)
       ccccATTATGCTTCATCGCGGCAATAACAGCTAAATCAGAGCGATCAATGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCTACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATgg                                                                                                                                              <  2:566127/149‑1 (MQ=255)
       ccccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATgg                                                                                                                                              <  2:565896/149‑1 (MQ=255)
         ccAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCACCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTGATGGAg                                                                                                                                            >  2:46205/1‑149 (MQ=255)
           aaTATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCa                                                                                                                                          <  2:225679/149‑1 (MQ=255)
            atatGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACAGACAGGCATGGAGCaa                                                                                                                                         <  2:402296/149‑1 (MQ=255)
                  ttCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACa                                                                                                                                   <  2:702858/149‑1 (MQ=255)
                   tCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCCACGTACAg                                                                                                                                  >  1:119962/1‑149 (MQ=255)
                   tCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAg                                                                                                                                  <  2:549034/149‑1 (MQ=255)
                   tCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAg                                                                                                                                  >  2:445564/1‑149 (MQ=255)
                   tCATCGCGGCAATAACAGCTAAATAAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAg                                                                                                                                  <  2:431441/149‑1 (MQ=255)
                     atgccgGCAACAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTg                                                                                                                                <  1:520976/145‑1 (MQ=255)
                     aTCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTg                                                                                                                                >  2:307949/1‑149 (MQ=255)
                     aTCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATGAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGACTGTCAGTCAATGAGCAACGTACAggg                                                                                                                                >  2:51741/1‑147 (MQ=255)
                      tCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGa                                                                                                                               <  1:743387/149‑1 (MQ=255)
                        gcgGCAAGAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACGATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGaca                                                                                                                             <  2:831049/149‑1 (MQ=255)
                             aaTAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGt                                                                                                                        <  1:205880/149‑1 (MQ=255)
                                aaCAGGTAAAGCACAACCACCAAAGCCAATTCCAAAAGATGTCGGCAAATCAACCGCATAATGGTGCGCCCCACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTaa                                                                                                                     <  1:90534/149‑1 (MQ=255)
                                aaCAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTaa                                                                                                                     <  2:366618/149‑1 (MQ=255)
                                 aCAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAAt                                                                                                                    >  2:372767/1‑149 (MQ=255)
                                  cAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAAtt                                                                                                                   <  2:818031/149‑1 (MQ=255)
                                  cAGCTAAAGCAGCACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAAtt                                                                                                                   <  1:627035/149‑1 (MQ=255)
                                   aGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTa                                                                                                                  >  1:683109/1‑149 (MQ=255)
                                    gCTAAATCACAACGCTCGACGCCAATGCCAGTAGAGGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTAc                                                                                                                 <  2:236172/149‑1 (MQ=255)
                                       aaaTCAGAACGATCAAAGCCAATGCCAGAAGATGTCGGCAAATTGACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCgt                                                                                                              <  1:250405/149‑1 (MQ=255)
                                          tCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGGGCAACGTACAGGGACACAGGTTAATTTCCGTGtt                                                                                                           >  2:813917/1‑149 (MQ=255)
                                          tCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGtt                                                                                                           >  2:868080/1‑149 (MQ=255)
                                           cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTg                                                                                                          <  2:513900/149‑1 (MQ=255)
                                           cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTg                                                                                                          <  2:621152/149‑1 (MQ=255)
                                           cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTg                                                                                                          <  2:81945/149‑1 (MQ=255)
                                           cAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTg                                                                                                          >  2:311330/1‑149 (MQ=255)
                                           aaGAACGATAAAAGCCAATTCTAGATGATGTCGGCACAGCAACCGCATACTGTTGCGCCCAACTATTCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTAAATTACCGTGTTg                                                                                                          <  1:27173/148‑1 (MQ=255)
                                             gAACTATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGt                                                                                                        <  1:865704/149‑1 (MQ=255)
                                              aaCGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAAACGCATAATGTTACGCCAAACTATGCCATTCATCCGTCCCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGtt                                                                                                       >  1:207459/1‑149 (MQ=255)
                                                cGACCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTaa                                                                                                     <  2:663727/149‑1 (MQ=255)
                                                   tgaaaTCCAATTCCAGAATATGTCGCCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATTCCTGTCTGTCATGGAGCAACGTACAGTGACACATGTTAATTACCGTGTTGGTTAAAGc                                                                                                  <  1:743810/147‑1 (MQ=255)
                                                              tCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAg                                                                                       >  2:360785/1‑149 (MQ=255)
                                                              tCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGGGTTGGTTAAAGCATGTCTTTGAg                                                                                       >  1:348181/1‑149 (MQ=255)
                                                               ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGc                                                                                      >  1:644514/1‑149 (MQ=255)
                                                               ccAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGc                                                                                      >  1:731988/1‑149 (MQ=255)
                                                                 agaagaTGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGt                                                                                    >  1:264895/1‑149 (MQ=255)
                                                                         tCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCAGCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCggg                                                                            <  1:302688/149‑1 (MQ=255)
                                                                           ggCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCATTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACCCAGGTTATGTTCCGTGGTGGTTAAAGCATGTGTTTGAGCGTATTCCGGGCt                                                                          >  2:285119/1‑149 (MQ=255)
                                                                             cAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTAc                                                                        >  1:761362/1‑149 (MQ=255)
                                                                             cAAATTAACCGCATAATGTTGCGCCCAACTATACCATTCATCCGTACCGGAGTCGTTGCTCCCTGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTACTTACCGTGTAGGTTAAAGCATGTCTTTGGGCGTATTCCGGGCtcc                                                                        >  1:346473/1‑147 (MQ=255)
                                                                              aaaTTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCGACa                                                                       >  1:799797/1‑149 (MQ=255)
                                                                                  tAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCGTGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACt                                                                   >  1:743930/1‑149 (MQ=255)
                                                                                  tAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCCTGGAGCAACGTACAGTGACCCAGGTTAATTACCGTGTTGGTTAAAGCATGACTTTGAGCGTATTCCGGGCTACATACt                                                                   >  2:339659/1‑149 (MQ=255)
                                                                                  tAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGAAGCAACGTACAGTGACACAGGTTAATTCCCGTGTTGGTTAAAGCATGTCTTTGCGCGTATTACGGGCTACATACt                                                                   >  1:299262/1‑149 (MQ=255)
                                                                                  tAACCGCATAATGTTGCGACCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGGCTTTGAGCGTATTCCGGGCTACATACt                                                                   >  1:560069/1‑149 (MQ=255)
                                                                                     ccGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGc                                                                >  1:828054/1‑149 (MQ=255)
                                                                                     ccGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCAGGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGAATTTGAGCGTATTCCGGGCTACATACTgag                                                                >  1:181400/1‑147 (MQ=255)
                                                                                     ccGCAGAATGTTGCGCCAAACTATGACATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGc                                                                >  1:181326/1‑149 (MQ=255)
                                                                                       gCATAATGTTGCGCCCAACTATGCCATTCATCCGTTCCGGAGTCGTTGCCCCATACCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTAAATTACCGTGATGGTTAAAGCGTGTCTTTGAGCGGATTCCGGACTACATACGGGCAg                                                              >  2:383244/1‑149 (MQ=255)
                                                                                       gCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTTTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAg                                                              <  2:296713/149‑1 (MQ=255)
                                                                                         aTAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGa                                                            >  1:769510/1‑149 (MQ=255)
                                                                                           aaTGTTGCGCCCAACTATGCCATTCATCCGTACCGGGGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGAtt                                                          <  1:82406/149‑1 (MQ=255)
                                                                                             cgtcgcgcCCATCTATGCCATTCATCCCTACCGCAGTCTTTGCTCCATGCCTATATGTATTTGACCATCGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATttc                                                        <  2:417635/145‑1 (MQ=255)
                                                                                                   gcCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGACATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatc                                                  >  2:816736/1‑149 (MQ=255)
                                                                                                     ccAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTcatcca                                                <  2:822088/149‑1 (MQ=255)
                                                                                                        agtctGCCATTCATCCCTACCTGAGTCGTTGCTCCATGCCTGTCTGTTATGGAAGAACGTACAGTGACACAAGTTAATTACCGTGTTGGTTAAAGAATGTCTTTCAGCGTTTTCCGGGCTACATACTGACAGGATTTCTTCATCCATCa                                             <  1:72717/145‑1 (MQ=255)
                                                                                                            tGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGaaa                                         <  2:273033/149‑1 (MQ=255)
                                                                                                            tGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGaaa                                         >  1:195419/1‑149 (MQ=255)
                                                                                                            tGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGaaa                                         >  1:465270/1‑149 (MQ=255)
                                                                                                            tGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTAACGTGTTGGTTAAAGCATGTCTTTGTGCGTATTCCGGGCTCCATACTGGCAAGATTTCTTCATCCTTCAGaaa                                         >  1:247551/1‑149 (MQ=255)
                                                                                                              ccATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGt                                       <  1:311330/149‑1 (MQ=255)
                                                                                                               cATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAACGACTTTTTGTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGtg                                      >  2:202572/1‑149 (MQ=255)
                                                                                                                    cctgtgtcactgtacGTTGCTCCATGACAGACAGGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgag                                 >  1:280774/15‑149 (MQ=25)
                                                                                                                       cGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGACAGTGAtgagtga                              >  1:627441/1‑149 (MQ=255)
                                                                                                                        gTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtgag                             <  1:307949/149‑1 (MQ=255)
                                                                                                                        gTACCGGAGTCGTTCGTCCCTGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGAtgagtgag                             <  1:51741/149‑1 (MQ=255)
                                                                                                                         tACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATgaaggagt                            >  1:177744/1‑149 (MQ=255)
                                                                                                                            cGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgc                         >  2:563329/1‑149 (MQ=255)
                                                                                                                             ggcgTCGTTGCTCCATGCCTGTCCGGCAGGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACACACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgcg                        <  1:649734/146‑1 (MQ=255)
                                                                                                                             ggAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgcg                        <  1:268487/149‑1 (MQ=255)
                                                                                                                              gAGTCTGTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTgcgc                       <  1:64930/149‑1 (MQ=255)
                                                                                                                                    ttGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACgg                 >  1:725040/1‑149 (MQ=255)
                                                                                                                                      gCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcg               <  2:352709/149‑1 (MQ=255)
                                                                                                                                       caccATGCCTGTCTGTCATGGAGCAGCGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc              <  1:762050/147‑1 (MQ=255)
                                                                                                                                       cTCCATGCCCGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGGGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGgcgc              <  1:808980/149‑1 (MQ=255)
                                                                                                                                          cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACGCGTTAATTGCCGTGGTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTCCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc           <  1:142392/149‑1 (MQ=255)
                                                                                                                                          cATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatc           <  2:617882/149‑1 (MQ=255)
                                                                                                                                           aTGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGcatca          >  2:360385/1‑149 (MQ=255)
                                                                                                                                             gCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc        <  1:680557/149‑1 (MQ=255)
                                                                                                                                             gCCTGTCTCGCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCatcatc        <  2:117630/149‑1 (MQ=255)
                                                                                                                                               ctgtctgtCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTCACGCGTGTCTATGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCCGAAAGTGATGAGTGAGTTGCGCTGCCGGCGCATCATCaa      >  1:653213/1‑149 (MQ=255)
                                                                                                                                                tgtctgtcATGGAGCAGCGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaaa     <  1:342739/149‑1 (MQ=255)
                                                                                                                                                tgtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaaa     <  1:540761/149‑1 (MQ=255)
                                                                                                                                                tgtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCaaa     <  2:242546/149‑1 (MQ=255)
                                                                                                                                                 gtctgtcATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAg    <  2:683109/149‑1 (MQ=255)
                                                                                                                                                   ctgtcATGGAGCCACGTATAGTTACTCCGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGTTACATACTGGCAGGATTTCTTCATCCATCAGAAATTGTGTAGTTAGTTGCGCTGACGGCGCATCATCAAAGTa  <  1:403089/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTCGGCAAATTAACCGCATAATGTTGCGCCCAACTATGCCATTCATCCGTACCGGAGTCGTTGCTCCATGCCTGTCTGTCATGGAGCAACGTACAGTGACACAGGTTAATTACCGTGTTGGTTAAAGCATGTCTTTGAGCGTATTCCGGGCTACATACTGGCAGGATTTCTTCATCCATCAGAAAGTGATGAGTGAGTTGCGCTGACGGCGCATCATCAAAGTA  >  CP009273/2470351‑2470646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑