Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,088,065 C→A 100% 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,0650CA98.8% 1100.0 / 2.2 338V125V (GTG→GTTyggRputative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (158/176);  total (161/177)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
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: Variant not supported by required number of reads on each strand.

Only 100 of 338 total aligned reads displayed.

TTCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCG  >  CP009273/3087918‑3088211
                                                                                                                                                   |                                                                                                                                                  
ttCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGATATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCAACTCCCCGTCGCCCCCGTAAc                                                                                                                                                  >  2:1979729/1‑150 (MQ=255)
  cccGCTGCTGGATCAAACATCGCTGGCTGGCATAGATATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCa                                                                                                                                                <  1:4298855/150‑1 (MQ=255)
  cccGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCa                                                                                                                                                >  1:812102/1‑150 (MQ=255)
  cccGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCa                                                                                                                                                >  1:3347938/1‑150 (MQ=255)
   ccGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCat                                                                                                                                               <  2:2139084/150‑2 (MQ=255)
   ccGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAg                                                                                                                                               <  2:2468741/150‑1 (MQ=255)
        gctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCa                                                                                                                                          >  1:3715664/1‑150 (MQ=255)
         ctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAg                                                                                                                                         >  2:3609209/1‑150 (MQ=255)
          tgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAgg                                                                                                                                        >  2:4024210/1‑150 (MQ=255)
            gATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGcc                                                                                                                                      <  1:1288032/150‑1 (MQ=255)
            gATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGcc                                                                                                                                      >  2:2520758/1‑150 (MQ=255)
            gATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGcc                                                                                                                                      >  2:5078889/1‑150 (MQ=255)
              tCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCAt                                                                                                                                    <  2:4840397/150‑1 (MQ=255)
                aaaCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATtc                                                                                                                                  >  1:3257483/1‑150 (MQ=255)
                 aaCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATtct                                                                                                                                 <  2:1874735/150‑1 (MQ=255)
                  aCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATtctc                                                                                                                                <  2:3217063/150‑1 (MQ=255)
                    aTCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGc                                                                                                                              >  2:4066815/1‑150 (MQ=255)
                    aTCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGc                                                                                                                              >  1:3116001/1‑150 (MQ=255)
                     tCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGct                                                                                                                             >  2:352333/1‑150 (MQ=255)
                     tCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGct                                                                                                                             <  2:969376/150‑1 (MQ=255)
                       gctggctgGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGctct                                                                                                                           <  2:748034/150‑1 (MQ=255)
                        ctggctggCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCtt                                                                                                                          <  2:1318216/150‑1 (MQ=255)
                        ctggctggCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCtt                                                                                                                          <  2:4288848/150‑1 (MQ=255)
                                cATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAAtt                                                                                                                  <  2:3457155/150‑1 (MQ=255)
                                  tAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGACTCAGGCCATTCTCGCTCTTGAGTAATTcc                                                                                                                >  1:3536726/1‑150 (MQ=255)
                                        tatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAAt                                                                                                          <  1:386963/150‑1 (MQ=255)
                                        tatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAAt                                                                                                          <  1:173250/150‑1 (MQ=255)
                                               cAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTg                                                                                                   <  2:867483/150‑1 (MQ=255)
                                                    tCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcg                                                                                              <  1:990180/150‑1 (MQ=255)
                                                     cTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcgc                                                                                             >  1:4795254/1‑150 (MQ=255)
                                                     cTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcgc                                                                                             >  1:1338576/1‑150 (MQ=255)
                                                      ttCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGcc                                                                                            <  2:1585740/150‑1 (MQ=255)
                                                      ttCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCAGTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGcc                                                                                            <  2:4993920/150‑1 (MQ=255)
                                                       tCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCa                                                                                           <  1:2012131/150‑1 (MQ=255)
                                                       tCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCa                                                                                           >  2:419548/1‑150 (MQ=255)
                                                         cAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAg                                                                                         <  1:2520758/150‑1 (MQ=255)
                                                             gTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGt                                                                                     <  2:2021245/150‑1 (MQ=255)
                                                               cAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTc                                                                                   <  1:2702504/150‑1 (MQ=255)
                                                               cAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAAGTCCGGCAATACCGTTGGTGCGCCAAGCTGTTc                                                                                   <  2:5128502/150‑1 (MQ=255)
                                                                aGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCg                                                                                  <  2:3831182/150‑1 (MQ=255)
                                                                  aaTATGCGCATCGGCATGTTGATTGAGATAGCCAACCCTCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAg                                                                                <  2:3430379/150‑1 (MQ=255)
                                                                   atatGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGc                                                                               >  1:2298425/1‑150 (MQ=255)
                                                                    tatGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCt                                                                              >  2:4780897/1‑150 (MQ=255)
                                                                    tatGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCt                                                                              <  2:2057103/150‑1 (MQ=255)
                                                                     atGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTg                                                                             >  1:2741849/1‑150 (MQ=255)
                                                                      tGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGc                                                                            <  1:4899991/150‑1 (MQ=255)
                                                                      tGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGc                                                                            <  1:4841651/150‑1 (MQ=255)
                                                                       gcgcATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCg                                                                           <  2:2784312/150‑1 (MQ=255)
                                                                       gcgcATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCg                                                                           <  1:2776690/150‑1 (MQ=255)
                                                                          cATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGc                                                                        >  2:308090/1‑150 (MQ=255)
                                                                          cATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGc                                                                        <  1:2485844/150‑1 (MQ=255)
                                                                           aTCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCa                                                                       >  1:4184671/1‑150 (MQ=255)
                                                                             cGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGt                                                                     >  2:2646852/1‑150 (MQ=255)
                                                                             cGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGt                                                                     >  2:683743/1‑150 (MQ=255)
                                                                             cGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGt                                                                     >  2:4503432/1‑150 (MQ=255)
                                                                             cGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGt                                                                     >  2:4662825/1‑150 (MQ=255)
                                                                                cATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCg                                                                  >  2:2168281/1‑150 (MQ=255)
                                                                                cATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCg                                                                  <  1:1094458/150‑1 (MQ=255)
                                                                                 aTGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGa                                                                 <  1:413075/150‑1 (MQ=255)
                                                                                   gTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAg                                                               >  1:4571787/1‑150 (MQ=255)
                                                                                    ttgattgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAgg                                                              >  1:2958309/1‑150 (MQ=255)
                                                                                    ttgattgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAgg                                                              >  2:2245012/1‑150 (MQ=255)
                                                                                    ttgattgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTAGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAgg                                                              <  1:4456829/150‑1 (MQ=255)
                                                                                      gattgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTa                                                            >  2:2342493/1‑150 (MQ=255)
                                                                                         tgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACa                                                         <  2:1061733/150‑1 (MQ=255)
                                                                                            gaTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGcc                                                      >  1:1181055/1‑150 (MQ=255)
                                                                                            gaTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCCGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGcc                                                      <  1:436119/150‑1 (MQ=255)
                                                                                              tAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGt                                                    <  1:2826005/150‑1 (MQ=255)
                                                                                               aGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTa                                                   >  1:902820/1‑150 (MQ=255)
                                                                                                  cAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACg                                                <  1:507017/150‑1 (MQ=255)
                                                                                                     ccATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCa                                             >  1:1802280/1‑150 (MQ=255)
                                                                                                     ccATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCa                                             <  2:587673/150‑1 (MQ=255)
                                                                                                       aTCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTGGAGTAATTCCGGCCATACCGTTGGTGCGCCAAGCTGTGCGCGCTGCGGGCAGGGCGACGGTACCAGCCGTCACGCCAGc                                           >  2:4356731/1‑150 (MQ=255)
                                                                                                        tCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCg                                          >  1:3000324/1‑150 (MQ=255)
                                                                                                         cgccgcCAGCGTGGTAGATGTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCCGCGa                                         <  2:2708969/150‑1 (MQ=255)
                                                                                                          gccgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAATAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGaa                                        >  2:4751243/1‑150 (MQ=255)
                                                                                                            cgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTACCCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAat                                      >  1:5085250/1‑150 (MQ=255)
                                                                                                            cgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAat                                      <  1:4234112/150‑1 (MQ=255)
                                                                                                            cgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAat                                      <  1:1001722/150‑1 (MQ=255)
                                                                                                                aGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgcca                                  <  2:1474754/150‑1 (MQ=255)
                                                                                                                aGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgcca                                  >  2:5054421/1‑150 (MQ=255)
                                                                                                                   gTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgccatgc                               <  1:3087058/150‑1 (MQ=255)
                                                                                                                       tAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGtt                           >  1:676927/1‑150 (MQ=255)
                                                                                                                        aGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTg                          <  2:2744099/150‑1 (MQ=255)
                                                                                                                         gATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgc                         <  1:1611041/150‑1 (MQ=255)
                                                                                                                         gATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgc                         <  1:1870253/150‑1 (MQ=255)
                                                                                                                         gATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgc                         <  1:587147/150‑1 (MQ=255)
                                                                                                                         gATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgc                         <  1:2600557/150‑1 (MQ=255)
                                                                                                                          atgtGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAAGTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgcg                        <  2:4046234/147‑1 (MQ=255)
                                                                                                                               ccACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAAt                   >  2:1720044/1‑150 (MQ=255)
                                                                                                                                cACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATg                  >  1:1927534/1‑150 (MQ=255)
                                                                                                                                 aCTCCCCGTCGCCCCCGTAACCAGAATCCGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATgc                 <  2:2335014/150‑1 (MQ=255)
                                                                                                                                 aCTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATgc                 <  1:315967/150‑1 (MQ=255)
                                                                                                                                   tCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATgcgc               <  2:4574061/150‑1 (MQ=255)
                                                                                                                                   tCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATgcgc               <  1:4526227/150‑1 (MQ=255)
                                                                                                                                     cccGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTg             <  1:1502830/150‑1 (MQ=255)
                                                                                                                                             ccccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCaa     <  1:232726/150‑1 (MQ=255)
                                                                                                                                              cccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAAt    <  1:4970990/150‑1 (MQ=255)
                                                                                                                                                cGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCg  >  1:597465/1‑150 (MQ=255)
                                                                                                                                                cGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCg  >  1:3583812/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCG  >  CP009273/3087918‑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: