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,0650CA99.7% 1067.4 / ‑4.6 330V125V (GTG→GTTyggRputative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (177/152);  total (177/153)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 330 total aligned reads displayed.

TTCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATC  >  CP009273/3087918‑3088210
                                                                                                                                                   |                                                                                                                                                 
ttCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAAc                                                                                                                                                 <  2:4433600/150‑1 (MQ=255)
  cccGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCa                                                                                                                                               <  2:1085794/150‑1 (MQ=255)
      ctgctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCTGAAt                                                                                                                                           <  1:2739420/150‑1 (MQ=255)
       tgctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATGGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATc                                                                                                                                          <  1:171136/150‑1 (MQ=255)
         ctgGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTACCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGACGCCCCCGTAACCAGAATCAg                                                                                                                                        >  1:937513/1‑150 (MQ=255)
           ggATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGc                                                                                                                                      <  2:3421091/150‑1 (MQ=255)
           ggATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGc                                                                                                                                      >  2:4042899/1‑150 (MQ=255)
           ggATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGc                                                                                                                                      >  2:576584/1‑150 (MQ=255)
              tCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCAt                                                                                                                                   >  1:2794287/1‑150 (MQ=255)
                     tCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGct                                                                                                                            <  1:319108/150‑1 (MQ=255)
                     tCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACCCCCCGTCGCCCCCGTAACCAGAAACAGGCCATTCTCGct                                                                                                                            <  1:674544/150‑1 (MQ=255)
                      cGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGctc                                                                                                                           <  2:209133/150‑1 (MQ=255)
                       gctggctgGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGTATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGctct                                                                                                                          <  2:4900452/150‑1 (MQ=255)
                        ctggctggCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCtt                                                                                                                         >  1:186124/1‑150 (MQ=255)
                         tggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGTATCAGGCCATTCTCGCTCTTg                                                                                                                        >  1:1933987/1‑150 (MQ=255)
                         tggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTg                                                                                                                        <  1:2270296/150‑1 (MQ=255)
                         tggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTg                                                                                                                        <  1:4915339/150‑1 (MQ=255)
                         tggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTg                                                                                                                        >  2:3118002/1‑150 (MQ=255)
                          ggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGa                                                                                                                       <  2:2573963/150‑1 (MQ=255)
                          ggctggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGa                                                                                                                       <  2:1780973/150‑1 (MQ=255)
                             tggcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTa                                                                                                                    <  2:3876119/150‑1 (MQ=255)
                             tggcATAGAGAGATTCAACCGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTa                                                                                                                    <  1:3452177/150‑1 (MQ=255)
                               gcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAAt                                                                                                                  <  1:5187011/150‑1 (MQ=255)
                               gcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAAt                                                                                                                  >  1:4041981/1‑150 (MQ=255)
                               gcATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAAt                                                                                                                  <  2:4214581/150‑1 (MQ=255)
                                cATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAAtt                                                                                                                 >  1:2988428/1‑150 (MQ=255)
                                cATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAAtt                                                                                                                 >  1:3234063/1‑150 (MQ=255)
                                cATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAAtt                                                                                                                 >  2:5251215/1‑150 (MQ=255)
                                 aTAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTc                                                                                                                >  1:1006048/1‑150 (MQ=255)
                                  tAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTcc                                                                                                               >  2:2464895/1‑150 (MQ=255)
                                      gaTATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCa                                                                                                           >  1:3784582/1‑150 (MQ=255)
                                       atatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCaa                                                                                                          >  1:3525518/1‑150 (MQ=255)
                                       atatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCaa                                                                                                          >  1:3166696/1‑150 (MQ=255)
                                        tatTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCATCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAAt                                                                                                         <  2:2703883/150‑1 (MQ=255)
                                         atTCATCAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGAGAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGTATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATa                                                                                                        <  2:4835318/150‑1 (MQ=255)
                                                aGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTgg                                                                                                 >  2:1916933/1‑150 (MQ=255)
                                                aGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTgg                                                                                                 >  2:959188/1‑150 (MQ=255)
                                                aGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGAGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGATCTTGAGTAATTCCGGCAATACCGTTgg                                                                                                 >  2:4630067/1‑150 (MQ=255)
                                                 ggATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGt                                                                                                >  2:1288870/1‑150 (MQ=255)
                                                   aTCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgc                                                                                              <  2:1418753/150‑1 (MQ=255)
                                                    tCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcg                                                                                             >  1:5003815/1‑150 (MQ=255)
                                                    tCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTgcg                                                                                             >  1:121332/1‑150 (MQ=255)
                                                         cAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAg                                                                                        <  2:1818195/150‑1 (MQ=255)
                                                         cAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAg                                                                                        >  2:3165308/1‑150 (MQ=255)
                                                         cAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAg                                                                                        <  2:3629735/150‑1 (MQ=255)
                                                          aGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGc                                                                                       >  1:4334149/1‑150 (MQ=255)
                                                           gCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCt                                                                                      <  1:4131696/150‑1 (MQ=255)
                                                             gTCAGAAGATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGt                                                                                    <  1:3916107/150‑1 (MQ=255)
                                                               cAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTc                                                                                  <  2:5040032/150‑1 (MQ=255)
                                                                 gAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGa                                                                                >  2:3480360/1‑150 (MQ=255)
                                                                 gAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGa                                                                                <  1:3117394/150‑1 (MQ=255)
                                                                  aaTATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAg                                                                               <  1:180659/150‑1 (MQ=255)
                                                                   atatGCGCATCGGCATGTTTATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGAATTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGc                                                                              >  1:4693712/1‑150 (MQ=255)
                                                                       gcgcATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCg                                                                          <  2:3004861/150‑1 (MQ=255)
                                                                              ggCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTg                                                                   >  1:4512890/1‑150 (MQ=255)
                                                                                cATGTTGATTGAGATAGCCAAACATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCg                                                                 <  1:1353299/150‑1 (MQ=255)
                                                                                 aTGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGa                                                                >  2:3604096/1‑150 (MQ=255)
                                                                                 aTGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGa                                                                >  2:4689062/1‑150 (MQ=255)
                                                                                   gTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAg                                                              <  2:4124029/150‑1 (MQ=255)
                                                                                   gTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAg                                                              >  2:1936604/1‑150 (MQ=255)
                                                                                    ttgattgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAgg                                                             >  1:3360624/1‑150 (MQ=255)
                                                                                       attgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTaa                                                          >  2:2873185/1‑150 (MQ=255)
                                                                                        ttgaGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAAc                                                         <  2:2794287/150‑1 (MQ=255)
                                                                                          gagaTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAg                                                       >  2:789683/1‑150 (MQ=255)
                                                                                            gaTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGCAGGTAACAGcc                                                     >  1:330285/1‑150 (MQ=255)
                                                                                            gaTAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGcc                                                     >  1:4190534/1‑150 (MQ=255)
                                                                                              tAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGt                                                   >  1:1155171/1‑150 (MQ=255)
                                                                                               aGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTa                                                  >  2:853336/1‑150 (MQ=255)
                                                                                               aGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTa                                                  <  2:3342798/150‑1 (MQ=255)
                                                                                                 ccAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAAc                                                <  2:2330988/150‑1 (MQ=255)
                                                                                                 ccAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAAc                                                <  1:481510/150‑1 (MQ=255)
                                                                                                   aaCCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGc                                              <  2:2386981/150‑1 (MQ=255)
                                                                                                   aaCCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGc                                              <  2:2210315/150‑1 (MQ=255)
                                                                                                    aCCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGcc                                             >  1:5151528/1‑150 (MQ=255)
                                                                                                     ccATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCa                                            <  1:1898011/150‑1 (MQ=255)
                                                                                                       aTCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGc                                          >  2:3274856/1‑150 (MQ=255)
                                                                                                       aTCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGc                                          <  1:1542452/150‑1 (MQ=255)
                                                                                                         cgccgcCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGa                                        <  1:4430292/150‑1 (MQ=255)
                                                                                                            cgccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAat                                     <  1:1478984/150‑1 (MQ=255)
                                                                                                             gccAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatg                                    <  1:1070607/150‑1 (MQ=255)
                                                                                                                aGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgcca                                 >  2:1601326/1‑150 (MQ=255)
                                                                                                                    tGGTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAatgccatgcc                             >  1:2436251/1‑150 (MQ=255)
                                                                                                                     ggTAGATTTGCCACTCCCCGTCGCCCCCGTAGCCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCg                            >  2:4361133/1‑150 (MQ=255)
                                                                                                                     ggTAGATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCg                            >  1:3737309/1‑150 (MQ=255)
                                                                                                                         gATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgc                        >  2:5315476/1‑150 (MQ=255)
                                                                                                                         gATTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgc                        >  1:259971/1‑150 (MQ=255)
                                                                                                                          aTTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgcg                       <  1:545738/150‑1 (MQ=255)
                                                                                                                          aTTTGCCACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTgcg                       >  1:4738186/1‑150 (MQ=255)
                                                                                                                               ccACTCCCCGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAAt                  <  2:5100161/150‑1 (MQ=255)
                                                                                                                                    ccccGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCt             <  2:2988428/150‑1 (MQ=255)
                                                                                                                                    ccccGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCt             >  2:3778619/1‑150 (MQ=255)
                                                                                                                                    caccGACGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTTCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACTCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCt             >  1:2925476/3‑150 (MQ=255)
                                                                                                                                     cccGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTg            >  1:815859/1‑150 (MQ=255)
                                                                                                                                     cccGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTg            >  1:2030608/1‑150 (MQ=255)
                                                                                                                                       cGTCGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGcc          <  1:178896/150‑1 (MQ=255)
                                                                                                                                          cGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACg       <  2:1013218/150‑1 (MQ=255)
                                                                                                                                          cGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACg       <  1:212591/150‑1 (MQ=255)
                                                                                                                                          cGCCCCCGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACg       >  1:5067378/1‑150 (MQ=255)
                                                                                                                                             ccccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCaa    >  1:3961721/1‑150 (MQ=255)
                                                                                                                                               ccGTAACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATc  <  2:4645852/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTCCCGCTGCTGGATCAAACATCGCTGGCTGGCATAGAGATATTCAACAGGATCTTCCAGCGTCAGAATATGCGCATCGGCATGTTGATTGAGATAGCCAACCATCGCCGCCAGCGTGGTAGATTTGCCACTCCCCGTCGCCCCCGTCACCAGAATCAGGCCATTCTCGCTCTTGAGTAATTCCGGCAATACCGTTGGTGCGCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATC  >  CP009273/3087918‑3088210

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: