Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,895,069 C→T 13.5% G258G (GGG→GGA yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,895,0690CT13.5% 750.7 / 83.2 254G258G (GGG→GGAyoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base C (112/105);  new base T (15/19);  total (129/125)
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.23e-01

Only 100 of 262 total aligned reads displayed.

CCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAAT  >  CP009273/1894920‑1895218
                                                                                                                                                     |                                                                                                                                                     
 cgcgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTcc                                                                                                                                                      <  2:1750411/150‑1 (MQ=255)
  gcgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTtcc                                                                                                                                                     <  2:1539026/150‑1 (MQ=255)
  gcgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTccc                                                                                                                                                     <  2:1296532/150‑1 (MQ=255)
      ggcgTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGccc                                                                                                                                                 <  1:1163899/150‑1 (MQ=255)
        agTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGTCATCGTGCTGAACATTGGCCTGACGTTTACCGCCCAt                                                                                                                                               <  1:1284983/149‑1 (MQ=255)
         gTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATc                                                                                                                                              >  1:4481994/1‑150 (MQ=255)
         gTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATc                                                                                                                                              <  2:1688079/150‑1 (MQ=255)
          tCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCa                                                                                                                                             <  1:3071722/150‑1 (MQ=255)
           cATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCaa                                                                                                                                            <  2:942942/150‑1 (MQ=255)
             tGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAAc                                                                                                                                          <  2:3678323/150‑1 (MQ=255)
               aTCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGc                                                                                                                                        >  1:284578/1‑150 (MQ=255)
                tCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCa                                                                                                                                       >  2:1616154/1‑150 (MQ=255)
                tCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCa                                                                                                                                       >  1:2461902/1‑150 (MQ=255)
                  ccGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAgg                                                                                                                                     <  2:341507/150‑1 (MQ=255)
                     cgcAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATg                                                                                                                                  <  2:3239904/150‑1 (MQ=255)
                      gcAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATgg                                                                                                                                 >  2:3532409/1‑150 (MQ=255)
                        agagAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCa                                                                                                                               >  1:1628709/1‑150 (MQ=255)
                         gagaACGCGACGCCAGCGTCAGTACGCCGTTAATCAGGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCCTCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCAt                                                                                                                              >  1:992250/1‑150 (MQ=255)
                           gaACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCg                                                                                                                            >  2:3324144/1‑150 (MQ=255)
                           gaACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCg                                                                                                                            <  2:1169929/150‑1 (MQ=255)
                           gaACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCg                                                                                                                            <  1:2769294/150‑1 (MQ=255)
                               gcgACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAg                                                                                                                        <  1:108939/150‑1 (MQ=255)
                                cgACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGt                                                                                                                       <  1:1315408/150‑1 (MQ=255)
                                  aCGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTAc                                                                                                                     <  1:2037293/150‑1 (MQ=255)
                                     ccAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTAcgcg                                                                                                                  >  2:3270490/1‑150 (MQ=255)
                                        gCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTc                                                                                                               <  2:356078/150‑1 (MQ=255)
                                         cGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCg                                                                                                              <  1:2703738/150‑1 (MQ=255)
                                          gTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGc                                                                                                             >  2:1411641/1‑150 (MQ=255)
                                                 cGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCa                                                                                                      >  1:2383930/1‑150 (MQ=255)
                                                    cGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAg                                                                                                   <  2:3096794/150‑1 (MQ=255)
                                                    cGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAg                                                                                                   <  2:4488313/150‑1 (MQ=255)
                                                    cGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAg                                                                                                   <  1:4366650/150‑1 (MQ=255)
                                                       tAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCg                                                                                                >  2:4312581/1‑150 (MQ=255)
                                                       tAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCg                                                                                                >  2:4532261/1‑150 (MQ=255)
                                                          tCATGTAACGTGCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACt                                                                                             <  2:2929246/150‑1 (MQ=255)
                                                           cATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACtg                                                                                            <  2:3094154/150‑1 (MQ=255)
                                                           cATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACtg                                                                                            >  2:2991701/1‑150 (MQ=255)
                                                            aTGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGTGTTATCGGCATCGTGCTTAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACtgg                                                                                           <  1:3762611/150‑1 (MQ=255)
                                                              gTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACtggtg                                                                                         <  2:2249620/150‑1 (MQ=255)
                                                               tAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACtggtgg                                                                                        <  1:2555850/150‑1 (MQ=255)
                                                                aaCGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGc                                                                                       >  2:1385925/1‑150 (MQ=255)
                                                                aaCGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGc                                                                                       <  1:1907197/150‑1 (MQ=255)
                                                                aaCGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGc                                                                                       <  1:2361625/150‑1 (MQ=255)
                                                                  cGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGa                                                                                     <  1:2653102/150‑1 (MQ=255)
                                                                     tcttcttcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATa                                                                                  <  2:4159099/150‑1 (MQ=255)
                                                                       ttcttcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATaaa                                                                                <  1:455836/150‑1 (MQ=255)
                                                                       ttcttcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATaaa                                                                                >  2:3809317/1‑150 (MQ=255)
                                                                           tcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTg                                                                            <  2:2342602/150‑1 (MQ=255)
                                                                           tcGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTg                                                                            <  2:1568606/150‑1 (MQ=255)
                                                                             ggCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTACGGTTACCGCCCATCAAACGCAGGATGGCCTCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACGGGTGGCGAATAAAGTTGCGa                                                                         >  2:2378060/1‑150 (MQ=255)
                                                                              gCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGa                                                                         >  2:3592195/1‑150 (MQ=255)
                                                                                aaaTGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGAcg                                                                       >  2:204784/1‑150 (MQ=255)
                                                                                aaaTGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGAcg                                                                       <  1:3621000/150‑1 (MQ=255)
                                                                                aaaTGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGAcg                                                                       >  2:1075570/1‑150 (MQ=255)
                                                                                    gCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCa                                                                   <  2:724659/150‑1 (MQ=255)
                                                                                     cACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCaa                                                                  >  1:2898983/1‑150 (MQ=255)
                                                                                      aCCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAAt                                                                 <  1:2338260/150‑1 (MQ=255)
                                                                                          tCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTgg                                                             <  2:4187146/150‑1 (MQ=255)
                                                                                           ccGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGt                                                            <  2:2090871/150‑1 (MQ=255)
                                                                                             gngATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTNCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCNAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTNGTTa                                                          >  2:480365/1‑150 (MQ=255)
                                                                                             gggATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTa                                                          >  1:3089804/1‑150 (MQ=255)
                                                                                             gggATCGGCATCGGCATCGGGGTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTAAGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTa                                                          >  1:2211653/1‑150 (MQ=255)
                                                                                                 tcggcatcggcaTCGGGTTATTGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAaca                                                      >  2:1748137/1‑150 (MQ=255)
                                                                                                   ggcatcggcatcGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACt                                                    >  2:4062683/1‑150 (MQ=255)
                                                                                                      atcggcatcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCg                                                 >  1:712267/1‑150 (MQ=255)
                                                                                                       tcggcatcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCga                                                <  1:1359710/150‑1 (MQ=255)
                                                                                                        cggcatcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgat                                               <  2:68855/150‑1 (MQ=255)
                                                                                                        cggcatcggGTTATCGGCATCGTGCTGAACATTGGACTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgat                                               <  1:183579/150‑1 (MQ=255)
                                                                                                         ggcatcggGTTATCGGCATCGTGCGGAACATTGCCCTGACGTTTAACCCACATCAAACGCAGGATGGCATCGGAAGTACGCGCTCGCAGCGGCAAAGTCGACTGGAGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgatg                                              >  2:78377/1‑150 (MQ=255)
                                                                                                          gcatcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgatga                                             <  1:57109/150‑1 (MQ=255)
                                                                                                           catcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgatgat                                            <  1:4605981/150‑1 (MQ=255)
                                                                                                            atcggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgatgatg                                           >  1:613860/1‑150 (MQ=255)
                                                                                                              cggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgatgatgat                                         >  1:1250254/1‑150 (MQ=255)
                                                                                                              cggGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCgatgatgat                                         <  1:1890746/150‑1 (MQ=255)
                                                                                                                ggTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCg                                       <  2:2540161/150‑1 (MQ=255)
                                                                                                                   tATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCgaga                                    >  1:764271/1‑150 (MQ=255)
                                                                                                                    aTCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGaa                                   >  1:3936322/1‑150 (MQ=255)
                                                                                                                       ggCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGcc                                >  1:1163296/1‑150 (MQ=255)
                                                                                                                        gCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCa                               <  2:4221265/150‑1 (MQ=255)
                                                                                                                         cATCGTGCTGAACATTGGCCTGACGTTTTCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCaa                              >  2:52588/1‑150 (MQ=255)
                                                                                                                         cATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCaa                              >  2:4340070/1‑150 (MQ=255)
                                                                                                                          aTCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAAt                             <  2:1789225/150‑1 (MQ=255)
                                                                                                                           tCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATc                            <  2:1101419/150‑1 (MQ=255)
                                                                                                                           tCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATc                            <  1:517698/150‑1 (MQ=255)
                                                                                                                            cGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATcg                           <  2:3874895/150‑1 (MQ=255)
                                                                                                                             gTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATcgc                          >  2:2734369/1‑150 (MQ=255)
                                                                                                                              tGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATcgcg                         >  2:1867749/1‑150 (MQ=255)
                                                                                                                                cTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGc                       >  1:4316006/1‑150 (MQ=255)
                                                                                                                                 tGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCa                      <  2:4452877/150‑1 (MQ=255)
                                                                                                                                      aTTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAg                 >  2:2286185/1‑150 (MQ=255)
                                                                                                                                         ggCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTa              <  1:842679/150‑1 (MQ=255)
                                                                                                                                           ccTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAAc            >  2:286194/1‑150 (MQ=255)
                                                                                                                                            cTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAAcc           <  2:1666327/150‑1 (MQ=255)
                                                                                                                                               aCGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCttt        >  1:1425916/1‑150 (MQ=255)
                                                                                                                                               aCGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCttt        >  1:1888061/1‑150 (MQ=255)
                                                                                                                                               aCGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCttt        <  1:1011203/150‑1 (MQ=255)
                                                                                                                                                  tttCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCgg     <  1:3185160/150‑1 (MQ=255)
                                                                                                                                                   ttCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGa    <  1:2060771/150‑1 (MQ=255)
                                                                                                                                                    tCCCCCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGaa   >  2:2527758/1‑150 (MQ=255)
                                                                                                                                                     cccGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAAt  <  1:1752522/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCAGGATGGCATCGGCAGTACGCGCTCGCAGCGGCAAAGTCGACTGGTGGCGAATAAAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAAT  >  CP009273/1894920‑1895218

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: