Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,013,366 (C)5→4 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,3620C.99.6% 1244.6 / ‑2.6 240coding (134/1659 nt)fliFflagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (109/130);  minor base G (1/0);  total (110/130)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

GAGTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCTTCAG  >  NC_000913/2013262‑2013458
                                                                                                    |                                                                                                
  gTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGccc                                                                                                >  2:2622217/1‑101 (MQ=255)
  gTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGccc                                                                                                <  2:2621854/101‑1 (MQ=255)
   tGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGcccc                                                                                               <  2:2621831/101‑1 (MQ=255)
   tGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGcccc                                                                                               <  2:2621906/101‑1 (MQ=255)
        tAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTa                                                                                         >  1:2622219/1‑101 (MQ=255)
         aaTCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTAc                                                                                        >  2:2622222/1‑101 (MQ=255)
          aTCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTAcc                                                                                       >  1:2622223/1‑101 (MQ=255)
            cGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTAc                                                                                        >  2:2622226/1‑98 (MQ=255)
            cGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGc                                                                                     <  1:2621664/101‑1 (MQ=255)
             gCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGca                                                                                    >  1:2622228/1‑101 (MQ=255)
             gCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGca                                                                                    <  2:2621803/101‑1 (MQ=255)
             gCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGca                                                                                    <  1:2621717/101‑1 (MQ=255)
              ccTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGcac                                                                                   <  1:2621946/101‑1 (MQ=255)
              ccTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGcac                                                                                   <  2:2621792/101‑1 (MQ=255)
               cTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGcaca                                                                                  <  2:2621877/101‑1 (MQ=255)
               cTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGcaca                                                                                  >  2:2622230/1‑101 (MQ=255)
                tgcgtgcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACat                                                                                 <  2:2622132/101‑1 (MQ=255)
                tgcgtgcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACat                                                                                 >  2:2622233/1‑101 (MQ=255)
                  cgtgcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACatta                                                                               <  1:2621785/101‑1 (MQ=255)
                    tgcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACattatt                                                                             <  2:2622046/101‑1 (MQ=255)
                     gcgAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTc                                                                            >  1:2622236/1‑101 (MQ=255)
                        aaTCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGc                                                                         <  1:2621865/101‑1 (MQ=255)
                        aaTCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGc                                                                         >  2:2622238/1‑101 (MQ=255)
                         aTCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCa                                                                        >  2:2622241/1‑101 (MQ=255)
                         aTCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCa                                                                        >  1:2622242/1‑101 (MQ=255)
                           ccGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAAt                                                                      >  1:2622246/1‑101 (MQ=255)
                              aaaaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTa                                                                                         <  2:2621576/79‑1 (MQ=255)
                              aaaaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCtt                                                                   >  1:2622249/1‑101 (MQ=255)
                               aaaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCttt                                                                  >  2:2622252/1‑101 (MQ=255)
                               aaaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCttt                                                                  <  1:2621631/101‑1 (MQ=255)
                               aaaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCttt                                                                  <  1:2621998/101‑1 (MQ=255)
                               aaaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCttt                                                                  <  1:2621812/101‑1 (MQ=255)
                                aaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTc                                                                 >  1:2622255/1‑101 (MQ=255)
                                aaTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTc                                                                 >  2:2622254/1‑101 (MQ=255)
                                 aTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTcc                                                                <  2:2621836/101‑1 (MQ=255)
                                 aTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTcc                                                                <  1:2621857/101‑1 (MQ=255)
                                 aTTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTcc                                                                >  1:2622257/1‑101 (MQ=255)
                                  ttCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCg                                                               <  1:2622045/101‑1 (MQ=255)
                                   tCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGa                                                              <  1:2622113/101‑1 (MQ=255)
                                    ccATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGAt                                                             <  2:2621816/101‑1 (MQ=255)
                                    ccATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGAt                                                             >  1:2622261/1‑101 (MQ=255)
                                      attgattgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCa                                                           >  2:2622266/1‑101 (MQ=255)
                                      attgattgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCa                                                           >  1:2622268/1‑101 (MQ=255)
                                        tgattgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAgg                                                         <  1:2621802/101‑1 (MQ=255)
                                        tgattgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAgg                                                         <  2:2621747/101‑1 (MQ=255)
                                        tgattgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAgg                                                         <  1:2621872/101‑1 (MQ=255)
                                         gattgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGa                                                        >  2:2622269/1‑101 (MQ=255)
                                          attgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAt                                                       <  2:2621965/101‑1 (MQ=255)
                                          attgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAt                                                       <  2:2622024/101‑1 (MQ=255)
                                          attgTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAt                                                       <  1:2621582/101‑1 (MQ=255)
                                           ttgttgCCTGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAtg                                                      <  2:2621783/101‑1 (MQ=255)
                                            tgttgCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAtgg                                                     <  2:2621993/101‑1 (MQ=255)
                                             gttgCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAtggt                                                    <  2:2622064/101‑1 (MQ=255)
                                             gttgCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAtggt                                                    >  2:2622274/1‑101 (MQ=255)
                                             gttgCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAtggt                                                    >  1:2622275/1‑101 (MQ=255)
                                              ttgCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGAtggtg                                                   <  1:2621891/101‑1 (MQ=255)
                                               tgCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGc                                                                                     >  1:2622277/1‑66 (MQ=255)
                                                   ggTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCa                                              >  2:2622281/1‑101 (MQ=255)
                                                     ttCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTcc                                                                >  1:2622284/1‑81 (MQ=255)
                                                     ttCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAAt                                            <  2:2621866/101‑1 (MQ=255)
                                                     ttCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAAt                                            <  2:2621942/101‑1 (MQ=255)
                                                      tCCGCGGCAGTGGCGGTCATGTTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAAtt                                           >  1:2622286/1‑101 (MQ=255)
                                                      tCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAAtt                                           <  1:2621741/101‑1 (MQ=255)
                                                      tCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAAtt                                           <  2:2622036/101‑1 (MQ=255)
                                                      tCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAAtt                                           <  2:2621954/101‑1 (MQ=255)
                                                       ccGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGc                                                                         <  1:2621765/70‑1 (MQ=255)
                                                       ccGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTg                                          >  2:2622287/1‑101 (MQ=255)
                                                       ccGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTg                                          <  2:2621870/101‑1 (MQ=255)
                                                        cgcgGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGt                                         >  1:2622292/1‑101 (MQ=255)
                                                        cgcgGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGt                                         >  2:2622290/1‑101 (MQ=255)
                                                         gcgGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTc                                        <  2:2621938/101‑1 (MQ=255)
                                                           ggCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAg                                      <  1:2622067/101‑1 (MQ=255)
                                                            gCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGc                                     >  2:2622295/1‑101 (MQ=255)
                                                             cAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGcc                                    >  1:2622299/1‑101 (MQ=255)
                                                             cAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGcc                                    >  2:2622298/1‑101 (MQ=255)
                                                              aGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGcc                                    <  2:2621764/100‑1 (MQ=255)
                                                                tGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAAc                                 <  1:2621878/101‑1 (MQ=255)
                                                                  gCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTg                               >  1:2622304/1‑101 (MQ=255)
                                                                   cGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGa                              >  2:2622307/1‑101 (MQ=255)
                                                                   cGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGa                              >  1:2622306/1‑101 (MQ=255)
                                                                    ggTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGAc                             >  2:2622309/1‑101 (MQ=255)
                                                                      tCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGc                           >  2:2622312/1‑101 (MQ=255)
                                                                       cATGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCa                          <  2:2621683/101‑1 (MQ=255)
                                                                        aTGGTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCaa                         >  1:2622315/1‑101 (MQ=255)
                                                                          ggTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAAt                       <  1:2621902/101‑1 (MQ=255)
                                                                          ggTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAAt                       >  1:2622318/1‑101 (MQ=255)
                                                                           gTCGCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATg                      <  1:2621977/101‑1 (MQ=255)
                                                                              gCACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAat                   <  1:2622145/101‑1 (MQ=255)
                                                                               cACTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAata                  <  2:2622002/101‑1 (MQ=255)
                                                                                aCTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAatat                 <  1:2621910/101‑1 (MQ=255)
                                                                                aCTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAatat                 >  2:2622320/1‑101 (MQ=255)
                                                                                 cTGATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGa                     <  2:2622160/96‑1 (MQ=255)
                                                                                   gATCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTcc              <  2:2621884/101‑1 (MQ=255)
                                                                                    aTCCTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCt             >  2:2622323/1‑101 (MQ=255)
                                                                                       cTGTGGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTAc          <  2:2621849/101‑1 (MQ=255)
                                                                                          tgGGCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGc       <  1:2621833/101‑1 (MQ=255)
                                                                                           gggCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCt      >  2:2622326/1‑101 (MQ=255)
                                                                                             gCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCTTc    >  1:2622327/1‑101 (MQ=255)
                                                                                             gCGAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCTTc    <  1:2621845/101‑1 (MQ=255)
                                                                                               gAAAGCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCTTCAg  <  2:2621921/101‑1 (MQ=255)
                                                                                                    |                                                                                                
GAGTGGCTTAATCGCCTGCGTGCGAATCCGAAAATTCCATTGATTGTTGCCGGTTCCGCGGCAGTGGCGGTCATGGTCGCACTGATCCTGTGGGCGAAAGCCCCCGACTACCGCACATTATTCAGCAATCTTTCCGATCAGGATGGTGGCGCAATTGTCAGCCAACTGACGCAAATGAATATTCCTTACCGCTTCAG  >  NC_000913/2013262‑2013458

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