Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,013,366 (C)5→6 coding (138/1659 nt) fliF → flagellar basal‑body MS‑ring and collar protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,013,3611.C99.6% 1099.2 / NA 271A45P (GCC→CCC) fliFflagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (126/143);  minor base A (1/0);  total (128/143)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

GAGTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCTTC  >  NC_000913/2013262‑2013456
                                                                                                    |                                                                                               
gAGTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGc                                                                                                <  2:3370414/101‑1 (MQ=255)
gAGTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGc                                                                                                >  2:3371017/1‑101 (MQ=255)
 aGTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGcc                                                                                               <  1:3370308/101‑1 (MQ=255)
  gTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGccc                                                                                              >  1:3371018/1‑101 (MQ=255)
   tGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGcccc                                                                                             <  2:3370480/101‑1 (MQ=255)
    ggCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGccccc                                                                                            <  2:3370855/101‑1 (MQ=255)
     gCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGcccccc                                                                                           <  2:3370213/101‑2 (MQ=255)
     gCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGcccccc                                                                                           <  2:3370383/101‑2 (MQ=255)
     gCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGcccccc                                                                                           <  1:3370581/101‑2 (MQ=255)
     gCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGcccccc                                                                                           <  1:3370483/101‑2 (MQ=255)
      cTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGccccccg                                                                                          <  2:3370423/101‑3 (MQ=255)
      cTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGccccccg                                                                                          >  2:3371023/1‑99 (MQ=255)
      cTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGccccccg                                                                                          <  1:3370458/101‑3 (MQ=255)
       ttAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGccccccga                                                                                         <  2:3370074/101‑4 (MQ=255)
          aTCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTa                                                                                       >  1:3371028/1‑101 (MQ=255)
          aTCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTa                                                                                       <  2:3370353/101‑1 (MQ=255)
            cGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTAcc                                                                                     <  2:3370472/101‑1 (MQ=255)
             gCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCg                                                                                    >  1:3371032/1‑101 (MQ=255)
             gCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCg                                                                                    <  1:3370673/101‑1 (MQ=255)
              ccTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGc                                                                                   <  2:3370439/101‑1 (MQ=255)
               cTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGca                                                                                  <  1:3370333/101‑1 (MQ=255)
               cTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGca                                                                                  >  2:3371035/1‑101 (MQ=255)
               cTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGca                                                                                  <  1:3370657/101‑1 (MQ=255)
                tgcgtgcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGcac                                                                                 <  1:3370395/101‑1 (MQ=255)
                 gcgtgcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGcaca                                                                                <  2:3370662/101‑1 (MQ=255)
                   gtgcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACatt                                                                              >  1:3371040/1‑101 (MQ=255)
                    tgcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACatta                                                                             <  1:3370430/101‑1 (MQ=255)
                    tgcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACatta                                                                             <  1:3370469/101‑1 (MQ=255)
                     gcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCCCTGATCCTGTGGGCGAAAGCCCCCCGACTACCGca                                                                                  <  2:3370323/95‑1 (MQ=255)
                     gcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACattat                                                                            >  2:3371043/1‑101 (MQ=255)
                       gAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTc                                                                          <  1:3370515/101‑1 (MQ=255)
                       gAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTc                                                                          >  2:3371047/1‑101 (MQ=255)
                       gAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTc                                                                          <  1:3370524/101‑1 (MQ=255)
                        aaTCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGACCCCCGACTACCGCACATTATTCa                                                                         >  1:3371049/1‑101 (MQ=255)
                          tCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGTCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGc                                                                       <  2:3370542/101‑1 (MQ=255)
                           ccGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCa                                                                      >  2:3371052/1‑101 (MQ=255)
                            cGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCaa                                                                     <  2:3370714/101‑1 (MQ=255)
                             gAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAAt                                                                    >  1:3371055/1‑101 (MQ=255)
                              aaaaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATc                                                                   >  1:3371058/1‑101 (MQ=255)
                               aaaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCt                                                                  >  2:3371062/1‑101 (MQ=255)
                                aaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCtt                                                                 <  1:3370378/101‑1 (MQ=255)
                                 aTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCttt                                                                >  1:3371065/1‑101 (MQ=255)
                                   tCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTcc                                                              >  1:3371069/1‑101 (MQ=255)
                                    ccATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCCCATTATTCAGCAATCTTTCCg                                                             >  2:3371071/1‑101 (MQ=255)
                                    ccATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCg                                                             <  1:3370391/101‑1 (MQ=255)
                                      attgattgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGAt                                                           <  1:3370567/101‑1 (MQ=255)
                                        tgattgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCa                                                         >  2:3371075/1‑101 (MQ=255)
                                         gattgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAg                                                        >  1:3371078/1‑101 (MQ=255)
                                          attgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAgg                                                       >  2:3371081/1‑101 (MQ=255)
                                           ttgttgCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGa                                                      >  1:3371083/1‑101 (MQ=255)
                                             gttgCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAtg                                                    >  1:3371086/1‑101 (MQ=255)
                                                gCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAtggtg                                                 >  1:3371089/1‑101 (MQ=255)
                                                  cGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGgc                                               <  2:3370596/101‑1 (MQ=255)
                                                  cGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGgc                                               >  1:3371093/1‑101 (MQ=255)
                                                   ggTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGgcg                                              <  1:3370641/101‑1 (MQ=255)
                                                    gTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGgcgc                                             >  1:3371097/1‑101 (MQ=255)
                                                    gTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGgcgc                                             >  2:3371098/1‑101 (MQ=255)
                                                     ttCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCa                                            <  2:3370968/101‑1 (MQ=255)
                                                      tCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCaa                                           >  2:3371102/1‑101 (MQ=255)
                                                       ccGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAAt                                          <  1:3370500/101‑1 (MQ=255)
                                                         gcgGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTg                                        <  2:3370444/101‑1 (MQ=255)
                                                          cgGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAAtt                                         >  1:3371107/1‑99 (MQ=255)
                                                          cgGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGt                                       >  2:3371106/1‑101 (MQ=255)
                                                            gCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCa                                     >  2:3371112/1‑101 (MQ=255)
                                                            gCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCa                                     <  2:3370585/101‑1 (MQ=255)
                                                             cAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAg                                    >  1:3371114/1‑101 (MQ=255)
                                                              aGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAtg                                                    >  1:3371117/1‑84 (MQ=255)
                                                                 ggCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCaa                                >  1:3371121/1‑101 (MQ=255)
                                                                  gCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAAc                               <  1:3370610/101‑1 (MQ=255)
                                                                    ggTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTg                             >  1:3371124/1‑101 (MQ=255)
                                                                     gTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGa                            <  2:3370459/101‑1 (MQ=255)
                                                                       cATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACg                          <  2:3370429/101‑1 (MQ=255)
                                                                       cATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACg                          >  2:3371126/1‑101 (MQ=255)
                                                                        aTGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGcaca                                                                                <  1:3370449/46‑1 (MQ=39)
                                                                         tGGTCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCa                        >  1:3371128/1‑101 (MQ=255)
                                                                            tCGCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAAt                     <  1:3370408/101‑1 (MQ=255)
                                                                              gCACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGa                   >  2:3371132/1‑101 (MQ=255)
                                                                               cACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGaa                  <  1:3370819/101‑1 (MQ=255)
                                                                               cACTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGaa                  >  1:3371136/1‑101 (MQ=255)
                                                                                aCTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAat                 <  2:3370291/101‑1 (MQ=255)
                                                                                aCTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAat                 <  1:3370221/101‑1 (MQ=255)
                                                                                 cTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAata                <  1:3370690/101‑1 (MQ=255)
                                                                                 cTGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAata                >  2:3371139/1‑101 (MQ=255)
                                                                                  tGATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAatat               >  1:3371142/1‑101 (MQ=255)
                                                                                   gATCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATAtt              >  2:3371144/1‑101 (MQ=255)
                                                                                    aTCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTc             >  1:3371148/1‑101 (MQ=255)
                                                                                    aTCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTc             <  2:3370703/101‑1 (MQ=255)
                                                                                     tCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTcc            <  2:3370766/101‑1 (MQ=255)
                                                                                     tCCTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTcc            <  2:3370521/101‑1 (MQ=255)
                                                                                      cgtgtgGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCaa                       <  1:3370792/87‑1 (MQ=255)
                                                                                      ccTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCa                                     <  1:3370547/75‑1 (MQ=255)
                                                                                      ccTGTGGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCt           <  2:3370616/101‑1 (MQ=255)
                                                                                        tgtgGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTa         >  2:3371151/1‑101 (MQ=255)
                                                                                         gtgGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTAc        <  2:3370636/101‑1 (MQ=255)
                                                                                          tgGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTAcc       <  1:3370856/101‑1 (MQ=255)
                                                                                          tgGGCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTAcc       >  2:3371155/1‑101 (MQ=255)
                                                                                           gggCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCg      >  1:3371156/1‑101 (MQ=255)
                                                                                             gCGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCt    <  2:3370563/101‑1 (MQ=255)
                                                                                              cGAAAGCCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCtt   <  1:3370628/101‑1 (MQ=255)
                                                                                                 acgCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGgc                                               <  2:3370508/51‑1 (MQ=255)
                                                                                                    |                                                                                               
GAGTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCTTC  >  NC_000913/2013262‑2013456

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