Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,459,603 A→T 40.4% intergenic (+24/+5) gspO → / ← bfr bifunctional prepilin leader peptidase/ methylase/bacterioferritin, iron storage and detoxification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,459,6030AT40.4% 132.5 / 233.5 181intergenic (+24/+5)gspO/bfrbifunctional prepilin leader peptidase/ methylase/bacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (57/51);  new base T (41/32);  total (98/83)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 189 total aligned reads displayed.

CATTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATC  >  CP009273/3459456‑3459748
                                                                                                                                                   |                                                                                                                                                 
cATTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGaca                                                                                                                                                   >  2:465543/1‑148 (MQ=255)
cATTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGaca                                                                                                                                                   >  1:151880/1‑148 (MQ=255)
 aTTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGacac                                                                                                                                                  >  1:1474769/1‑148 (MQ=255)
    ggTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGc                                                                                                                                               >  1:1817931/1‑148 (MQ=255)
     gTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCt                                                                                                                                              >  1:2082058/1‑148 (MQ=255)
     gTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCt                                                                                                                                              >  1:2971216/1‑148 (MQ=255)
     gTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCt                                                                                                                                              >  1:2971290/1‑148 (MQ=255)
       tACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCa                                                                                                                                            <  1:857517/148‑1 (MQ=255)
       tACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCa                                                                                                                                            >  2:633289/1‑148 (MQ=255)
       tACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCa                                                                                                                                            <  2:2709058/148‑1 (MQ=255)
            cATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACctt                                                                                                                                       <  2:2442348/148‑1 (MQ=255)
             aTCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACcttc                                                                                                                                      <  2:2973380/148‑1 (MQ=255)
              tCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACcttct                                                                                                                                     <  1:210401/148‑1 (MQ=255)
                  ttCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTcgc                                                                                                                                 >  2:2810698/1‑148 (MQ=255)
                      caACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGAt                                                                                                                             >  2:329383/1‑148 (MQ=255)
                       aaCTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATc                                                                                                                            >  2:2555221/1‑148 (MQ=255)
                        aCTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCt                                                                                                                           >  1:300766/1‑148 (MQ=255)
                          ttCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCtgt                                                                                                                         >  1:1302232/1‑148 (MQ=255)
                             tGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCt                                                                                                                      <  1:239236/148‑1 (MQ=255)
                               cAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTg                                                                                                                    >  2:344680/1‑148 (MQ=255)
                                aaCGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGc                                                                                                                   <  1:2754508/148‑1 (MQ=255)
                                 aCGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCa                                                                                                                  >  1:1905497/1‑148 (MQ=255)
                                 aCGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCa                                                                                                                  >  1:133024/1‑148 (MQ=255)
                                 aCGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCa                                                                                                                  >  2:2066071/1‑148 (MQ=255)
                                  cGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAg                                                                                                                 >  2:1115735/1‑148 (MQ=255)
                                   gCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGa                                                                                                                <  2:2882478/148‑1 (MQ=255)
                                      ggACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATaa                                                                                                             >  1:1785041/1‑148 (MQ=255)
                                       gACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAt                                                                                                            <  2:2684251/148‑1 (MQ=255)
                                        aCTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAtt                                                                                                           <  1:1870234/148‑1 (MQ=255)
                                         cTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAttt                                                                                                          <  1:786834/148‑1 (MQ=255)
                                         cTATAACCACTATCGCAGTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCGGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCGCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAttt                                                                                                          <  2:2397508/148‑1 (MQ=255)
                                          tataACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCCTCTTCGCGGATCTGTGCTTGCAGATAAtttt                                                                                                         >  2:2510880/1‑148 (MQ=255)
                                               ccACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAgg                                                                                                    <  1:286935/148‑1 (MQ=255)
                                                cACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGc                                                                                                   >  2:1376119/1‑148 (MQ=255)
                                                cACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGc                                                                                                   <  1:228075/148‑1 (MQ=255)
                                                   tATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCa                                                                                                >  1:3258205/1‑148 (MQ=255)
                                                           tCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGa                                                                                        >  1:2504873/1‑148 (MQ=255)
                                                           tCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGa                                                                                        >  1:1473152/1‑148 (MQ=255)
                                                           tCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGa                                                                                        >  1:1255091/1‑148 (MQ=255)
                                                              gCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATc                                                                                     <  2:1920404/148‑1 (MQ=255)
                                                               cccGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCa                                                                                    <  2:1503684/148‑1 (MQ=255)
                                                                ccGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAg                                                                                   <  1:1088019/148‑1 (MQ=255)
                                                                   tGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATc                                                                                >  2:2521390/1‑148 (MQ=255)
                                                                    ggCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCa                                                                               >  2:435558/1‑148 (MQ=255)
                                                                           gcTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCg                                                                        <  1:524381/148‑1 (MQ=255)
                                                                               ggCAGCATGATAAACCTGGGGTATCTGGCGTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGTTCATGTTCCGTTTc                                                                    <  2:3018970/148‑1 (MQ=255)
                                                                               ggCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGATCCGTTTc                                                                    >  1:447748/1‑148 (MQ=255)
                                                                                gcagcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTcc                                                                   <  1:465543/148‑1 (MQ=255)
                                                                                gcagcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTcc                                                                   <  2:403660/148‑1 (MQ=255)
                                                                                 cagcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGATGGTTGACCCTCCCGACTCCGCTCAACCTTCGTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTcca                                                                  <  2:1803549/148‑1 (MQ=255)
                                                                                   gcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagc                                                                >  2:1449002/1‑148 (MQ=255)
                                                                                   gcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagc                                                                <  2:2382219/148‑1 (MQ=255)
                                                                                    caTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagcc                                                               >  1:3193507/1‑148 (MQ=255)
                                                                                    caTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagcc                                                               >  1:2066731/1‑148 (MQ=255)
                                                                                    caTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAAATTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTccagcc                                                               >  1:2021825/1‑148 (MQ=255)
                                                                                         tAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGAAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCg                                                          >  1:3046608/1‑148 (MQ=255)
                                                                                          aaaCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGa                                                         >  2:2308348/1‑148 (MQ=255)
                                                                                           aaCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGat                                                        >  2:926045/1‑148 (MQ=255)
                                                                                           aaCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGat                                                        <  2:3290590/148‑1 (MQ=255)
                                                                                            aCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGata                                                       >  2:860587/1‑148 (MQ=255)
                                                                                            aCCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGata                                                       >  2:3216081/1‑148 (MQ=255)
                                                                                               tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATgg                                                    >  1:405608/1‑148 (MQ=255)
                                                                                               tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATgg                                                    <  2:1243139/148‑1 (MQ=255)
                                                                                               tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATgg                                                    <  2:1243146/148‑1 (MQ=255)
                                                                                                   gTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCTATATGGCctt                                                <  2:1474106/148‑1 (MQ=255)
                                                                                                   gTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCctt                                                <  1:1548637/148‑1 (MQ=255)
                                                                                                   gTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCctt                                                <  1:147351/148‑1 (MQ=255)
                                                                                                             tGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                      <  1:1140645/148‑1 (MQ=255)
                                                                                                              ggATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGc                                     >  2:2560707/1‑148 (MQ=255)
                                                                                                                aTCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaa                                   >  2:2256686/1‑148 (MQ=255)
                                                                                                                aTCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATATTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaa                                   >  2:1196560/1‑148 (MQ=255)
                                                                                                                 tctcTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaa                                  >  1:2718805/1‑148 (MQ=255)
                                                                                                                  ctctTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaaa                                 >  2:497483/1‑148 (MQ=255)
                                                                                                                  ctctTATTAACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaaa                                 <  2:1762877/148‑1 (MQ=255)
                                                                                                                     ttattaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAAttt                              <  1:1425310/148‑1 (MQ=255)
                                                                                                                      tattaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                             <  1:2266883/148‑1 (MQ=255)
                                                                                                                        ttaACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                           >  1:2769689/1‑148 (MQ=255)
                                                                                                                        ttaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTa                           >  2:82033/1‑148 (MQ=255)
                                                                                                                         taACCGGGAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTat                          >  1:104403/1‑148 (MQ=255)
                                                                                                                         taACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTat                          <  1:479577/148‑1 (MQ=255)
                                                                                                                               ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATcatcat                    <  2:1817931/148‑1 (MQ=255)
                                                                                                                                gAGGGTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCata                   <  1:2689042/148‑1 (MQ=255)
                                                                                                                                gAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCata                   >  1:843956/1‑148 (MQ=255)
                                                                                                                                    gTTCTCCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcgc               >  2:2943189/1‑148 (MQ=255)
                                                                                                                                      tctcCCTCCCGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCgg             <  1:1592958/148‑1 (MQ=255)
                                                                                                                                      tctcCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCgg             <  1:2543295/148‑1 (MQ=255)
                                                                                                                                       ctccctccCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc            <  1:11274/148‑1 (MQ=255)
                                                                                                                                        tccctcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCt           >  1:2772454/1‑148 (MQ=255)
                                                                                                                                          cctcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGa         <  2:300766/148‑1 (MQ=255)
                                                                                                                                          cctcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGa         >  2:350079/1‑148 (MQ=255)
                                                                                                                                          cctcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGa         >  1:2939348/1‑148 (MQ=255)
                                                                                                                                          cctcccGACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGa         <  2:1583371/148‑1 (MQ=255)
                                                                                                                                          cctcccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGa         >  1:1214822/1‑148 (MQ=255)
                                                                                                                                            tcccGACACGGCTCAACCTTCTTCGCGGATCTGTGCGTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACg       <  2:2199329/148‑1 (MQ=255)
                                                                                                                                              ccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTa     <  1:329383/148‑1 (MQ=255)
                                                                                                                                               cGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTaa    <  2:447748/148‑1 (MQ=255)
                                                                                                                                                gACTCGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAAt   <  2:1390959/148‑1 (MQ=255)
                                                                                                                                                 acacGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATc  <  2:1057077/148‑1 (MQ=255)
                                                                                                                                                 acacGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATc  <  2:1302232/148‑1 (MQ=255)
                                                                                                                                                 acacGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGNTCCGTTTGCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATc  >  1:1683119/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CATTGGTTACGCCATCGTTTCACAACTTCTGCAACGCCGGACTATAACCACTATCGCATTCGGCCCGTGGCTTGCGCTCGGCAGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATC  >  CP009273/3459456‑3459748

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: