Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,980,472 C→A 18.0% L23I (CTC→ATC)  ftnB → ferritin B, putative ferrous iron reservoir

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,980,4720CA18.0% 733.7 / 99.9 238L23I (CTC→ATC) ftnBferritin B, putative ferrous iron reservoir
Reads supporting (aligned to +/- strand):  ref base C (95/100);  new base A (19/24);  total (114/124)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 247 total aligned reads displayed.

AAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCC  >  CP009273/1980324‑1980621
                                                                                                                                                    |                                                                                                                                                     
aaaGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtct                                                                                                                                                      >  2:915250/1‑150 (MQ=255)
aaaGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtct                                                                                                                                                      >  2:2208167/1‑150 (MQ=255)
 aaGGAATTATTCTGCAAACAGTACGTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGACTGCTTCTCAAACTCAACTCTCAAAAGAACCGCGAGTTTTACGCATCCAAtatc                                                                                                                                                     >  2:524418/1‑150 (MQ=255)
   ggAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTa                                                                                                                                                   <  2:3964342/150‑1 (MQ=255)
    gAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATTAACCGCGAGTTTTACGCATCCAATCTCTAc                                                                                                                                                  >  1:2157572/1‑150 (MQ=255)
    gAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTAc                                                                                                                                                  >  2:4583924/1‑150 (MQ=255)
     aaTTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTAcc                                                                                                                                                 >  2:1290048/1‑150 (MQ=255)
        tattCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTc                                                                                                                                              <  1:106678/150‑1 (MQ=255)
         attCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCa                                                                                                                                             <  2:4137234/150‑1 (MQ=255)
         attCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCa                                                                                                                                             <  1:1254120/150‑1 (MQ=255)
             tGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                                                                                                                                         >  1:4012738/1‑150 (MQ=255)
             tGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTg                                                                                                                                         <  1:2430424/150‑1 (MQ=255)
             tGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCACTCCACTTCAGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCAACCTTCACCTg                                                                                                                                         >  2:1715255/1‑150 (MQ=255)
               cAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAg                                                                                                                                       <  2:3849779/150‑1 (MQ=255)
                aaaCAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGt                                                                                                                                      >  1:672439/1‑150 (MQ=255)
                aaaCAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGt                                                                                                                                      >  1:2871375/1‑150 (MQ=255)
                  aCAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGTTATTTTGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCATATATCTACCTTCACCTGAGTaa                                                                                                                                    <  2:1944630/150‑1 (MQ=255)
                  aCAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTaa                                                                                                                                    <  2:2558236/150‑1 (MQ=255)
                     gTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTg                                                                                                                                 <  2:4230491/150‑1 (MQ=255)
                     gTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTg                                                                                                                                 >  2:4344168/1‑150 (MQ=255)
                      tAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg                                                                                                                                >  2:4235038/1‑150 (MQ=255)
                             ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTg                                                                                                                         <  1:1919651/150‑1 (MQ=255)
                             ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTg                                                                                                                         >  1:1300673/1‑150 (MQ=255)
                              gtgtTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGa                                                                                                                        >  1:1076033/1‑150 (MQ=255)
                                gtTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCAAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAAc                                                                                                                      <  1:3617897/150‑1 (MQ=255)
                                    tGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACagag                                                                                                                  <  2:710834/150‑1 (MQ=255)
                                      aTTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGTCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTc                                                                                                                <  2:3509691/150‑1 (MQ=255)
                                        ttATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTg                                                                                                              >  1:859808/1‑150 (MQ=255)
                                            cTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACg                                                                                                          <  1:2086377/150‑1 (MQ=255)
                                              tGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGc                                                                                                        >  1:4402851/1‑150 (MQ=255)
                                              tGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGc                                                                                                        <  1:3463355/150‑1 (MQ=255)
                                                 aCCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCAcc                                                                                                     <  1:3180949/150‑1 (MQ=255)
                                                 aCCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCAcc                                                                                                     <  2:1786884/150‑1 (MQ=255)
                                                    tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGcc                                                                                                  >  1:478677/1‑150 (MQ=255)
                                                    tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGcc                                                                                                  >  1:216411/1‑150 (MQ=255)
                                                        cACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACtt                                                                                              >  2:1657499/1‑150 (MQ=255)
                                                        cACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACtt                                                                                              >  2:176100/1‑150 (MQ=255)
                                                        cACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACtt                                                                                              >  2:4290676/1‑150 (MQ=255)
                                                           ttCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTcc                                                                                           <  1:294123/150‑1 (MQ=255)
                                                            tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCt                                                                                          >  2:4119888/1‑150 (MQ=255)
                                                              tGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTc                                                                                        <  1:1592171/150‑1 (MQ=255)
                                                              tGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTc                                                                                        <  2:1733548/150‑1 (MQ=255)
                                                              tGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTc                                                                                        <  2:3594915/150‑1 (MQ=255)
                                                              tGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTc                                                                                        >  1:2231007/1‑150 (MQ=255)
                                                              tGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTc                                                                                        <  1:2569590/150‑1 (MQ=255)
                                                                gATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTcgc                                                                                      <  2:28964/150‑1 (MQ=255)
                                                                gATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTcgc                                                                                      >  2:1377725/1‑150 (MQ=255)
                                                                  tataAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCgcgc                                                                                    >  1:1525723/1‑150 (MQ=255)
                                                                      aGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAg                                                                                <  1:1378929/150‑1 (MQ=255)
                                                                        gATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGc                                                                              <  2:1184938/150‑1 (MQ=255)
                                                                         attttAGGTATGGCAACAGATGGAATGATTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACATTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTACTTCGCGCCCAGGca                                                                             <  1:1900585/147‑1 (MQ=255)
                                                                         atatTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAANGGCACCGCCACTTTCCTTCGCGCCCAGGca                                                                             <  1:4285506/150‑1 (MQ=255)
                                                                         atatTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCGACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGca                                                                             <  2:1592145/150‑1 (MQ=255)
                                                                           atTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGcaca                                                                           <  2:2437535/150‑1 (MQ=255)
                                                                           atTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGcaca                                                                           <  1:3020874/150‑1 (MQ=255)
                                                                           atTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGcaca                                                                           <  1:2872907/150‑1 (MQ=255)
                                                                            ttAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTATGAACATAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGTCACag                                                                          >  1:1693078/1‑150 (MQ=255)
                                                                             tAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACaga                                                                         >  2:3305382/1‑150 (MQ=255)
                                                                             tAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACaga                                                                         <  2:4012738/150‑1 (MQ=255)
                                                                              aGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACagag                                                                        >  1:2744931/1‑150 (MQ=255)
                                                                               ggTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGt                                                                       <  2:2921026/150‑1 (MQ=255)
                                                                               ggTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGt                                                                       <  2:3672774/150‑1 (MQ=255)
                                                                                gTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTa                                                                      <  1:3317819/150‑1 (MQ=255)
                                                                                gTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTa                                                                      <  2:4012457/150‑1 (MQ=255)
                                                                                   tGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAAtg                                                                   <  1:943660/150‑1 (MQ=255)
                                                                                     gCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAAtgtg                                                                 <  2:2648692/150‑1 (MQ=255)
                                                                                             tGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAat                                                         >  1:4096545/1‑150 (MQ=255)
                                                                                                 aTGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAatgatg                                                     <  1:2034593/150‑1 (MQ=255)
                                                                                                  tGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                    <  1:3864459/150‑1 (MQ=255)
                                                                                                   gCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgcg                                                   <  2:1300673/150‑1 (MQ=255)
                                                                                                   gCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgcg                                                   <  1:327005/150‑1 (MQ=255)
                                                                                                   gCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgcg                                                   <  1:168151/150‑1 (MQ=255)
                                                                                                       ctcAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATg                                               <  2:2389732/150‑1 (MQ=255)
                                                                                                          aaaCTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGttt                                            <  2:636640/150‑1 (MQ=255)
                                                                                                             ctcaactcTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAAc                                         >  2:3716491/1‑150 (MQ=255)
                                                                                                               caactcTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACtt                                       <  2:737680/150‑1 (MQ=255)
                                                                                                                aactcTCAAATGAACCGCGAGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACttt                                      <  1:282708/150‑1 (MQ=255)
                                                                                                                 actcTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTa                                     <  1:3767534/150‑1 (MQ=255)
                                                                                                                  ctctCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTAt                                    <  2:2868015/150‑1 (MQ=255)
                                                                                                                     tcAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGaa                                 >  2:1839490/1‑150 (MQ=255)
                                                                                                                     tcAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGaa                                 <  1:275937/150‑1 (MQ=255)
                                                                                                                      cAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAag                                >  1:1557630/1‑150 (MQ=255)
                                                                                                                        aaTGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAagag                              <  1:299974/150‑1 (MQ=255)
                                                                                                                         aTGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAgt                             >  1:2157402/1‑150 (MQ=255)
                                                                                                                         aTGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAgt                             >  1:3138397/1‑150 (MQ=255)
                                                                                                                         aTGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAgt                             >  1:2455890/1‑150 (MQ=255)
                                                                                                                          tGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAgtg                            <  2:1008043/150‑1 (MQ=255)
                                                                                                                           gAACCGCGAGTTTTACGCATCCAATCTCGACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAgtgt                           <  2:2478945/150‑1 (MQ=255)
                                                                                                                            aaCCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTc                          <  1:2616212/150‑1 (MQ=255)
                                                                                                                            aaCCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTATCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTc                          >  1:2584677/1‑150 (MQ=255)
                                                                                                                                   aGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAACTGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTa                   <  2:2007519/150‑1 (MQ=255)
                                                                                                                                   aGTTTTACGCATCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTa                   >  2:2090998/1‑150 (MQ=255)
                                                                                                                                           gCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATc           >  2:2776389/1‑150 (MQ=255)
                                                                                                                                             aTCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGt         >  2:2970708/1‑150 (MQ=255)
                                                                                                                                              tCCAATATCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGtt        >  2:1563738/1‑150 (MQ=255)
                                                                                                                                                 aaTCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTaaa     >  2:2880236/1‑150 (MQ=255)
                                                                                                                                                  aTCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAg    <  2:817196/150‑1 (MQ=255)
                                                                                                                                                  aTCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAg    >  1:1047559/1‑150 (MQ=255)
                                                                                                                                                   tctcTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGc   >  1:3694962/1‑150 (MQ=255)
                                                                                                                                                   tctcTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGc   <  2:859808/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCC  >  CP009273/1980324‑1980621

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