Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,586,960 C→A 15.2% pseudogene (5285/5323 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,586,9600CA15.2% 1018.6 / 104.0 342pseudogene (5285/5323 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base C (145/144);  new base A (28/25);  total (173/169)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

TGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAG  >  CP009273/1586811‑1587108
                                                                                                                                                     |                                                                                                                                                    
   gATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCt                                                                                                                                                   <  2:2592759/150‑1 (MQ=255)
   gATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCt                                                                                                                                                   >  1:1640214/1‑150 (MQ=255)
   gATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCt                                                                                                                                                   >  1:2365273/1‑150 (MQ=255)
   gATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCt                                                                                                                                                   >  2:4148817/1‑150 (MQ=255)
       tCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCGGAGTATCGCTATAg                                                                                                                                               <  2:2005908/150‑1 (MQ=255)
         cgcgATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCTATAGcc                                                                                                                                             <  1:1437121/150‑1 (MQ=255)
          gcgATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCt                                                                                                                                            <  1:1112305/150‑1 (MQ=255)
            gATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCttt                                                                                                                                          >  2:79734/1‑150 (MQ=255)
             aTCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTa                                                                                                                                         <  1:3990565/150‑1 (MQ=255)
              tCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCTATAGCCTTTAt                                                                                                                                        <  2:892763/150‑1 (MQ=255)
               cgcgGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATc                                                                                                                                       <  2:4318966/150‑1 (MQ=255)
                 cgGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGTGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCAc                                                                                                                                     <  2:2395334/150‑1 (MQ=255)
                 cgGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCTATAGCCTTTATCAc                                                                                                                                     <  2:3004731/150‑1 (MQ=255)
                    aTAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACGAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTa                                                                                                                                  <  2:129770/150‑1 (MQ=255)
                    aTAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTa                                                                                                                                  >  2:1569805/1‑150 (MQ=255)
                    aTAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTa                                                                                                                                  >  1:5245348/1‑150 (MQ=255)
                     tAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCTATAGCCTTTATCACCTAg                                                                                                                                 >  2:30575/1‑150 (MQ=255)
                       aaaaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGtt                                                                                                                               <  2:2542346/150‑1 (MQ=255)
                         aaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGc                                                                                                                             <  1:5055455/150‑1 (MQ=255)
                         aaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGc                                                                                                                             >  1:3982206/1‑150 (MQ=255)
                         aaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGc                                                                                                                             >  1:3097367/1‑150 (MQ=255)
                          aggaggATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGca                                                                                                                            >  1:777291/1‑150 (MQ=255)
                           ggaggaTTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCCTCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGcac                                                                                                                           >  2:2300819/1‑150 (MQ=255)
                            gatgaTTTAAATAAGGCGTTGATGCGACGATTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGcaca                                                                                                                          <  1:1990855/147‑1 (MQ=255)
                                aTTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGa                                                                                                                      >  1:3884694/1‑150 (MQ=255)
                                 tttAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGAc                                                                                                                     <  2:1688922/150‑1 (MQ=255)
                                  ttAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCCGCATGCCCTGAGTATAGCTATAGCCTTTATCACCTAGTTGCACACCGACa                                                                                                                    >  1:3703695/1‑150 (MQ=255)
                                  ttAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCTATAGCCTTTATCACCTAGTTGCACACCGACa                                                                                                                    <  2:3852310/150‑1 (MQ=255)
                                   tAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACAt                                                                                                                   >  2:4965044/1‑150 (MQ=255)
                                     aaTAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTc                                                                                                                 <  1:3407168/150‑1 (MQ=255)
                                       tAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTccc                                                                                                               >  1:1977446/1‑150 (MQ=255)
                                             gTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAg                                                                                                         >  2:1114410/1‑150 (MQ=255)
                                             gTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTAATACGGTTTACCAGCTATATTTCACTCCCACCATGCCCTGAGTATCGCTATAGCCTTAATCACCTAGTAGCACACCGACATTCCCCCACAg                                                                                                         >  2:3694049/1‑150 (MQ=255)
                                              ttGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAgg                                                                                                        >  1:849271/1‑150 (MQ=255)
                                                   gcgcCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACTGGCTAAg                                                                                                   <  1:4807448/150‑1 (MQ=255)
                                                    cgcCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAgg                                                                                                  <  2:2680832/150‑1 (MQ=255)
                                                      ccGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGgtt                                                                                                >  2:3312767/1‑150 (MQ=255)
                                                      ccGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGgtt                                                                                                <  2:3178075/150‑1 (MQ=255)
                                                       cGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGgttg                                                                                               >  1:5179089/1‑150 (MQ=255)
                                                        gTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGgttgt                                                                                              >  1:1475967/1‑150 (MQ=255)
                                                         ttttAAGTCGGCAATATTAAGCCGCATGCCATCTCGACCTGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCGTTCTCACCTAGTTGAACACCGACATTCCCCCACAGGCTAAGgttgtt                                                                                             >  2:2958206/1‑150 (MQ=255)
                                                           tttaGTCGGCAAGATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTAt                                                                                           <  1:2722232/147‑1 (MQ=255)
                                                           ttaaGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTAt                                                                                           <  1:3015085/148‑1 (MQ=255)
                                                            taaGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTAtt                                                                                          >  1:2304872/2‑150 (MQ=255)
                                                              aGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTa                                                                                        >  1:1515322/1‑150 (MQ=255)
                                                                tCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACt                                                                                      <  1:1088284/150‑1 (MQ=255)
                                                                  ggCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACttt                                                                                    >  1:1716033/1‑150 (MQ=255)
                                                                  ggCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACttt                                                                                    >  1:2578867/1‑150 (MQ=255)
                                                                   gCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGTTTGTTATTTACTTTc                                                                                   >  1:65177/1‑150 (MQ=255)
                                                                    cAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCg                                                                                  >  1:2195446/1‑150 (MQ=255)
                                                                     aaTATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGc                                                                                 >  1:2339102/1‑150 (MQ=255)
                                                                       tatTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCt                                                                               >  1:1154620/1‑150 (MQ=255)
                                                                       tatTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCt                                                                               <  1:4491925/150‑1 (MQ=255)
                                                                        atTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTc                                                                              >  2:5202660/1‑150 (MQ=255)
                                                                        atTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTc                                                                              <  1:4730211/150‑1 (MQ=255)
                                                                           aaGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAAc                                                                           >  2:2835996/1‑150 (MQ=255)
                                                                                   tGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGgtac                                                                   <  2:3157738/150‑1 (MQ=255)
                                                                                   tGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGgtac                                                                   <  1:4327088/150‑1 (MQ=255)
                                                                                      cATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGgtacgta                                                                <  1:3217241/150‑1 (MQ=255)
                                                                                      cATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGgtacgta                                                                >  2:5032811/1‑150 (MQ=255)
                                                                                       aTCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGgtacgtac                                                               <  2:4581543/150‑1 (MQ=255)
                                                                                        tctcGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACt                                                              <  1:2829382/150‑1 (MQ=255)
                                                                                         ctcGACAGGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTTAGTATAGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGGACtt                                                             >  1:735249/1‑150 (MQ=255)
                                                                                           cGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCa                                                           >  2:5256623/1‑150 (MQ=255)
                                                                                            gACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCAc                                                          >  2:4664153/1‑150 (MQ=255)
                                                                                               aTGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCga                                                       >  1:1555551/1‑150 (MQ=255)
                                                                                               aTGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATAGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCga                                                       >  1:5081812/1‑150 (MQ=255)
                                                                                                   ggCTTATACGGTGTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACGTCACCGAGAtt                                                   <  1:4127117/150‑1 (MQ=255)
                                                                                                     cTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGCACTTCACCGAGATTAc                                                 >  2:1039233/1‑150 (MQ=255)
                                                                                                     cTTATACGGTATACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTAc                                                 <  1:4391460/150‑1 (MQ=255)
                                                                                                        ataCGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTAcgcg                                              >  1:559112/1‑150 (MQ=255)
                                                                                                             gTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCt                                         <  2:310968/150‑1 (MQ=255)
                                                                                                             gTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCt                                         >  1:5039077/1‑150 (MQ=255)
                                                                                                               ttACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTc                                       <  2:589846/150‑1 (MQ=255)
                                                                                                               ttACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTc                                       >  2:185624/1‑150 (MQ=255)
                                                                                                                tcccAGCTATATTTCACTCCCCGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACAGTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCa                                      <  2:2720822/148‑1 (MQ=255)
                                                                                                                 aCCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCAc                                     <  1:1776387/150‑1 (MQ=255)
                                                                                                                 aCCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCAc                                     <  1:2425398/150‑1 (MQ=255)
                                                                                                                  ccAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACg                                    <  2:3738256/150‑1 (MQ=255)
                                                                                                                   cAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACgg                                   >  1:3806093/1‑150 (MQ=255)
                                                                                                                   cAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGAAAACACGCAccaacacaa                                   >  2:3828998/1‑143 (MQ=255)
                                                                                                                     gCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGcc                                 >  2:4821654/1‑150 (MQ=255)
                                                                                                                      cTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCt                                <  2:1260014/150‑1 (MQ=255)
                                                                                                                      cTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCt                                >  1:2705101/1‑150 (MQ=255)
                                                                                                                             tCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGttt                         <  2:4924554/150‑1 (MQ=255)
                                                                                                                             tCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGttt                         <  2:67050/150‑1 (MQ=255)
                                                                                                                             tCACGCCGAGCTTGCCCTGAGTATCGCTAGAGCCTGTATCACCTAGGTGCACACCGACAGGCCCCCCCCGGCTCAGGTTGGTATTTACTGTCGCCTCACCCCCGGGACGGCCGTCACCGAGAGGACGCGCCCCTTCACGGCCTACGGttt                         <  2:4534720/150‑1 (MQ=255)
                                                                                                                                 tCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGAcc                     <  1:940307/150‑1 (MQ=255)
                                                                                                                                  cccAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCa                    >  2:3239403/1‑150 (MQ=255)
                                                                                                                                     aGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTc                 <  2:4385310/150‑1 (MQ=255)
                                                                                                                                     aGCATGCCCTGAGTATAGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTc                 >  2:2087487/1‑150 (MQ=255)
                                                                                                                                       cATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCAt               <  1:3576928/150‑1 (MQ=255)
                                                                                                                                          gCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCtt            <  1:4677894/150‑1 (MQ=255)
                                                                                                                                           cccTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTc           <  1:2978620/150‑1 (MQ=255)
                                                                                                                                           cccTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTc           <  2:1359967/150‑1 (MQ=255)
                                                                                                                                                aGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGc      >  1:4840091/1‑150 (MQ=255)
                                                                                                                                                aGTATAGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGc      >  2:3357128/1‑150 (MQ=255)
                                                                                                                                                 gTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCg     >  2:3930719/1‑150 (MQ=255)
                                                                                                                                                  tATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGt    >  1:4638927/1‑150 (MQ=255)
                                                                                                                                                    tCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAg  <  2:5293863/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAG  >  CP009273/1586811‑1587108

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: