Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,088,065 C→A 37.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,0650CA37.6% 413.6 / inf 302V125V (GTG→GTTyggRputative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base C (105/83);  new base A (56/58);  total (161/141)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

ATTTCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTG  >  CP009273/3087916‑3088214
                                                                                                                                                     |                                                                                                                                                     
aTTTCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTc                                                                                                                                                       <  1:4978360/150‑1 (MQ=255)
  ttCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAAc                                                                                                                                                     >  2:78939/1‑150 (MQ=255)
    cccGCTGCTGGATCACACATCGCGGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCa                                                                                                                                                   <  2:1116389/150‑1 (MQ=255)
      cGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGa                                                                                                                                                 >  2:931888/1‑150 (MQ=255)
      cGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGa                                                                                                                                                 <  1:2530989/150‑1 (MQ=255)
       gctgctGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGaa                                                                                                                                                >  2:2531788/1‑150 (MQ=255)
       gctgctGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGaa                                                                                                                                                >  1:4652930/1‑150 (MQ=255)
          gctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCa                                                                                                                                             <  1:608494/150‑1 (MQ=255)
           ctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAg                                                                                                                                            >  1:2858150/1‑150 (MQ=255)
           ctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAg                                                                                                                                            <  1:1022681/150‑1 (MQ=255)
            tgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAgg                                                                                                                                           <  1:3147821/150‑1 (MQ=255)
               aTCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCa                                                                                                                                        >  2:4314005/1‑150 (MQ=255)
               aTCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCa                                                                                                                                        <  1:3549530/150‑1 (MQ=255)
                 cAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCAtt                                                                                                                                      <  1:769486/150‑1 (MQ=255)
                  aaaCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATtc                                                                                                                                     <  1:4544493/150‑1 (MQ=255)
                    aCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATtctc                                                                                                                                   <  1:3487342/150‑1 (MQ=255)
                    aCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATtctc                                                                                                                                   <  1:1448890/150‑1 (MQ=255)
                    aCATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATtctc                                                                                                                                   <  1:1448998/150‑1 (MQ=255)
                     cATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCg                                                                                                                                  >  1:2964636/1‑150 (MQ=255)
                     cATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCg                                                                                                                                  <  2:2483323/150‑1 (MQ=255)
                       tCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGct                                                                                                                                <  1:3928187/150‑1 (MQ=255)
                       tCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGct                                                                                                                                >  2:3180287/1‑150 (MQ=255)
                        cGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCAAGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGctc                                                                                                                               >  1:4268815/1‑150 (MQ=255)
                            ggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGa                                                                                                                           <  2:2758072/150‑1 (MQ=255)
                            ggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGa                                                                                                                           <  1:3268282/150‑1 (MQ=255)
                             gctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAg                                                                                                                          >  1:2159154/1‑150 (MQ=255)
                                  cATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCTCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAAtt                                                                                                                     <  1:2613083/150‑1 (MQ=255)
                                   aTAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTc                                                                                                                    >  2:3041700/1‑150 (MQ=255)
                                     agagATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCg                                                                                                                  <  1:4464044/150‑1 (MQ=255)
                                     agagATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCg                                                                                                                  <  1:1081571/150‑1 (MQ=255)
                                     agagATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCg                                                                                                                  <  2:1432185/150‑1 (MQ=255)
                                          tatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAAt                                                                                                             >  1:4347176/1‑150 (MQ=255)
                                           atTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATa                                                                                                            >  2:4479222/1‑150 (MQ=255)
                                           atTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATa                                                                                                            <  2:2055180/150‑1 (MQ=255)
                                             tCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATaac                                                                                                          >  2:3243582/1‑148 (MQ=255)
                                             tCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATAcc                                                                                                          <  1:3642895/150‑1 (MQ=255)
                                               aaCAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGt                                                                                                        >  2:3577040/1‑150 (MQ=255)
                                               aaCAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGt                                                                                                        >  2:843071/1‑150 (MQ=255)
                                                aCAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGtt                                                                                                       <  2:4718663/150‑1 (MQ=255)
                                                 cAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTg                                                                                                      >  1:565310/1‑150 (MQ=255)
                                                    gATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTg                                                                                                   >  1:328037/1‑150 (MQ=255)
                                                      tCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcg                                                                                                 >  2:688588/1‑150 (MQ=255)
                                                      tCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcg                                                                                                 <  1:5046020/150‑1 (MQ=255)
                                                      tCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcg                                                                                                 >  1:4947117/1‑150 (MQ=255)
                                                        ttCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGcc                                                                                               >  1:1137685/1‑150 (MQ=255)
                                                            aGCGTCAGAATATGCGCATCGGCATTTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCACCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATAACGTTGGTGCGCCaatc                                                                                           >  1:135117/1‑148 (MQ=255)
                                                            aGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGc                                                                                           <  1:231706/150‑1 (MQ=255)
                                                             gCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCt                                                                                          >  2:2290897/1‑150 (MQ=255)
                                                              cGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTg                                                                                         >  2:1922240/1‑150 (MQ=255)
                                                                 cAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTc                                                                                      <  1:2432344/150‑1 (MQ=255)
                                                                  aGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCg                                                                                     >  1:1964036/1‑150 (MQ=255)
                                                                    aaTATGCGCATCGGCAGGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAg                                                                                   >  2:3858594/1‑150 (MQ=255)
                                                                     atatGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTAGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGc                                                                                  >  1:2551440/1‑150 (MQ=255)
                                                                      tatGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCt                                                                                 <  2:523050/150‑1 (MQ=255)
                                                                       atGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTg                                                                                >  1:4604737/1‑150 (MQ=255)
                                                                        tGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGc                                                                               <  1:1255478/150‑1 (MQ=255)
                                                                        tGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGc                                                                               >  2:1692994/1‑150 (MQ=255)
                                                                           gcATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCggg                                                                            >  2:379161/1‑150 (MQ=255)
                                                                              tCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAg                                                                         <  2:2875893/150‑1 (MQ=255)
                                                                               cGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGt                                                                        >  1:205661/1‑150 (MQ=255)
                                                                                 gCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGc                                                                      <  2:3910775/150‑1 (MQ=255)
                                                                                  cATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCg                                                                     >  1:1572410/1‑150 (MQ=255)
                                                                                     gTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAg                                                                  <  1:4757891/150‑1 (MQ=255)
                                                                                         attgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTaa                                                              >  1:4080600/1‑150 (MQ=255)
                                                                                         attgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTaa                                                              >  2:5127374/1‑150 (MQ=255)
                                                                                          ttgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAAc                                                             >  2:2622200/1‑150 (MQ=255)
                                                                                              gaTAGCCTACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCACTACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGcc                                                         <  1:1541236/150‑1 (MQ=255)
                                                                                                tAGCAAACCATCGACGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCAGGCAATACCGTAGGTCCGCCAAGCTGTTCGAGCTTCGGGCAGTGCGAAGGTAACAGCCGt                                                       >  1:2186480/1‑150 (MQ=255)
                                                                                                      aCCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGcc                                                 <  2:746760/150‑1 (MQ=255)
                                                                                                           cgccgcCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGa                                            >  2:3701620/1‑150 (MQ=255)
                                                                                                            gccgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGaa                                           >  1:4940746/1‑150 (MQ=255)
                                                                                                            gccgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGaa                                           >  1:2236372/1‑150 (MQ=255)
                                                                                                              cgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAat                                         >  1:1845231/1‑150 (MQ=255)
                                                                                                               gccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGACatg                                        >  2:1654079/1‑150 (MQ=255)
                                                                                                                ccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgc                                       >  1:897078/1‑150 (MQ=255)
                                                                                                                 cAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgcc                                      >  2:5006195/1‑150 (MQ=255)
                                                                                                                    cGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgccatg                                   <  2:41509/150‑1 (MQ=255)
                                                                                                                    cGTGGGAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCAGTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgccatg                                   <  2:4825445/150‑1 (MQ=255)
                                                                                                                          aGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTg                             >  2:3193114/1‑150 (MQ=255)
                                                                                                                          aGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTg                             <  1:3419867/150‑1 (MQ=255)
                                                                                                                              ttGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGcgcg                         <  2:1310910/150‑1 (MQ=255)
                                                                                                                              ttGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGcgcg                         <  2:1631713/150‑1 (MQ=255)
                                                                                                                                gCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGaa                       <  2:4250376/150‑1 (MQ=255)
                                                                                                                                 ccACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAAt                      <  1:377237/150‑1 (MQ=255)
                                                                                                                                  cACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATg                     <  1:1650000/150‑1 (MQ=255)
                                                                                                                                  cACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTCACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATg                     >  2:3514434/1‑150 (MQ=255)
                                                                                                                                     tCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATgcgc                  >  1:807223/1‑150 (MQ=255)
                                                                                                                                     tCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATgcgc                  <  1:322684/150‑1 (MQ=255)
                                                                                                                                     tCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATgcgc                  >  2:4017648/1‑150 (MQ=255)
                                                                                                                                      ccccGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCt                 <  2:1796204/150‑1 (MQ=255)
                                                                                                                                       cccGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTg                <  2:329976/150‑1 (MQ=255)
                                                                                                                                        ccGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGc               <  1:4237768/150‑1 (MQ=255)
                                                                                                                                         cGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGcc              <  2:2985423/150‑1 (MQ=255)
                                                                                                                                           tCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCGGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCAc            <  2:2547286/150‑1 (MQ=255)
                                                                                                                                            cGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACg           >  2:2130236/1‑150 (MQ=255)
                                                                                                                                              cccccGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCa         <  1:3063288/150‑1 (MQ=255)
                                                                                                                                               ccccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCaa        >  1:2384877/1‑150 (MQ=255)
                                                                                                                                                cccGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAAt       <  2:3434736/150‑1 (MQ=255)
                                                                                                                                                     cACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCCGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTg  >  2:4515303/1‑150 (MQ=255)
                                                                                                                                                     cACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTg  <  1:140074/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTTCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTG  >  CP009273/3087916‑3088214

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: