Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,478,971 G→A 100% E552K (GAG→AAG)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,478,9710GA96.0% 640.5 / 13.4 200E552K (GAG→AAG) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base A (97/94);  total (102/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 200 total aligned reads displayed.

GTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACC  >  CP009273/1478828‑1479114
                                                                                                                                               |                                                                                                                                               
  aTGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTg                                                                                                                                           >  2:2303045/1‑148 (MQ=255)
   tGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGa                                                                                                                                          <  1:1654403/148‑1 (MQ=255)
    ggTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGAc                                                                                                                                         <  1:1506443/148‑1 (MQ=255)
     gTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACt                                                                                                                                        <  1:2175363/148‑1 (MQ=255)
       gCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACttt                                                                                                                                      <  1:1209803/148‑1 (MQ=255)
       gCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACttt                                                                                                                                      >  2:2671169/1‑148 (MQ=255)
       gCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACttt                                                                                                                                      <  1:3182735/148‑1 (MQ=255)
       gCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACttt                                                                                                                                      <  1:3182234/148‑1 (MQ=255)
        cTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtc                                                                                                                                     <  2:674607/148‑1 (MQ=255)
         tGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGAACAAAAGTCAGACtttat                                                                                                                                    >  1:1356572/1‑146 (MQ=255)
           gAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcg                                                                                                                                  >  1:1181665/1‑148 (MQ=255)
           gAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcg                                                                                                                                  >  1:1529853/1‑148 (MQ=255)
           gAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcg                                                                                                                                  >  1:3086889/1‑148 (MQ=255)
            aaGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgc                                                                                                                                 <  1:374355/148‑1 (MQ=255)
            aaGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgc                                                                                                                                 <  1:2588056/148‑1 (MQ=255)
            aaGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgc                                                                                                                                 >  1:1766507/1‑148 (MQ=255)
             aGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgcg                                                                                                                                >  2:1034049/1‑148 (MQ=255)
              gcgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGt                                                                                                                               <  2:2888711/148‑1 (MQ=255)
              gcgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGt                                                                                                                               <  2:1288921/148‑1 (MQ=255)
                  aaaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTtgt                                                                                                                           >  1:2621228/1‑148 (MQ=255)
                  aaaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTtgt                                                                                                                           >  2:3248751/1‑148 (MQ=255)
                  aaaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTtgt                                                                                                                           >  2:1346086/1‑148 (MQ=255)
                   aaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTtgtg                                                                                                                          >  1:2946900/1‑148 (MQ=255)
                    aaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGa                                                                                                                         <  2:3592245/148‑1 (MQ=255)
                    aaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGa                                                                                                                         >  2:1437534/1‑148 (MQ=255)
                     aaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGaa                                                                                                                        <  2:326280/148‑1 (MQ=255)
                      aCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAAt                                                                                                                       <  2:3492728/148‑1 (MQ=255)
                         tGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGTTGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAACAGTCTGACTTTCTCGCGTGTGGGACTCtg                                                                                                                    >  2:1006810/1‑148 (MQ=255)
                         tGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCtg                                                                                                                    <  2:3434216/148‑1 (MQ=255)
                          ggCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtgt                                                                                                                   <  2:166166/148‑1 (MQ=255)
                                gtgcgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAAtt                                                                                                             <  2:3461603/148‑1 (MQ=255)
                                   cgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATc                                                                                                          <  1:2809431/148‑1 (MQ=255)
                                    gtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCt                                                                                                         <  2:429211/148‑1 (MQ=255)
                                    gtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCt                                                                                                         <  1:3468811/148‑1 (MQ=255)
                                    gtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCAGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCt                                                                                                         <  2:2351321/148‑1 (MQ=255)
                                     tgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCtt                                                                                                        <  1:715658/148‑1 (MQ=255)
                                         gCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCg                                                                                                    >  2:1564069/1‑148 (MQ=255)
                                             tgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagca                                                                                                >  1:3811554/1‑148 (MQ=255)
                                              gatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagcag                                                                                               >  1:1405544/1‑148 (MQ=255)
                                              gatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTCCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagcag                                                                                               >  1:3308898/1‑148 (MQ=255)
                                                ttATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACACAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAgcagca                                                                                             >  2:2040104/1‑148 (MQ=255)
                                                ttATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAgcagca                                                                                             >  2:2039937/1‑148 (MQ=255)
                                                ttATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAgcagca                                                                                             >  2:3674211/1‑148 (MQ=255)
                                                ttATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAgcagca                                                                                             >  2:3673709/1‑148 (MQ=255)
                                                 tATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCaa                                                                                            >  2:732230/1‑148 (MQ=255)
                                                        tCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGgcgc                                                                                     <  1:2684075/148‑1 (MQ=255)
                                                        tCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGgcgc                                                                                     >  1:1963353/1‑148 (MQ=255)
                                                          cgcgCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCtt                                                                                   <  1:3346724/148‑1 (MQ=255)
                                                            cgcTCTCCATTCAGGATCCGCGCGAACGTCAGATGGACAAACAGCAGGCATAGGACTAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttc                                                                                 <  1:3722071/148‑1 (MQ=255)
                                                              ctctCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttctt                                                                               >  1:1041404/1‑148 (MQ=255)
                                                              ctctCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttctt                                                                               >  1:2903804/1‑148 (MQ=255)
                                                               tctcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttcttc                                                                              <  1:774547/148‑1 (MQ=255)
                                                               tctcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttcttc                                                                              <  2:1775905/148‑1 (MQ=255)
                                                                ctcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTcc                                                                             >  2:1719116/1‑148 (MQ=255)
                                                                 tcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCa                                                                            >  1:3440981/1‑148 (MQ=255)
                                                                 tcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCa                                                                            >  2:2113524/1‑148 (MQ=255)
                                                                  ccATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCaa                                                                           >  2:1338890/1‑148 (MQ=255)
                                                                      tCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGcc                                                                       <  2:3793066/148‑1 (MQ=255)
                                                                         ggATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTc                                                                    <  2:79036/148‑1 (MQ=255)
                                                                         ggATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTc                                                                    <  2:2859899/148‑1 (MQ=255)
                                                                          gATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTcc                                                                   <  1:3374734/148‑1 (MQ=255)
                                                                          gATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTcc                                                                   >  2:3333992/1‑148 (MQ=255)
                                                                          gATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTcc                                                                   >  2:1885982/1‑148 (MQ=255)
                                                                             ccGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTc                                                                >  1:2923609/1‑148 (MQ=255)
                                                                                 gcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCGTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCt                                                            <  2:1493932/148‑1 (MQ=255)
                                                                                  cgAACGTCCGATGGACAAACAGTAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCtg                                                           <  2:84856/148‑1 (MQ=255)
                                                                                  cgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCAtg                                                           <  1:641179/148‑1 (MQ=255)
                                                                                   gAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCtgt                                                          <  2:2452944/148‑1 (MQ=255)
                                                                                    aaCGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCtgtg                                                         <  2:3007534/148‑1 (MQ=255)
                                                                                     aCGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTgtgt                                                        >  1:2045198/1‑148 (MQ=255)
                                                                                      cGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTc                                                       >  2:33823/1‑148 (MQ=255)
                                                                                      cGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTc                                                       >  1:2721302/1‑148 (MQ=255)
                                                                                        tCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGt                                                     >  1:1266533/1‑148 (MQ=255)
                                                                                         ccGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTa                                                    >  1:2346169/1‑148 (MQ=255)
                                                                                            aTGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCg                                                 <  2:1464374/148‑1 (MQ=255)
                                                                                              ggACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGAt                                               >  2:285164/1‑148 (MQ=255)
                                                                                               gACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGAtt                                              <  2:932362/148‑1 (MQ=255)
                                                                                                   aaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtct                                          >  2:966755/1‑148 (MQ=255)
                                                                                                    aCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtctc                                         >  1:356518/1‑148 (MQ=255)
                                                                                                    aCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtctc                                         >  1:947814/1‑148 (MQ=255)
                                                                                                    aCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAAGTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtctc                                         >  2:16151/1‑148 (MQ=255)
                                                                                                        cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                     <  1:222017/148‑1 (MQ=255)
                                                                                                        cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                     <  2:1208215/148‑1 (MQ=255)
                                                                                                        cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                     <  2:1165361/148‑1 (MQ=255)
                                                                                                          ggCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTAt                                   >  1:475314/1‑148 (MQ=255)
                                                                                                                      aaaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcgc                       >  2:1133063/1‑148 (MQ=255)
                                                                                                                      aaaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcgc                       >  2:2952787/1‑148 (MQ=255)
                                                                                                                       aaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcgcg                      >  2:1927256/1‑148 (MQ=255)
                                                                                                                        aaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGa                     <  1:2303745/148‑1 (MQ=255)
                                                                                                                        aaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGa                     <  1:618596/148‑1 (MQ=255)
                                                                                                                         aCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGaa                    >  2:624636/1‑148 (MQ=255)
                                                                                                                          cATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGGGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAAt                   >  1:3000772/1‑148 (MQ=255)
                                                                                                                             cgtcgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCATACAAAGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACAATCTGCGCGTGCGCGAATGGc                >  2:3545705/1‑148 (MQ=255)
                                                                                                                                cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAgg             <  1:127563/148‑1 (MQ=255)
                                                                                                                                cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAgg             >  1:232080/1‑148 (MQ=255)
                                                                                                                                 gCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGa            <  1:3673709/148‑1 (MQ=255)
                                                                                                                                 gCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGa            <  1:3674211/148‑1 (MQ=255)
                                                                                                                                  cTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGat           >  2:1400517/1‑148 (MQ=255)
                                                                                                                                  cTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGat           >  2:2141758/1‑148 (MQ=255)
                                                                                                                                   ttCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGata          <  2:1654499/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACC  >  CP009273/1478828‑1479114

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: ‑