Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,203 T→A 6.5% intergenic (+17/‑169) ydiK → / → rprA putative transporter YdiK/small regulatory RNA RprA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,770,2030TA 6.5% 269.7 / 7.9 92intergenic (+17/‑169)ydiK/rprAputative transporter YdiK/small regulatory RNA RprA
Reads supporting (aligned to +/- strand):  ref base T (41/45);  new base A (3/3);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-01

TGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAG  >  NC_000913/1770060‑1770349
                                                                                                                                               |                                                                                                                                                  
tGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccg                                                                                                                                              >  2:134024/1‑150 (MQ=255)
  gTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                                                                            >  1:609777/1‑150 (MQ=255)
  gTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                                                                            >  2:1256940/1‑150 (MQ=255)
  gTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                                                                            >  2:511585/1‑150 (MQ=255)
   tCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                                                                           >  1:365882/1‑150 (MQ=255)
    ccGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTg                                                                                                                                          <  1:1124644/150‑1 (MQ=255)
    ccGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTg                                                                                                                                          <  2:243737/150‑1 (MQ=255)
    ccGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCAGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTg                                                                                                                                          <  1:194301/150‑1 (MQ=255)
     cGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGGGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGa                                                                                                                                         <  2:381053/150‑1 (MQ=255)
      ggTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAt                                                                                                                                        <  1:275844/150‑1 (MQ=255)
      ggTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                                                                        <  1:448996/150‑1 (MQ=255)
       gTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAtt                                                                                                                                       <  1:1528218/150‑1 (MQ=255)
        ttCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                                                                      <  2:1247282/150‑1 (MQ=255)
          cTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGtt                                                                                                                                    <  1:287681/150‑1 (MQ=255)
           tGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTa                                                                                                                                   <  2:1039429/150‑1 (MQ=255)
           tGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTa                                                                                                                                   <  1:581565/150‑1 (MQ=255)
            gTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTaa                                                                                                                                  >  2:334794/1‑150 (MQ=255)
                gCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccac                                                                                                                              >  2:1415920/1‑150 (MQ=255)
                 ccGTTTCCTGGCGTCTTTTTGCCGCGTGGGTTGAAGAAGTCCCGCCGCCGACGGACCAACCTGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAATTAATTTCTGATCAGGCGGCTCTGCCGCCTGATTGTTAaccacc                                                                                                                             <  2:617031/150‑1 (MQ=255)
                 ccGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccacc                                                                                                                             <  1:1162889/150‑1 (MQ=255)
                   gTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCGCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                                                                           <  1:767024/150‑1 (MQ=255)
                    tttCCTGGCGTCTGGTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAAACACCGCt                                                                                                                          <  1:159403/150‑1 (MQ=255)
                      tCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTaa                                                                                                                        >  2:958618/1‑150 (MQ=255)
                        cTGGGGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                                                                      >  2:848623/1‑150 (MQ=255)
                         tGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTg                                                                                                                     <  1:199922/150‑1 (MQ=255)
                               cTGTTTGCCTCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATAttt                                                                                                               <  1:1196030/150‑1 (MQ=255)
                                     gCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttgc                                                                                                         >  1:620462/1‑150 (MQ=255)
                                      ccGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTtgcttgct                                                                                                        <  1:380310/150‑1 (MQ=255)
                                       cgcgTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctc                                                                                                       >  1:275929/1‑150 (MQ=255)
                                        gcgTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                                                                      <  1:180712/150‑1 (MQ=255)
                                           tgggtggAAGAAGTCCCGCCGCCGACGGCCCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTcc                                                                                                   >  1:1269172/1‑150 (MQ=255)
                                               tggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCat                                                                                               <  1:646424/150‑1 (MQ=255)
                                                ggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCata                                                                                              >  1:71090/1‑150 (MQ=255)
                                                   agaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAttt                                                                                           >  2:735243/1‑150 (MQ=255)
                                                    gaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGt                                                                                                                                     <  1:1390427/107‑1 (MQ=255)
                                                    gaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGt                                                                                                                                     >  2:1390427/1‑107 (MQ=255)
                                                    gaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtttt                                                                                          >  1:1428761/1‑150 (MQ=255)
                                                      agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTcc                                                                                        <  1:1317283/150‑1 (MQ=255)
                                                      agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGAGAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTcc                                                                                        >  1:301488/1‑150 (MQ=255)
                                                       gTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                                                                       >  1:919093/1‑150 (MQ=255)
                                                        tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                                                                      <  1:904140/150‑1 (MQ=255)
                                                         cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                                                                     >  1:1510634/1‑150 (MQ=255)
                                                           cgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                                                                       >  2:1456566/1‑146 (MQ=255)
                                                           cgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                                                                       <  1:1456566/146‑1 (MQ=255)
                                                             ccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTAc                                                                                 <  1:134024/150‑1 (MQ=255)
                                                                ccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGc                                                                              >  2:1406648/1‑150 (MQ=255)
                                                                 cgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGcc                                                                             >  1:1154682/1‑150 (MQ=255)
                                                                 cgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGcc                                                                             <  2:919093/150‑1 (MQ=255)
                                                                  gACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCa                                                                            >  2:597326/1‑150 (MQ=255)
                                                                   ccGGACCACCCGCAAGAACTGCTCGAAGCACAGGCCCAAATCGCGAAACCGCATAAGTAATATCTCATCAGGCGCCTCTGCCGCCTGATTGTTAACCACCGCTAATTCATATTTGCTTGCTGTTCCCCATATTTTCCTGCTTACCGCCAt                                                                           <  2:988756/149‑1 (MQ=255)
                                                                      gACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCa                                                                        <  1:1406648/150‑1 (MQ=255)
                                                                        ccAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTCAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGc                                                                      >  2:563254/1‑150 (MQ=255)
                                                                          aaCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTg                                                                    >  2:662295/1‑150 (MQ=255)
                                                                          aaCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTCATTGTTAACCACCGCTAATTGATATTTGATTGCTCTTCCCCATATTTTCCTGCTTACAGCCAATCAGCTg                                                                    >  1:1266822/1‑150 (MQ=255)
                                                                             cGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgta                                                                 >  2:283101/1‑150 (MQ=255)
                                                                               gaaganATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagt                                                               >  1:576137/1‑150 (MQ=255)
                                                                                  gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagtggc                                                            <  1:298971/150‑4 (MQ=255)
                                                                                  gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc                                                            <  2:620462/150‑1 (MQ=255)
                                                                                      ttCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTc                                                        >  1:79867/1‑150 (MQ=255)
                                                                                       tctcGAAGAACTTGGCGAAATCGATAAACCGAATAATTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACGGCTAATTGATATTTGCTTGGTCTTCCCCATTTTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTct                                                       <  2:543927/150‑2 (MQ=255)
                                                                                                aCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCtt                                              <  2:71090/150‑1 (MQ=255)
                                                                                                aCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCtt                                              >  1:133805/1‑150 (MQ=255)
                                                                                                 cTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCttt                                             <  2:1291563/150‑1 (MQ=255)
                                                                                                   tGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCCCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTAc                                           >  2:1243615/1‑150 (MQ=255)
                                                                                                       gAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGtt                                       >  2:764442/1‑150 (MQ=255)
                                                                                                           tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTa                                   <  2:275929/150‑1 (MQ=255)
                                                                                                            cGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTaa                                  <  1:371031/150‑1 (MQ=255)
                                                                                                            cGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTaa                                  >  2:327555/1‑150 (MQ=255)
                                                                                                               gaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACt                               <  2:834468/150‑1 (MQ=255)
                                                                                                               gaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACt                               <  2:609777/150‑1 (MQ=255)
                                                                                                                   ccGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATg                           <  2:365882/150‑1 (MQ=255)
                                                                                                                   ccGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATg                           <  2:779377/150‑1 (MQ=255)
                                                                                                                    cGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATga                          >  1:433302/1‑150 (MQ=255)
                                                                                                                     gAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATgag                         >  1:258459/1‑150 (MQ=255)
                                                                                                                     gAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGATAATTGATATTTGCTTTCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGTTAGTACCTGTCGCAAATTCTTTACATTTTTTAAACTAATgag                         <  2:846646/150‑1 (MQ=255)
                                                                                                                        tAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACg                      <  2:576137/150‑1 (MQ=255)
                                                                                                                          aGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGaa                    >  1:30930/1‑150 (MQ=255)
                                                                                                                            tAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCttt                                             >  1:1064065/1‑123 (MQ=255)
                                                                                                                            tAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCttt                                             <  2:1064065/123‑1 (MQ=255)
                                                                                                                             aaTTTCTCATCAGGCGGATCTGCCGCCTGATTGTTAACCACCGCTATTTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGGAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCt                 <  2:1266822/150‑1 (MQ=255)
                                                                                                                                ttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGAt              <  1:764442/150‑1 (MQ=255)
                                                                                                                                ttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGAt              >  2:741022/1‑150 (MQ=255)
                                                                                                                                  ctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCg            >  2:14079/1‑150 (MQ=255)
                                                                                                                                    catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGAc          >  2:393634/1‑150 (MQ=255)
                                                                                                                                    catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGAc          <  2:1154682/150‑1 (MQ=255)
                                                                                                                                     atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCTATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATga                          <  2:689734/133‑1 (MQ=255)
                                                                                                                                     atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCGATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATga                          >  1:689734/1‑133 (MQ=255)
                                                                                                                                     atcaGGCGGAAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTATTCCCCATATTTTACTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACg         >  2:224332/1‑150 (MQ=255)
                                                                                                                                       caGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCa       <  1:1169066/150‑1 (MQ=255)
                                                                                                                                        aGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCaa      >  1:1351470/1‑150 (MQ=255)
                                                                                                                                        aGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCaa      >  1:14929/1‑150 (MQ=255)
                                                                                                                                            ggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAg  <  1:1256940/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TGGTCCGGTTCTGTTAGCCGTTTCCTGGCGTCTGTTTGCCGCGTGGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAG  >  NC_000913/1770060‑1770349

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