Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 571,745 G→T 43.2% pseudogene (537/1040 nt) nmpC ← DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273571,7450GT43.2% 251.3 / inf 257pseudogene (537/1040 nt)nmpCDLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base G (78/68);  new base T (63/48);  total (141/116)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 261 total aligned reads displayed.

GAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGCAGGACGTC  >  CP009273/571597‑571894
                                                                                                                                                    |                                                                                                                                                     
gaagaaCTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGt                                                                                                                                                      >  2:4090508/1‑150 (MQ=255)
gaagaaCTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGt                                                                                                                                                      >  2:4141849/1‑150 (MQ=255)
gaagaaCTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGt                                                                                                                                                      >  1:2775245/1‑150 (MQ=255)
   gaaCTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGtttt                                                                                                                                                   <  1:1343684/150‑1 (MQ=255)
   gaaCTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGtttt                                                                                                                                                   <  1:3663091/150‑1 (MQ=255)
     aCTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTg                                                                                                                                                 >  2:3597346/1‑150 (MQ=255)
     aCTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTg                                                                                                                                                 >  1:2831967/1‑150 (MQ=255)
     aCTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTg                                                                                                                                                 >  1:739739/1‑150 (MQ=255)
      cTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGc                                                                                                                                                >  2:3257892/1‑150 (MQ=255)
       tttCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGcc                                                                                                                                               >  1:4427400/1‑150 (MQ=255)
       tttCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGcc                                                                                                                                               <  1:3194010/150‑1 (MQ=255)
        ttCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCt                                                                                                                                              >  2:460354/1‑150 (MQ=255)
         tCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCtt                                                                                                                                             <  1:2611850/150‑1 (MQ=255)
         tCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCtt                                                                                                                                             <  1:963209/150‑1 (MQ=255)
         tCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCtt                                                                                                                                             >  1:402120/1‑150 (MQ=255)
          cccTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTg                                                                                                                                            >  1:2283058/1‑150 (MQ=255)
           ccTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTgg                                                                                                                                           <  1:438713/150‑1 (MQ=255)
            cTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGt                                                                                                                                          <  1:2500041/150‑1 (MQ=255)
                aTTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTg                                                                                                                                      >  1:2430952/1‑150 (MQ=255)
                aTTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTg                                                                                                                                      >  1:2358961/1‑150 (MQ=255)
                    aCTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagca                                                                                                                                  <  1:1718103/150‑1 (MQ=255)
                    aCTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagca                                                                                                                                  <  1:1021619/150‑1 (MQ=255)
                    aCTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagca                                                                                                                                  >  1:4173250/1‑150 (MQ=255)
                     cTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagcag                                                                                                                                 <  1:3307978/150‑1 (MQ=255)
                      ttGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagcagc                                                                                                                                >  1:2960855/1‑150 (MQ=255)
                      ttGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagcagc                                                                                                                                >  1:3106653/1‑150 (MQ=255)
                      ttGAGAGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagcagc                                                                                                                                <  2:2001252/150‑1 (MQ=255)
                        gAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCaa                                                                                                                              <  2:1486110/150‑1 (MQ=255)
                         aGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCaaa                                                                                                                             >  1:2558154/1‑150 (MQ=255)
                          gtgtCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAg                                                                                                                            <  2:2178110/150‑1 (MQ=255)
                           tgtCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGt                                                                                                                           <  1:1447310/150‑1 (MQ=255)
                              cGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCa                                                                                                                        <  2:2886957/150‑1 (MQ=255)
                                gTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGa                                                                                                                      <  1:376215/150‑1 (MQ=255)
                                gTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGa                                                                                                                      >  2:736495/1‑150 (MQ=255)
                                  aCGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGAcc                                                                                                                    >  1:1773259/1‑150 (MQ=255)
                                   cGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTTCCTTGGTACTGAGCAGCAAAGTTCAGACCa                                                                                                                   >  1:539677/1‑150 (MQ=255)
                                       cAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCaa                                                                                                               >  1:4562724/1‑150 (MQ=255)
                                       cAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCaa                                                                                                               >  1:1900714/1‑150 (MQ=255)
                                        aGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAAc                                                                                                              >  2:2220500/1‑150 (MQ=255)
                                          aattCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCa                                                                                                            <  2:3476589/148‑1 (MQ=255)
                                          aTTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCa                                                                                                            >  1:4252535/1‑150 (MQ=255)
                                             tCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGAc                                                                                                         <  2:2230833/150‑1 (MQ=255)
                                              cGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGAcc                                                                                                        >  1:2552779/1‑150 (MQ=255)
                                                 aTAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCaaa                                                                                                     >  1:3565533/1‑150 (MQ=255)
                                                  tAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAaag                                                                                                    >  2:2100815/1‑150 (MQ=255)
                                                   aaGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAaaga                                                                                                   >  1:2903052/1‑150 (MQ=255)
                                                    aGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAaagaa                                                                                                  >  2:2281824/1‑150 (MQ=255)
                                                     gTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAaagaag                                                                                                 >  1:115664/1‑150 (MQ=255)
                                                     gTTGCACCGATACAGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTTTCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCTAAGTTCAGACATTCAACCAGACCAaagaag                                                                                                 <  2:1904076/150‑1 (MQ=255)
                                                      ttGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGt                                                                                                >  2:1182899/1‑150 (MQ=255)
                                                      ttGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGt                                                                                                <  2:4100334/150‑1 (MQ=255)
                                                       tGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTc                                                                                               >  1:1332051/1‑150 (MQ=255)
                                                        gCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCg                                                                                              <  2:1050370/150‑1 (MQ=255)
                                                          aCCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCgtt                                                                                            <  1:46831/150‑1 (MQ=255)
                                                          aCCGATACAGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCgtt                                                                                            <  1:2706281/150‑1 (MQ=255)
                                                              aTACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTa                                                                                        >  1:3272797/1‑150 (MQ=255)
                                                                  cGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGAt                                                                                    <  1:3286290/150‑1 (MQ=255)
                                                                    aaTCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAg                                                                                  >  1:3308365/1‑150 (MQ=255)
                                                                       ccTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGtt                                                                               <  2:3053997/150‑1 (MQ=255)
                                                                            gTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAac                                                                          >  1:1855153/1‑150 (MQ=255)
                                                                             tATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAaca                                                                         <  2:3367550/150‑1 (MQ=255)
                                                                             tATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAaca                                                                         <  2:1049784/150‑1 (MQ=255)
                                                                                tCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCa                                                                      <  1:757677/150‑1 (MQ=255)
                                                                                  aTAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTGGTTACGATAGGTTGCAACACCAGt                                                                    >  2:2964991/1‑150 (MQ=255)
                                                                                  aTAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGt                                                                    >  1:930370/1‑150 (MQ=255)
                                                                                   tAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGtt                                                                   >  2:1572384/1‑150 (MQ=255)
                                                                                      gTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCa                                                                >  2:1863573/1‑150 (MQ=255)
                                                                                        aGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACg                                                              >  1:4068910/1‑150 (MQ=255)
                                                                                         gCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGt                                                             >  1:3385218/1‑150 (MQ=255)
                                                                                          cAGAGAAACCGAAGCCATCACCGTTTCCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGtt                                                            <  2:2436321/150‑1 (MQ=255)
                                                                                          cAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGtt                                                            <  2:124876/150‑1 (MQ=255)
                                                                                                    gAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGaa                                                  >  1:3479475/1‑150 (MQ=255)
                                                                                                    gAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGaa                                                  >  2:2945154/1‑150 (MQ=255)
                                                                                                     aaGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAac                                                 >  1:946925/1‑150 (MQ=255)
                                                                                                            cACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCgg                                          >  2:1753955/1‑150 (MQ=255)
                                                                                                                gTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTg                                      >  2:96323/1‑150 (MQ=255)
                                                                                                                gTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTg                                      >  2:4460523/1‑150 (MQ=255)
                                                                                                                  tACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAg                                    <  2:1318721/150‑1 (MQ=255)
                                                                                                                  tACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAg                                    <  2:2443461/150‑1 (MQ=255)
                                                                                                                   aCCTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGt                                   >  1:3436024/1‑150 (MQ=255)
                                                                                                                    ccTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTc                                  >  2:1342815/1‑150 (MQ=255)
                                                                                                                     cTTCAGTGTAGTTATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTcc                                 <  1:4211702/150‑1 (MQ=255)
                                                                                                                     cTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCACAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTcc                                 <  2:353948/150‑1 (MQ=255)
                                                                                                                          gtgtAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGCTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAgtg                            <  2:642618/150‑1 (MQ=255)
                                                                                                                               gTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTcacc                       >  2:2448953/1‑150 (MQ=255)
                                                                                                                                 tATCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTcaccac                     >  1:4568107/1‑150 (MQ=255)
                                                                                                                                   tCGAAATCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCg                   <  2:2727380/150‑1 (MQ=255)
                                                                                                                                   tCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCg                   <  2:3920587/150‑1 (MQ=255)
                                                                                                                                      aaaTCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAAt                >  2:3077230/1‑150 (MQ=255)
                                                                                                                                      aaaTCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAAt                >  2:1387383/1‑150 (MQ=255)
                                                                                                                                       aaTCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAAtt               <  2:3723201/150‑1 (MQ=255)
                                                                                                                                       aaTCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAAtt               >  1:549594/1‑150 (MQ=255)
                                                                                                                                        aTCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTc              >  1:3009783/1‑150 (MQ=255)
                                                                                                                                         tCGCTACGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCt             <  2:356357/150‑1 (MQ=255)
                                                                                                                                            cTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGc          <  2:2276193/150‑1 (MQ=255)
                                                                                                                                             tACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGCa         >  1:3230196/1‑150 (MQ=255)
                                                                                                                                               cGATCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGCAgg       <  2:2679893/150‑1 (MQ=255)
                                                                                                                                                 aTCTTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGCAGGAc     >  2:1099034/1‑150 (MQ=255)
                                                                                                                                                   cGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGCAGGACGt   <  1:3948599/150‑1 (MQ=255)
                                                                                                                                                    gTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGCAGGACGTc  >  2:3063750/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACACCAGTTGCACGTTGAGTCATGAACACGTCGGTTTGAGTCCAAGTGTCACCACCGAATTCTGGCAGGACGTC  >  CP009273/571597‑571894

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: