Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,079,618 A→G 10.8% intergenic (‑288/+26) cmtB ← / ← tktA mannitol‑specific PTS enzyme IIA component CmtB/transketolase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,079,6180AG10.8% 320.7 / 25.1 120intergenic (‑288/+26)cmtB/tktAmannitol‑specific PTS enzyme IIA component CmtB/transketolase 1
Reads supporting (aligned to +/- strand):  ref base A (57/50);  new base G (3/10);  total (60/60)
Fisher's exact test for biased strand distribution p-value = 7.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

Only 100 of 120 total aligned reads displayed.

GATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAG  >  NC_000913/3079472‑3079763
                                                                                                                                                  |                                                                                                                                                 
gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCg                                                                                                                                                >  1:79397/1‑150 (MQ=255)
  tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGAc                                                                                                                                              >  1:273771/1‑150 (MQ=255)
  tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGAc                                                                                                                                              <  2:76315/150‑1 (MQ=255)
  tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAAAGGACATATCAAGGTAATAAAAAAGGTCGCAGAAGCGAc                                                                                                                                              >  2:1255896/1‑150 (MQ=255)
       caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACACCACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtttt                                                                                                                                         >  2:1441658/1‑150 (MQ=255)
       caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtttt                                                                                                                                         >  2:1588763/1‑150 (MQ=255)
       caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtttt                                                                                                                                         >  2:822343/1‑150 (MQ=255)
        aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCttttt                                                                                                                                        <  1:1099942/150‑1 (MQ=255)
        aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGAc                                                                                                                                              <  2:666056/144‑1 (MQ=255)
        aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGAc                                                                                                                                              >  1:666056/1‑144 (MQ=255)
         acaaATGTCAATACGCATATCGTGATTCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTCTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttttt                                                                                                                                       <  1:1365711/150‑2 (MQ=255)
         acaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttttt                                                                                                                                       <  1:962752/150‑2 (MQ=255)
         acaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTa                                                                                                                                       >  1:1347178/1‑150 (MQ=255)
          caaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttttta                                                                                                                                      <  2:409783/150‑3 (MQ=255)
          caaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttttta                                                                                                                                      <  2:883849/150‑3 (MQ=255)
          caaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTAc                                                                                                                                      >  1:1017018/1‑150 (MQ=255)
          caaATCTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttttta                                                                                                                                      <  1:688233/150‑3 (MQ=255)
           aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCttttttat                                                                                                                                     <  1:322446/150‑4 (MQ=255)
            aaTGTCAATACGCATATCCTGATGCATAATACATAATTGATACTTATAGCATATAAATACTTCTGATCTACAACACGCCTTATCTAATGCTTCTCCGCAAACGCACATATCAAGGTAATAAAACACGGCGCCGAAGCGACCTTTTTAccc                                                                                                                                    <  2:851260/150‑1 (MQ=255)
                 cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                               >  2:371260/1‑150 (MQ=255)
                 cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                               >  2:152027/1‑150 (MQ=255)
                   aTACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGc                                                                                                                             <  1:737225/150‑1 (MQ=255)
                    tACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCt                                                                                                                            <  1:794455/150‑1 (MQ=255)
                      cGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTaa                                                                                                                          <  2:277882/150‑1 (MQ=255)
                      cGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTCATCTACAACACGCGTTATCTATTGCTTGTCCGCAAACCGACATATCAAGGTAATAAAAAAGGTAGCCGAAGCGACCTTTTTACCCGAAATGCTaa                                                                                                                          <  2:884908/150‑1 (MQ=255)
                       tcATATCGTGATGCATATTACATAATTGATATATAATGAATATAAATATTTGTGATCTACAACACTCCTTATCTATTGCTTGTCCCCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAt                                                                                                                         <  1:1296074/149‑1 (MQ=255)
                       gCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAt                                                                                                                         <  2:169475/150‑1 (MQ=255)
                        cATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAtt                                                                                                                        <  1:303005/150‑1 (MQ=255)
                         atatCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTa                                                                                                                       <  2:3789/150‑1 (MQ=255)
                         atatCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTa                                                                                                                       <  1:1454130/150‑1 (MQ=255)
                          tatCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGc                                                                                                                                                 >  2:590063/1‑123 (MQ=255)
                                gATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCAAAATGCTAATTACAGCAGt                                                                                                                >  1:899006/1‑150 (MQ=255)
                                  tGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTc                                                                                                              >  2:514594/1‑150 (MQ=255)
                                       atTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTg                                                                                                         >  2:608199/1‑150 (MQ=255)
                                       atTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTg                                                                                                         <  1:1256292/150‑1 (MQ=255)
                                       atTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGAAAACGGACATATCAAGGTAATAAAAAAGGTCGACGAAGCGACATTTTTACACGAAATGCTAATTACAGCAGTTCTTTTg                                                                                                         >  1:1088443/1‑150 (MQ=255)
                                           cATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCttt                                                                                                     >  2:1049109/1‑150 (MQ=255)
                                             tAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCg                                                                                                   <  2:577212/150‑1 (MQ=255)
                                             tAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCg                                                                                                   <  2:62080/150‑1 (MQ=255)
                                               atttatttTTATTGCATTTAAATATTTGTGATCTACAACACGGCTTATCTATTGCTTGTCCTCAAACTGACATATCAAGGTATTAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTGTTTTGCTTTctga                                                                                                 <  2:1167593/150‑4 (MQ=255)
                                               aTTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCGGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGca                                                                                                 <  2:105393/150‑1 (MQ=255)
                                               aTTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGca                                                                                                 <  1:1032215/150‑1 (MQ=255)
                                               aTTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGca                                                                                                 <  2:1017018/150‑1 (MQ=255)
                                                    tatttattGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCaacaac                                                                                            >  2:243365/1‑150 (MQ=255)
                                                         attGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTAt                                                                                       >  1:1469680/1‑150 (MQ=255)
                                                             cATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAAc                                                                                   <  1:558441/150‑1 (MQ=255)
                                                             cATATAAATATTTGTGATATACAACACGCCTTATCTATTCCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAAAGCAGTTCTTTTGCTTTCGCAACAACGTTATCaaa                                                                                   >  2:575437/1‑149 (MQ=255)
                                                              atatAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACa                                                                                  <  2:933917/150‑1 (MQ=255)
                                                                ataAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGt                                                                                >  1:1430490/1‑150 (MQ=255)
                                                                  aaaTATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCGGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGGTCTTTTGCTTTCGCAACAACGGTATCACCAGTgc                                                                              >  2:1406366/1‑149 (MQ=255)
                                                                    atatTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAg                                                                            <  1:1703/150‑1 (MQ=255)
                                                                     tatTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGc                                                                           <  1:1242250/150‑1 (MQ=255)
                                                                       tttGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCg                                                                         >  2:181303/1‑150 (MQ=255)
                                                                        ttGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGa                                                                        <  1:1120154/150‑1 (MQ=255)
                                                                          gtgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                      <  1:459925/150‑1 (MQ=255)
                                                                          gtgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                      >  2:1523491/1‑150 (MQ=255)
                                                                           tgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAct                                                                     <  2:1250820/150‑1 (MQ=255)
                                                                            gATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActc                                                                    >  2:218655/1‑150 (MQ=255)
                                                                            gATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActc                                                                    >  2:1264387/1‑150 (MQ=255)
                                                                                  caacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACg                                                                                           <  1:1517327/121‑1 (MQ=255)
                                                                                     cacGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCa                                                                >  2:819902/1‑145 (MQ=255)
                                                                                     cacGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCa                                                                <  1:819902/145‑1 (MQ=255)
                                                                                         ccTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACagcagc                                                       >  2:905666/1‑150 (MQ=255)
                                                                                          cTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGct                                                      <  2:184769/150‑1 (MQ=255)
                                                                                                tATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCg                                                >  1:647403/1‑150 (MQ=255)
                                                                                                 aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCTAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCgg                                               <  2:899006/150‑1 (MQ=255)
                                                                                                 aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCgg                                               <  1:608199/150‑1 (MQ=255)
                                                                                                   tGCTTGTCCGCAAACTGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAAAAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAg                                             <  1:1012204/150‑1 (MQ=255)
                                                                                                   tGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATCCTAATTACAGCAGTTCTTTTGCCTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCggcg                                             <  2:1599451/150‑3 (MQ=255)
                                                                                                    gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                               >  1:501596/1‑115 (MQ=255)
                                                                                                    gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                               <  2:501596/115‑1 (MQ=255)
                                                                                                    gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGc                                            <  2:79397/150‑1 (MQ=255)
                                                                                                    gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGc                                            >  2:1516657/1‑150 (MQ=255)
                                                                                                    gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGc                                            >  2:1510975/1‑150 (MQ=255)
                                                                                                    gCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGc                                            >  2:499167/1‑150 (MQ=255)
                                                                                                     cTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                           <  2:1469680/150‑1 (MQ=255)
                                                                                                       tGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                               <  1:1001930/112‑1 (MQ=255)
                                                                                                       tGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                               >  2:1001930/1‑112 (MQ=255)
                                                                                                         tCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGAtt                                       >  1:1615299/1‑150 (MQ=255)
                                                                                                          ccGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTc                                      >  2:30667/1‑150 (MQ=255)
                                                                                                          ccGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTc                                      <  2:764102/150‑1 (MQ=255)
                                                                                                           cGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCa                                     >  1:1472102/1‑150 (MQ=255)
                                                                                                           cGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGACGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCa                                     >  1:371955/1‑150 (MQ=255)
                                                                                                              aaaCGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCg                                  >  1:89447/1‑150 (MQ=255)
                                                                                                                 cGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAg                               >  2:706890/1‑150 (MQ=255)
                                                                                                                  ggACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAgg                              <  2:2013/150‑1 (MQ=255)
                                                                                                                   gACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGT‑‑TCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtgg                           <  1:1508811/150‑1 (MQ=255)
                                                                                                                     cATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtgg                           >  1:544081/1‑150 (MQ=255)
                                                                                                                     cATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtgg                           >  2:46417/1‑150 (MQ=255)
                                                                                                                       tatCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTc                         >  1:361519/1‑150 (MQ=255)
                                                                                                                             ggTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACcg                   >  2:1155840/1‑150 (MQ=255)
                                                                                                                               taataaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACcgac                 >  1:171434/1‑150 (MQ=255)
                                                                                                                                  taaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGAt              <  1:1024914/150‑1 (MQ=255)
                                                                                                                                   aaaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATa             >  2:661408/1‑150 (MQ=255)
                                                                                                                                   aaaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATa             >  1:151592/1‑150 (MQ=255)
                                                                                                                                    aaaaaGGTCGCTTCGGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAg            >  1:1159115/1‑150 (MQ=255)
                                                                                                                                     aaaaGGTCGCTTCGGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGc           >  2:791981/1‑150 (MQ=255)
                                                                                                                                          gTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGt      >  1:1453966/1‑150 (MQ=255)
                                                                                                                                             gCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACCCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCa   >  2:264731/1‑150 (MQ=255)
                                                                                                                                              ccGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAg  <  1:243365/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAG  >  NC_000913/3079472‑3079763

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: