Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,088,065 C→A 90.6% V125V (GTG→GTT yggR ← putative PilT family AAA+ ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,088,0650CA90.6% 807.5 / 77.8 285V125V (GTG→GTTyggRputative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base C (16/11);  new base A (139/119);  total (155/130)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

TCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCG  >  CP009273/3087919‑3088211
                                                                                                                                                  |                                                                                                                                                  
  ccGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAg                                                                                                                                               <  1:952602/150‑1 (MQ=255)
    gctgctGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGaa                                                                                                                                             >  2:4497725/1‑150 (MQ=255)
     ctgctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAAt                                                                                                                                            >  1:4138683/1‑150 (MQ=255)
      tgctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATc                                                                                                                                           >  1:3914711/1‑150 (MQ=255)
       gctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGTGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCa                                                                                                                                          <  1:1143599/150‑1 (MQ=255)
       gctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCa                                                                                                                                          <  1:4037345/150‑1 (MQ=255)
        ctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAg                                                                                                                                         <  2:4431505/150‑1 (MQ=255)
           gATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGcc                                                                                                                                      >  2:3334008/1‑150 (MQ=255)
               aaaCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATtc                                                                                                                                  <  1:463325/150‑1 (MQ=255)
                aaCATCGCTTGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATtct                                                                                                                                 >  1:747223/1‑150 (MQ=255)
                aaCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATtct                                                                                                                                 >  1:1439612/1‑150 (MQ=255)
                 aCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGACAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATtctc                                                                                                                                <  1:4616873/150‑1 (MQ=255)
                  cATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCg                                                                                                                               <  1:3021892/150‑1 (MQ=255)
                  cATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCg                                                                                                                               >  1:4137309/1‑150 (MQ=255)
                     cGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGctc                                                                                                                            >  2:2270852/1‑150 (MQ=255)
                     cGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGctc                                                                                                                            <  1:673568/150‑1 (MQ=255)
                         ggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGa                                                                                                                        <  1:4339902/150‑1 (MQ=255)
                         ggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGa                                                                                                                        <  1:3575565/150‑1 (MQ=255)
                             ggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTaa                                                                                                                    >  2:768747/1‑150 (MQ=255)
                                  agagATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCg                                                                                                               <  2:3844025/150‑1 (MQ=255)
                                   gagaTATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCgg                                                                                                              >  1:3959468/1‑150 (MQ=255)
                                      atatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCaa                                                                                                           >  1:2660005/1‑150 (MQ=255)
                                      atatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCaa                                                                                                           <  2:1992842/150‑1 (MQ=255)
                                      atatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCaa                                                                                                           >  2:1698730/1‑150 (MQ=255)
                                       tatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAAt                                                                                                          <  2:1451487/150‑1 (MQ=255)
                                        atTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATa                                                                                                         >  1:3800700/1‑150 (MQ=255)
                                            aaCAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATAAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGt                                                                                                     >  2:3722833/1‑150 (MQ=255)
                                               aGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTgg                                                                                                  >  1:4144282/1‑150 (MQ=255)
                                               aGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTAGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAAACAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACAGTTgg                                                                                                  >  1:2214751/1‑150 (MQ=255)
                                                 gATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTg                                                                                                >  1:2367053/1‑150 (MQ=255)
                                                  aTCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgc                                                                                               <  1:2832483/150‑1 (MQ=255)
                                                   tCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcg                                                                                              <  1:3951440/150‑1 (MQ=255)
                                                    cTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcgc                                                                                             >  2:481853/1‑150 (MQ=255)
                                                       ccAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCaa                                                                                          <  1:2453521/150‑1 (MQ=255)
                                                       ccAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCaa                                                                                          <  1:717179/150‑1 (MQ=255)
                                                            gTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGt                                                                                     >  1:568790/1‑150 (MQ=255)
                                                             tCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGtt                                                                                    <  2:182617/150‑1 (MQ=255)
                                                              cAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTc                                                                                   >  1:2769962/1‑150 (MQ=255)
                                                               aGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCg                                                                                  >  2:3491180/1‑150 (MQ=255)
                                                               aGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCg                                                                                  >  2:4479459/1‑150 (MQ=255)
                                                                      gcgcATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCg                                                                           <  2:4282604/150‑1 (MQ=255)
                                                                      gcgcATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCg                                                                           <  1:2543287/150‑1 (MQ=255)
                                                                         cATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGc                                                                        >  2:311268/1‑150 (MQ=255)
                                                                          aTCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCa                                                                       >  1:2079222/1‑150 (MQ=255)
                                                                          aTCGGCATGTTGATTGAGATAGCCAACCATCGACGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCa                                                                       <  2:4367739/150‑1 (MQ=255)
                                                                            cGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGt                                                                     >  2:2331776/1‑150 (MQ=255)
                                                                             ggCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTg                                                                    >  2:314412/1‑150 (MQ=255)
                                                                              gCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGc                                                                   >  1:3390012/1‑150 (MQ=255)
                                                                              gCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGc                                                                   <  1:1292485/150‑1 (MQ=255)
                                                                               cATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCg                                                                  >  1:1019138/1‑150 (MQ=255)
                                                                                aTGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGa                                                                 >  2:2113644/1‑150 (MQ=255)
                                                                                aTGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGa                                                                 >  1:2902862/1‑150 (MQ=255)
                                                                                aTGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGa                                                                 >  1:1564263/1‑150 (MQ=255)
                                                                                 tGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGaa                                                                >  1:993886/1‑150 (MQ=255)
                                                                                  gTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAg                                                               <  1:3193215/150‑1 (MQ=255)
                                                                                   ttgattgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAgg                                                              >  1:3000088/1‑150 (MQ=255)
                                                                                        tgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACa                                                         >  2:4530033/1‑150 (MQ=255)
                                                                                          agaTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGc                                                       <  2:1711231/150‑1 (MQ=255)
                                                                                           gaTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGcc                                                      <  1:3745542/150‑1 (MQ=255)
                                                                                           gaTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGcc                                                      <  1:2271139/150‑1 (MQ=255)
                                                                                           gaTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGcc                                                      <  1:1735373/150‑1 (MQ=255)
                                                                                            aTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCg                                                     >  2:2985874/1‑150 (MQ=255)
                                                                                              aGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTa                                                   <  1:3000838/150‑1 (MQ=255)
                                                                                              aGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTGGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTa                                                   >  1:1898940/1‑150 (MQ=255)
                                                                                               gCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTaa                                                  <  1:429930/150‑1 (MQ=255)
                                                                                                 cAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACg                                                <  1:4213620/150‑1 (MQ=255)
                                                                                                   aCCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGcc                                              >  1:1923855/1‑150 (MQ=255)
                                                                                                    ccATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCATCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCa                                             <  2:4443450/150‑1 (MQ=255)
                                                                                                     cATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAg                                            >  2:2944670/1‑150 (MQ=255)
                                                                                                     cATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAg                                            <  2:128848/150‑1 (MQ=255)
                                                                                                      aTCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGc                                           >  2:1350259/1‑150 (MQ=255)
                                                                                                           cgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAat                                      >  2:4141546/1‑150 (MQ=255)
                                                                                                           cgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAAACAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGATCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAat                                      >  2:1402475/1‑150 (MQ=255)
                                                                                                           cgccAGCGTGGTAGATTAGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTAGGTTCGCCAAGCTGTTCTAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAat                                      >  2:40781/1‑150 (MQ=255)
                                                                                                             ccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgc                                    >  2:1084188/1‑150 (MQ=255)
                                                                                                              cAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgcc                                   >  2:888109/1‑150 (MQ=255)
                                                                                                              cAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgcc                                   >  2:2827269/1‑150 (MQ=255)
                                                                                                                gCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgccat                                 <  2:1898370/150‑1 (MQ=255)
                                                                                                                 cGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgccatg                                <  1:200791/150‑1 (MQ=255)
                                                                                                                  gTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgccatgc                               >  2:1814945/1‑150 (MQ=255)
                                                                                                                  gTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgccatgc                               >  2:2091514/1‑150 (MQ=255)
                                                                                                                    ggTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCg                             <  2:74398/150‑1 (MQ=255)
                                                                                                                        gATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgc                         >  2:6330/1‑150 (MQ=255)
                                                                                                                               cACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATg                  <  1:4513796/150‑1 (MQ=255)
                                                                                                                               cACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATg                  <  1:1003764/150‑1 (MQ=255)
                                                                                                                               cACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATg                  <  1:872991/150‑1 (MQ=255)
                                                                                                                                aCTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATgc                 >  1:1397322/1‑150 (MQ=255)
                                                                                                                                  tCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATgcgc               >  2:2420704/1‑150 (MQ=255)
                                                                                                                                   ccccGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCt              <  2:2715331/150‑1 (MQ=255)
                                                                                                                                      cGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGcc           >  1:146593/1‑150 (MQ=255)
                                                                                                                                         cGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACg        <  2:3047384/150‑1 (MQ=255)
                                                                                                                                         cGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACg        <  1:1595385/150‑1 (MQ=255)
                                                                                                                                         cGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACg        >  1:129997/1‑150 (MQ=255)
                                                                                                                                          gCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGc       >  2:2511820/1‑150 (MQ=255)
                                                                                                                                             cccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAAt    >  1:2064840/1‑150 (MQ=255)
                                                                                                                                             cccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAAt    >  1:306702/1‑150 (MQ=255)
                                                                                                                                              ccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATc   <  2:3799177/150‑1 (MQ=255)
                                                                                                                                              ccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATc   <  2:3322217/150‑1 (MQ=255)
                                                                                                                                              ccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATACCGGCAATACCGTTGGTTCGCCAAGCTGTTCGAGCTGCTGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCTAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATc   >  1:1986177/1‑150 (MQ=255)
                                                                                                                                               cGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCg  >  2:2012851/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCG  >  CP009273/3087919‑3088211

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: