Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 634,078 T→A 22.7% intergenic (‑49/‑323) dsbG ← / → ahpC thiol:disulfide interchange protein, periplasmic/alkyl hydroperoxide reductase, C22 subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273634,0780TA22.7% 559.6 / 175.5 264intergenic (‑49/‑323)dsbG/ahpCthiol:disulfide interchange protein, periplasmic/alkyl hydroperoxide reductase, C22 subunit
Reads supporting (aligned to +/- strand):  ref base T (95/109);  new base A (30/30);  total (125/139)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

Only 100 of 266 total aligned reads displayed.

TGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCAGATTA  >  CP009273/633931‑634221
                                                                                                                                                   |                                                                                                                                               
tGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCt                                                                                                                                                 >  2:1313807/1‑148 (MQ=255)
  tAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTAt                                                                                                                                               <  2:2878886/148‑1 (MQ=255)
  tAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAAt                                                                                                                                               >  1:2381280/1‑148 (MQ=255)
    aTGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCa                                                                                                                                             <  2:1121662/148‑1 (MQ=255)
    aTGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCa                                                                                                                                             <  1:870257/148‑1 (MQ=255)
        ccTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGt                                                                                                                                         <  1:164290/148‑1 (MQ=255)
          tGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCa                                                                                                                                       >  1:1166403/1‑148 (MQ=255)
          tGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCa                                                                                                                                       >  1:1260605/1‑148 (MQ=255)
            tttttCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATa                                                                                                                                     >  2:1276610/1‑148 (MQ=255)
            tttttCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATa                                                                                                                                     <  1:3370542/148‑1 (MQ=255)
               ttCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGt                                                                                                                                  <  1:656383/148‑1 (MQ=255)
                tCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTa                                                                                                                                 <  1:1870402/148‑1 (MQ=255)
                   aTCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTCaga                                                                                                                              >  2:2147845/1‑148 (MQ=255)
                   aTCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAaga                                                                                                                              >  1:3297859/1‑148 (MQ=255)
                    tCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAagag                                                                                                                             >  1:2833208/1‑148 (MQ=255)
                     cGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGc                                                                                                                            <  1:437773/148‑1 (MQ=255)
                     cGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGc                                                                                                                            <  1:2359459/148‑1 (MQ=255)
                     cGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGc                                                                                                                            <  1:2077832/148‑1 (MQ=255)
                     cGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGc                                                                                                                            <  2:1418457/148‑1 (MQ=255)
                       cTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCtt                                                                                                                          <  1:3516431/148‑1 (MQ=255)
                        ttttACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTa                                                                                                                         >  2:193259/1‑148 (MQ=255)
                          ttACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGa                                                                                                                       <  1:3512204/148‑1 (MQ=255)
                          ttACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGa                                                                                                                       <  2:1622328/148‑1 (MQ=255)
                             cTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCa                                                                                                                    <  1:106758/148‑1 (MQ=255)
                              tGGATCAGGAAGTTCCTCTGCGAAGGCGATTGAAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAg                                                                                                                   >  2:1845391/1‑148 (MQ=255)
                              tGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCAGTGTCCTAGGTAAGAGCTTAGATCAg                                                                                                                   >  2:1844764/1‑148 (MQ=255)
                               ggAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAgg                                                                                                                  <  2:2738811/148‑1 (MQ=255)
                                gAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGt                                                                                                                 >  1:2591690/1‑148 (MQ=255)
                                   cAGGAAGGTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGAt                                                                                                              <  2:60258/148‑1 (MQ=255)
                                       aaGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGcc                                                                                                          <  2:1646443/148‑1 (MQ=255)
                                       aaGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGcc                                                                                                          <  2:2652000/148‑1 (MQ=255)
                                        aGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGccc                                                                                                         <  2:1026967/148‑1 (MQ=255)
                                        aGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGccc                                                                                                         <  1:3621503/148‑1 (MQ=255)
                                        aGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGccc                                                                                                         <  1:1729007/148‑1 (MQ=255)
                                           tCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCttt                                                                                                      >  2:1231756/1‑148 (MQ=255)
                                             ctctGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGt                                                                                                    <  1:2181587/148‑1 (MQ=255)
                                             ctctGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGt                                                                                                    <  2:2078630/148‑1 (MQ=255)
                                                          aTTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGtgtt                                                                                       <  1:342369/148‑1 (MQ=255)
                                                          aTTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGtgtt                                                                                       <  1:2972907/148‑1 (MQ=255)
                                                           ttGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGtgttg                                                                                      <  2:3297859/148‑1 (MQ=255)
                                                                ggAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATc                                                                                 >  2:2262112/1‑148 (MQ=255)
                                                                ggAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATc                                                                                 >  1:3145032/1‑148 (MQ=255)
                                                                 gAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATcc                                                                                >  1:430061/1‑148 (MQ=255)
                                                                 gAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATcc                                                                                <  1:1150571/148‑1 (MQ=255)
                                                                 gAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATcc                                                                                <  1:1353412/148‑1 (MQ=255)
                                                                 gAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATcc                                                                                >  2:3144650/1‑148 (MQ=255)
                                                                 gAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATcc                                                                                >  2:1721445/1‑148 (MQ=255)
                                                                  aaGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCa                                                                               <  2:2930629/148‑1 (MQ=255)
                                                                          ccAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCgttg                                                                       <  1:2537684/148‑1 (MQ=255)
                                                                           cAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCgttgt                                                                      <  2:3498440/148‑1 (MQ=255)
                                                                              tAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGc                                                                   >  1:156224/1‑148 (MQ=255)
                                                                               aaaaGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCa                                                                  >  1:2746820/1‑148 (MQ=255)
                                                                                  aGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCAttt                                                               >  2:3807129/1‑148 (MQ=255)
                                                                                   gTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTg                                                              <  2:2109528/148‑1 (MQ=255)
                                                                                   gTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTg                                                              >  1:1658473/1‑148 (MQ=255)
                                                                                    tATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGt                                                             <  1:2415275/148‑1 (MQ=255)
                                                                                    tATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGt                                                             >  2:3595117/1‑148 (MQ=255)
                                                                                     aTCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTa                                                            <  1:3600383/148‑1 (MQ=255)
                                                                                     aTCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTa                                                            >  2:2635355/1‑148 (MQ=255)
                                                                                     aTCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTa                                                            <  2:3503636/148‑1 (MQ=255)
                                                                                      tCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTaa                                                           <  2:197710/148‑1 (MQ=255)
                                                                                      tCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTaa                                                           <  2:748881/148‑1 (MQ=255)
                                                                                       cTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAg                                                          <  1:1828685/148‑1 (MQ=255)
                                                                                       cTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAg                                                          >  2:576975/1‑148 (MQ=255)
                                                                                       cTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAg                                                          <  1:2711127/148‑1 (MQ=255)
                                                                                         tttttAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAggg                                                        <  1:3261403/148‑1 (MQ=255)
                                                                                         tttttAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAggg                                                        >  2:1770460/1‑148 (MQ=255)
                                                                                               aCATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAacac                                                  <  2:3018740/148‑1 (MQ=255)
                                                                                               aCATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAacac                                                  >  1:489837/1‑148 (MQ=255)
                                                                                                cATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACAcc                                                 <  1:194023/148‑1 (MQ=255)
                                                                                                 aTTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCt                                                >  1:2507136/1‑148 (MQ=255)
                                                                                                     aTTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGc                                            >  1:388218/1‑148 (MQ=255)
                                                                                                     aTTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGc                                            >  1:316542/1‑148 (MQ=255)
                                                                                                      tttGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGcc                                           >  1:3181960/1‑148 (MQ=255)
                                                                                                      tttGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGcc                                           <  1:962472/148‑1 (MQ=255)
                                                                                                          tCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTgca                                       >  2:1609652/1‑148 (MQ=255)
                                                                                                          tCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTgca                                       <  1:2286766/148‑1 (MQ=255)
                                                                                                           ccTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTgcag                                      >  2:3682597/1‑148 (MQ=255)
                                                                                                           ccTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTgcag                                      >  2:1791419/1‑148 (MQ=255)
                                                                                                           ccTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTgcag                                      >  2:2061875/1‑148 (MQ=255)
                                                                                                             ttttCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTgcaggc                                    <  2:1708177/148‑1 (MQ=255)
                                                                                                             ttttCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTgcaggc                                    <  2:1708230/148‑1 (MQ=255)
                                                                                                                 cagtcagtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCaggcaggc                                >  1:1471476/1‑148 (MQ=255)
                                                                                                                  agtcagtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAggcaggca                               <  1:1376106/148‑1 (MQ=255)
                                                                                                                  agtcagtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAggcaggca                               <  1:2016387/148‑1 (MQ=255)
                                                                                                                    tcagtGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACt                             >  2:2447796/1‑148 (MQ=255)
                                                                                                                    tcagtGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACt                             >  2:526748/1‑148 (MQ=255)
                                                                                                                      agtGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGa                           <  2:6718/148‑1 (MQ=255)
                                                                                                                       gtGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGaa                          <  2:2381213/148‑1 (MQ=255)
                                                                                                                       gtGCAAAAGTCGAGTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGaa                          <  2:456145/148‑1 (MQ=255)
                                                                                                                        tGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAg                         >  2:2591583/1‑148 (MQ=255)
                                                                                                                               gTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCaa                  >  1:1501966/1‑148 (MQ=255)
                                                                                                                                    gTAAAAGGCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGt             <  2:115778/148‑1 (MQ=255)
                                                                                                                                     tAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTa            >  1:1013949/1‑148 (MQ=255)
                                                                                                                                      aaaaGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAg           <  2:3438194/148‑1 (MQ=255)
                                                                                                                                      aaaaGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAg           >  1:3487962/1‑148 (MQ=255)
                                                                                                                                        aaGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGtt         >  1:290007/1‑148 (MQ=255)
                                                                                                                                           gCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCAg      >  2:1779377/1‑148 (MQ=255)
                                                                                                                                           gCATAACCAATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCAg      <  2:1486320/148‑1 (MQ=255)
                                                                                                                                             aTAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCAGAt    >  2:931898/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGGAAGTTCCTCTGCGAAGGCGATTGCAGGAAGCAGAGCCAGTAAAAGTATCTTTTTTAACATTAATTTGTCCTTTTCAGTCAGTGCAAAAGTCGAGTAAAAGGCATAACCTATCACTGTCATAGGTAAGAGCTTAGATCAGGTGATTGCCCTTTGTTTATGAGGGTGTTGTAATCCATGTCGTTGTTGCATTTGTAAGGGCAACACCTCAGCCTGCAGGCAGGCACTGAAGATACCAAAGGGTAGTTCAGATTA  >  CP009273/633931‑634221

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: