Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,259,550 A→T 6.6% intergenic (‑254/+169) psuK ← / ← fruA putative pseudouridine kinase/fructose‑specific PTS multiphosphoryl transfer protein FruA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,5500AT 6.6% 434.5 / 14.0 136intergenic (‑254/+169)psuK/fruAputative pseudouridine kinase/fructose‑specific PTS multiphosphoryl transfer protein FruA
Reads supporting (aligned to +/- strand):  ref base A (67/60);  new base T (5/4);  total (72/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 140 total aligned reads displayed.

TGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT  >  NC_000913/2259403‑2259696
                                                                                                                                                   |                                                                                                                                                  
tGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCATCCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGATGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaat                                                                                                                                                  <  2:1283984/150‑1 (MQ=255)
tGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaat                                                                                                                                                  >  1:1902150/1‑150 (MQ=255)
tGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaat                                                                                                                                                  <  2:1747214/150‑1 (MQ=255)
  cATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAaataa                                                                                                                                                >  1:1825596/1‑150 (MQ=255)
    tGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATc                                                                                                                                              <  1:102082/150‑1 (MQ=255)
    tGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATc                                                                                                                                              <  1:481785/150‑1 (MQ=255)
     gTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAGGGGCAGAAAATAATCt                                                                                                                                             >  2:1418058/1‑150 (MQ=255)
      ttttACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTg                                                                                                                                            <  1:1672957/150‑1 (MQ=255)
         tACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGccc                                                                                                                                         >  1:1422108/1‑150 (MQ=255)
           cGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGCCCtt                                                                                                                                       <  1:1431581/150‑1 (MQ=255)
            gtaaagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGCCCttt                                                                                                                                      <  2:1680468/150‑2 (MQ=255)
            gtaaagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTg                                                                                                                                      <  1:1115868/150‑1 (MQ=255)
               aagtaaaTTACCTCCAGAGAGGCAAAATGGGCAGATTTTAGCAGTCTCACAAAATTAGATTAAATTCCTTAATTACAAAAGCGAGGAATTTGTGACGGTATTCACAGAACAGAACAATAAAAAATGGCAGATTATTTTCTGACCtttttt                                                                                                                                   <  2:1531816/150‑5 (MQ=255)
                        aCCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCt                                                                                                                          >  1:728668/1‑150 (MQ=255)
                        aCCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCt                                                                                                                          <  1:701905/150‑1 (MQ=255)
                           gCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAg                                                                                                                       <  2:1825596/150‑1 (MQ=255)
                           gCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAg                                                                                                                       >  1:1102736/1‑150 (MQ=255)
                              agagagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAgtgt                                                                                                                    >  2:781375/1‑150 (MQ=255)
                               gagagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTACCCTTGAGGAAATAACCTGAGTGtt                                                                                                                   <  2:1577696/150‑1 (MQ=255)
                                 gagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtg                                                                                                                 >  2:1991684/1‑150 (MQ=255)
                                 gagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtg                                                                                                                 <  1:1663766/150‑1 (MQ=255)
                                   ggCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                                               <  1:836819/150‑1 (MQ=255)
                                    gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgatg                                                                                                              <  2:1827593/150‑1 (MQ=255)
                                             gCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATa                                                                                                                              >  2:335325/1‑125 (MQ=255)
                                             gCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGa                                                                                                     >  1:820807/1‑150 (MQ=255)
                                             gCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGa                                                                                                     >  1:1251100/1‑150 (MQ=255)
                                              cAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGaa                                                                                                    <  1:1370626/150‑1 (MQ=255)
                                                   ttGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTc                                                                                               <  1:1944501/150‑1 (MQ=255)
                                                   ttGAGCAGGCTCACAAAATGCGATTAAATGAATTAATTACAAAAGCGAGTAATTTGTGACGGGATGCACAGAACATAAAAATAAAAAAGGGCAGAAAATAATCTGCACTTGATGAAATAACCTGAGTGTTTGATTAACTTTTGAAGttta                                                                                               <  2:394572/150‑2 (MQ=255)
                                                         aGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATAAACGTTTGAAGTTTCTACTGa                                                                                         <  2:809399/150‑1 (MQ=255)
                                                            cTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATACCCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAg                                                                                      >  1:776512/1‑150 (MQ=255)
                                                             tCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGa                                                                                     >  1:592024/1‑150 (MQ=255)
                                                              cacaAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGAt                                                                                    >  1:1448987/1‑150 (MQ=255)
                                                                caAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATaa                                                                                  >  2:1166102/1‑150 (MQ=255)
                                                                  aaaTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGCTGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATaaaa                                                                                >  2:1871625/1‑150 (MQ=255)
                                                                        gATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAttt                                                                          >  2:494734/1‑150 (MQ=255)
                                                                         aTTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTa                                                                         >  1:622158/1‑150 (MQ=255)
                                                                           tAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTaaa                                                                       >  2:1419092/1‑150 (MQ=255)
                                                                            aaaTGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAAt                                                                      <  2:592024/150‑1 (MQ=255)
                                                                              aTGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTc                                                                    <  2:622158/150‑1 (MQ=255)
                                                                               tGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCa                                                                   <  1:629855/150‑1 (MQ=255)
                                                                               tGCCAGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCa                                                                   >  2:651720/1‑150 (MQ=255)
                                                                                gCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCAt                                                                  >  1:1610945/1‑150 (MQ=255)
                                                                                gCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCAt                                                                  <  2:1543556/150‑1 (MQ=255)
                                                                                 ccTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATa                                                                 <  2:492828/150‑1 (MQ=255)
                                                                                 ccTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATa                                                                 >  1:1822014/1‑150 (MQ=255)
                                                                                 ccTGAATTACAAAACCGACGAATTTGTGCCGGGATTCAAATAACATAACAATAAAAAACGGCACAAAATAATCTGCCCTTGATGAAAAAACCTTCGTGTTTGCTGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATa                                                                 <  2:915412/150‑1 (MQ=255)
                                                                                  cTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATaa                                                                <  2:1266209/150‑1 (MQ=255)
                                                                                    taaTTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATaaaa                                                              <  1:414816/149‑1 (MQ=255)
                                                                                    gAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATaaaa                                                              >  2:1037944/1‑150 (MQ=255)
                                                                                       ttACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTa                                                           <  2:1822014/150‑1 (MQ=255)
                                                                                       ttACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTa                                                           <  2:728668/150‑1 (MQ=255)
                                                                                         aCAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACTATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACa                                                         <  2:476082/150‑1 (MQ=255)
                                                                                         aCAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACa                                                         >  2:23102/1‑150 (MQ=255)
                                                                                         aCAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACa                                                         >  2:768364/1‑150 (MQ=255)
                                                                                          cAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACAt                                                        <  1:412792/150‑1 (MQ=255)
                                                                                          cAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACAt                                                        >  1:284983/1‑150 (MQ=255)
                                                                                             aaGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTca                                                     <  1:1037944/150‑1 (MQ=255)
                                                                                              aGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTcac                                                    <  2:453191/150‑1 (MQ=255)
                                                                                                   ggAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATa                                                                                   >  2:2292/1‑114 (MQ=255)
                                                                                                   ggAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATa                                                                                   <  1:2292/114‑1 (MQ=255)
                                                                                                     aaTTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCgg                                             <  2:1251100/150‑1 (MQ=255)
                                                                                                           tgACGGGATGCACATAACAGAACAATAAAAAATGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATtc                                       <  1:482522/150‑1 (MQ=255)
                                                                                                            gACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATtct                                      >  1:295644/1‑150 (MQ=255)
                                                                                                                 gATGCACAGAACAGAACAATAAAAAGGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTg                                 <  1:1418058/150‑1 (MQ=255)
                                                                                                                    tcacaGAACAGAACAATAAAAAAGGGCAGAAAATAAACTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTgcgc                              >  1:369227/2‑150 (MQ=255)
                                                                                                                     cacaGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                             >  2:481597/1‑150 (MQ=255)
                                                                                                                     cacaGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                             >  1:1454734/1‑150 (MQ=255)
                                                                                                                     cacaGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                             >  2:356000/1‑150 (MQ=255)
                                                                                                                       cagaacagaaCAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtca                           >  2:1673454/1‑150 (MQ=255)
                                                                                                                       cagaacagaaCAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtca                           >  2:271069/1‑150 (MQ=255)
                                                                                                                         gaacagaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtcatc                         <  2:735161/150‑1 (MQ=255)
                                                                                                                          aacagaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtcatca                        <  1:356000/150‑1 (MQ=255)
                                                                                                                          aacagaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtcatca                        <  2:1000430/150‑1 (MQ=255)
                                                                                                                          aacagaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtcatca                        <  2:1102736/150‑1 (MQ=255)
                                                                                                                            cagaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCaaa                      >  2:117480/1‑150 (MQ=255)
                                                                                                                              gaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAAtt                    >  1:1904265/1‑149 (MQ=255)
                                                                                                                              gaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATg                    >  2:63655/1‑150 (MQ=255)
                                                                                                                              gaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATg                    >  2:1435042/1‑150 (MQ=255)
                                                                                                                                 caATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTa                 >  1:797127/1‑150 (MQ=255)
                                                                                                                                 caATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTa                 >  2:1648049/1‑150 (MQ=255)
                                                                                                                                 aataaaaaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTa                 >  2:428272/5‑150 (MQ=255)
                                                                                                                                  aaTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTac                <  2:369227/150‑1 (MQ=255)
                                                                                                                                  aaTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTac                <  1:651720/150‑1 (MQ=255)
                                                                                                                                   ttAAAATAGGGGAGAAAATAATCTTCCCTTGAGGAACTAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTTCTGAAAGATAAACCCATTTAAATTCATAAAAATACATTCAGACACCGGGCATTCTCTTTGCGCTCATCAAATGTTaca               <  2:697674/149‑1 (MQ=255)
                                                                                                                                   aTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAACCCGGTCATTCTCTTTGCGCTCATCAAATGTTaca               >  2:1671053/1‑150 (MQ=255)
                                                                                                                                   aTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTaca               >  1:381765/1‑150 (MQ=255)
                                                                                                                                   aTAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTaca               >  2:1863139/1‑150 (MQ=255)
                                                                                                                                    tAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacag              >  2:293231/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagg             >  2:1503860/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagg             >  2:936937/1‑150 (MQ=255)
                                                                                                                                      aaaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga            >  1:1700737/1‑150 (MQ=255)
                                                                                                                                      aaaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagca            >  2:807360/1‑148 (MQ=255)
                                                                                                                                        caaGGGCAGATTATTTTCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggaca          >  2:1031388/2‑150 (MQ=255)
                                                                                                                                        aaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAAGCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggaca          <  2:820807/150‑1 (MQ=255)
                                                                                                                                        aaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggaca          >  1:867201/1‑150 (MQ=255)
                                                                                                                                         aaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggacag         <  1:1991684/150‑1 (MQ=255)
                                                                                                                                             gCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGaaa     >  1:489924/1‑150 (MQ=255)
                                                                                                                                             gCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGaaa     >  1:1371740/1‑150 (MQ=255)
                                                                                                                                                gAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAAttt  >  2:625704/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTT  >  NC_000913/2259403‑2259696

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: