Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,569,865 G→T 5.4% intergenic (+137/+20) symR → / ← hsdS sRNA antisense regulator destabilizing divergent and overlapping symE mRNA/specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,569,8650GT 5.4% 1053.3 / 23.6 270intergenic (+137/+20)symR/hsdSsRNA antisense regulator destabilizing divergent and overlapping symE mRNA/specificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base G (159/96);  new base T (7/8);  total (166/104)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

TCCGCACCGCGCTAAATACCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTA  >  CP009273/4569716‑4570014
                                                                                                                                                     |                                                                                                                                                     
tCCGCACCGCGCTAAATACCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAg                                                                                                                                                       <  2:4046046/150‑1 (MQ=255)
     aCCGCGCTAAATACCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATCTTTATGTGGGTTATGAGTGGCGTGAAATAATAAAAAGCGCACCGGACAGGTgcg                                                                                                                                                  >  2:3459701/1‑150 (MQ=255)
     aCCGCGCTAAATACCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTgcg                                                                                                                                                  <  2:1210463/150‑1 (MQ=255)
     aCCGCGCTAAATACCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTgcg                                                                                                                                                  >  1:131075/1‑150 (MQ=255)
               aTACCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAATGTGCGCCAGAAaata                                                                                                                                        <  2:2809831/150‑1 (MQ=255)
                    tGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAATTATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGtt                                                                                                                                   >  2:3867443/1‑150 (MQ=255)
                        tatatCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAgg                                                                                                                               <  1:10098/150‑1 (MQ=255)
                          tatCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGAt                                                                                                                             >  2:1881026/1‑150 (MQ=255)
                           atcatcTGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATTTTTTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGAtt                                                                                                                            <  1:155079/150‑1 (MQ=255)
                           atcatcAGTAAATAGGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTACATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGAtt                                                                                                                            <  2:1297837/150‑1 (MQ=255)
                             catcaGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAATGTGCGCCAGAAAATAATGTTCAGGAtttt                                                                                                                          >  2:2876553/1‑150 (MQ=255)
                              atcaGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGAttttt                                                                                                                         >  1:2781727/1‑150 (MQ=255)
                                caGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTa                                                                                                                       >  1:1244048/1‑150 (MQ=255)
                                 aGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTAc                                                                                                                      >  2:2928751/1‑150 (MQ=255)
                                 aGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTAc                                                                                                                      <  1:2319782/150‑1 (MQ=255)
                                  gTAAATATGGGGAAAGTCCAGCTAAAAATAGACTAAAATGGTCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACg                                                                                                                     <  1:3306779/150‑1 (MQ=255)
                                  gTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAACGGTGCGCCAGACACTAATGTTCAGGATTTTTTACg                                                                                                                     >  2:3055899/1‑150 (MQ=255)
                                    aaaTATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGACTGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTg                                                                                                                   >  1:4016296/1‑150 (MQ=255)
                                    aaaTATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTg                                                                                                                   >  1:3932247/1‑150 (MQ=255)
                                     aaTATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGa                                                                                                                  >  2:366677/1‑150 (MQ=255)
                                     aaTATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGa                                                                                                                  >  2:863083/1‑150 (MQ=255)
                                        atGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGc                                                                                                               >  2:2572826/1‑150 (MQ=255)
                                            ggAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCtttt                                                                                                           >  2:4128105/1‑150 (MQ=255)
                                             gAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCttttt                                                                                                          >  1:1729154/1‑150 (MQ=255)
                                               aaGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAATGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTa                                                                                                        <  2:2823541/150‑1 (MQ=255)
                                                 gTCCAGCTAAAAATAGAATAACATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAcc                                                                                                      <  2:559699/150‑1 (MQ=255)
                                                 gTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAcc                                                                                                      <  1:125201/150‑1 (MQ=255)
                                                  tCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAccc                                                                                                     <  2:1560667/150‑1 (MQ=255)
                                                  tCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAccc                                                                                                     <  2:3295436/150‑1 (MQ=255)
                                                    cAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAccccc                                                                                                   >  2:2534696/1‑150 (MQ=255)
                                                     aGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTATACGTGAGGCTTTTTTACCCCCg                                                                                                  >  1:2094201/1‑150 (MQ=255)
                                                        tAAAAATAGACTAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTa                                                                                               <  1:2564943/150‑1 (MQ=255)
                                                         aaaaaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAg                                                                                              <  1:492289/150‑1 (MQ=255)
                                                          aaaaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGc                                                                                             >  1:1546899/1‑150 (MQ=255)
                                                          aaaaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGc                                                                                             >  1:3771195/1‑150 (MQ=255)
                                                            aaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAATGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTg                                                                                           <  2:1245200/150‑1 (MQ=255)
                                                             aTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAATGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACACCCGCTAGCTgc                                                                                          >  1:3437984/1‑150 (MQ=255)
                                                               aGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTgcgc                                                                                        >  1:3136066/1‑150 (MQ=255)
                                                                gAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGcgcg                                                                                       >  2:3965976/1‑150 (MQ=255)
                                                                    aacaTGGGCAAGTTCTGGAATGATTTAAAGATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCa                                                                                   <  2:2664091/147‑1 (MQ=255)
                                                                    aaaaTGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCa                                                                                   >  2:1531254/1‑150 (MQ=255)
                                                                        tGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCtt                                                                               >  1:2539569/1‑150 (MQ=255)
                                                                        tGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCtt                                                                               >  1:870812/1‑150 (MQ=255)
                                                                          ggCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTg                                                                             <  1:1800011/150‑1 (MQ=255)
                                                                            cAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGAt                                                                           >  2:2229381/1‑150 (MQ=255)
                                                                             aaTTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGAtt                                                                          <  2:1730842/150‑1 (MQ=255)
                                                                             aaTTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGAtt                                                                          >  1:2920306/1‑150 (MQ=255)
                                                                               tttCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGAtttt                                                                        <  1:1980271/150‑1 (MQ=255)
                                                                                ttCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGAttttt                                                                       <  1:2456378/150‑1 (MQ=255)
                                                                                   tGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTcc                                                                    <  2:1359025/150‑1 (MQ=255)
                                                                                   tGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTcc                                                                    >  2:1010042/1‑150 (MQ=255)
                                                                                    ggAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCa                                                                   >  1:899208/1‑150 (MQ=255)
                                                                                    ggAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCa                                                                   >  1:2224791/1‑150 (MQ=255)
                                                                                     gAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAg                                                                  >  1:2925028/1‑150 (MQ=255)
                                                                                      aaTGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGc                                                                 >  1:2081821/1‑150 (MQ=255)
                                                                                       aTGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCa                                                                <  1:1374953/150‑1 (MQ=255)
                                                                                       aTGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCa                                                                >  2:1670263/1‑150 (MQ=255)
                                                                                       aTGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCa                                                                >  2:1671277/1‑150 (MQ=255)
                                                                                       aTGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGATTTGATTTTTTCCAGCa                                                                >  2:1798639/1‑150 (MQ=255)
                                                                                        tGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCaa                                                               <  2:3538796/150‑1 (MQ=255)
                                                                                           tttAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAAcgc                                                            >  2:3584426/1‑150 (MQ=255)
                                                                                             tAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcgg                                                          >  1:132335/1‑150 (MQ=255)
                                                                                              aaaTATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggc                                                         <  1:3630056/150‑1 (MQ=255)
                                                                                              aaaTATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggc                                                         >  1:567359/1‑150 (MQ=255)
                                                                                                atatatTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTATTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggcgg                                                       >  2:1577264/1‑150 (MQ=255)
                                                                                                 tagatTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGCAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggcggc                                                      <  2:2308348/147‑1 (MQ=255)
                                                                                                  atatTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAAGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggcggcg                                                     >  1:297676/1‑150 (MQ=255)
                                                                                                   tatTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGgcgc                                                    >  1:4131017/1‑150 (MQ=255)
                                                                                                    atTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCt                                                   >  2:3926140/1‑150 (MQ=255)
                                                                                                            gggTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTtctc                                           >  1:916005/1‑150 (MQ=255)
                                                                                                              gTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATTATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg                                         <  1:3503617/150‑1 (MQ=255)
                                                                                                              gTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg                                         >  2:4226724/1‑150 (MQ=255)
                                                                                                               ttATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGc                                        >  2:2167770/1‑150 (MQ=255)
                                                                                                                 aTGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTg                                      <  1:4183763/150‑1 (MQ=255)
                                                                                                                 aTGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTg                                      >  1:515197/1‑150 (MQ=255)
                                                                                                                 aTGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTg                                      <  1:4044678/150‑1 (MQ=255)
                                                                                                                   gATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGAt                                    >  1:4352839/1‑150 (MQ=255)
                                                                                                                    aTTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATc                                   <  2:2085744/150‑1 (MQ=255)
                                                                                                                          gTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGCTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATcc                             >  2:4304971/1‑150 (MQ=255)
                                                                                                                           tGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCg                            >  1:389840/1‑150 (MQ=255)
                                                                                                                             aaaTAATAAAAAGCGCACCGGAAATGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCggg                          >  1:2339572/1‑150 (MQ=255)
                                                                                                                             aaaTAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCggg                          >  2:1438267/1‑150 (MQ=255)
                                                                                                                               ataataAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGtt                        <  1:1139966/150‑1 (MQ=255)
                                                                                                                                taataaAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGttt                       <  1:3825699/150‑1 (MQ=255)
                                                                                                                                taataaAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGttt                       >  1:3585425/1‑150 (MQ=255)
                                                                                                                                taataaAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGttt                       >  2:519513/1‑150 (MQ=255)
                                                                                                                                  ataaAAAGCGCACCGGAAATGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTc                     >  2:3470825/1‑150 (MQ=255)
                                                                                                                                  ataaAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTc                     <  2:44263/150‑1 (MQ=255)
                                                                                                                                  ataaAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTc                     >  2:3242732/1‑150 (MQ=255)
                                                                                                                                   taaAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCg                    <  2:1644096/150‑1 (MQ=255)
                                                                                                                                      aaaGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGcc                 >  2:1263761/1‑150 (MQ=255)
                                                                                                                                       aaGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGccc                <  1:3955011/150‑1 (MQ=255)
                                                                                                                                       aaGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGccc                >  2:2390767/1‑150 (MQ=255)
                                                                                                                                          cgcACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGcc             >  1:646158/1‑150 (MQ=255)
                                                                                                                                              ccGGAAAGGTGCGACAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTg         <  1:1660104/150‑1 (MQ=255)
                                                                                                                                               cGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg        >  2:4226492/1‑150 (MQ=255)
                                                                                                                                                ggAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTggg       <  1:1495948/150‑1 (MQ=255)
                                                                                                                                                  aaaGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCg     <  2:4198313/150‑1 (MQ=255)
                                                                                                                                                   aaGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCgg    >  1:2838234/1‑150 (MQ=255)
                                                                                                                                                     ggTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTa  >  2:1849815/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCGCACCGCGCTAAATACCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTA  >  CP009273/4569716‑4570014

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