Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,218,808 C→T G90D (GGT→GAT)  aer ← aerotaxis sensor receptor, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,218,8080CT100.0% 1364.7 / NA 315G90D (GGT→GAT) aeraerotaxis sensor receptor, flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (169/146);  total (169/146)

Only 100 of 324 total aligned reads displayed.

TCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTG  >  NC_000913/3218710‑3218905
                                                                                                  |                                                                                                 
tctcTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGtcatca                                                                                                 >  2:5261206/1‑101 (MQ=255)
  tctTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCAtt                                                                                               >  1:5261210/1‑101 (MQ=255)
     tCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCAttttt                                                                                            <  2:5260976/101‑1 (MQ=255)
      cATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTg                                                                                           >  2:5261215/1‑101 (MQ=255)
       aTCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCCTTGGCCCGCACCCAATAATGGTCATCATTTTTgc                                                                                          >  1:5261217/1‑101 (MQ=255)
         ccGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTgcgg                                                                                        <  2:5260909/101‑1 (MQ=255)
          cGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTgcggc                                                                                       <  2:5260508/101‑1 (MQ=255)
          cGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTgcggc                                                                                       <  1:5260818/101‑1 (MQ=255)
            tCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGa                                                                                     >  1:5261221/1‑101 (MQ=255)
              gCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGAtt                                                                                   >  2:5261223/1‑101 (MQ=255)
               cccGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGAttt                                                                                  >  1:5261228/1‑101 (MQ=255)
                  gggTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTc                                                                               >  1:5261232/1‑101 (MQ=255)
                  gggTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTc                                                                               <  2:5260628/101‑1 (MQ=255)
                   ggTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCa                                                                              <  1:5260670/101‑1 (MQ=255)
                    gTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCAc                                                                             <  1:5260688/101‑1 (MQ=255)
                    gTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCAc                                                                             <  1:5260680/101‑1 (MQ=255)
                     tACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACg                                                                            <  2:5261002/101‑1 (MQ=255)
                     tACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACg                                                                            >  1:5261236/1‑101 (MQ=255)
                     tACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTACGGCGATTTTTCACg                                                                            >  2:5261237/1‑101 (MQ=255)
                      aCGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGa                                                                           <  1:5260923/101‑1 (MQ=255)
                        gAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATg                                                                         <  1:5260831/101‑1 (MQ=255)
                        gAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATg                                                                         >  2:5261243/1‑101 (MQ=255)
                         aaTCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGc                                                                        >  1:5261246/1‑101 (MQ=255)
                           tCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCg                                                                      >  2:5261250/1‑101 (MQ=255)
                            cGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACg                                                                            >  2:5261253/1‑94 (MQ=255)
                             gACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCt                                                                    <  2:5260876/101‑1 (MQ=255)
                              aCATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTc                                                                   <  1:5260620/101‑1 (MQ=255)
                              aCATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTc                                                                   >  1:5261255/1‑101 (MQ=255)
                              aCATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTc                                                                   <  1:5260744/101‑1 (MQ=255)
                               cATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTcc                                                                  <  2:5260596/101‑1 (MQ=255)
                                atatAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCa                                                                 <  1:5260809/101‑1 (MQ=255)
                                atatAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCa                                                                 >  1:5261259/1‑101 (MQ=255)
                                atatAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCa                                                                 >  1:5261263/1‑101 (MQ=255)
                                 tataGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAg                                                                >  2:5261266/1‑101 (MQ=255)
                                  ataGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAgg                                                               <  1:5260659/101‑1 (MQ=255)
                                    aGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGc                                                             >  1:5261270/1‑101 (MQ=255)
                                    aGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGc                                                             >  2:5261271/1‑101 (MQ=255)
                                    aGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGc                                                             <  1:5261011/101‑1 (MQ=255)
                                     gCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCt                                                            >  1:5261275/1‑101 (MQ=255)
                                      ccACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTc                                                           <  2:5260771/101‑1 (MQ=255)
                                      ccACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTc                                                           >  2:5261280/1‑101 (MQ=255)
                                       cACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTcc                                                          >  1:5261281/1‑101 (MQ=255)
                                       cACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTcc                                                          <  1:5260718/101‑1 (MQ=255)
                                         cTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTcccc                                                        >  2:5261286/1‑101 (MQ=255)
                                          tGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCt                                                       <  1:5260889/101‑1 (MQ=255)
                                             tttttCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCtttt                                                    <  1:5260755/101‑1 (MQ=255)
                                              ttttCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCttttt                                                   >  2:5261290/1‑101 (MQ=255)
                                              ttttCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCttttt                                                   <  2:5260701/101‑1 (MQ=255)
                                                ttCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTc                                                 >  1:5261291/1‑101 (MQ=255)
                                                 tCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCa                                                >  2:5261294/1‑101 (MQ=255)
                                                   ccTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAgg                                              <  1:5260863/101‑1 (MQ=255)
                                                   ccTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAgg                                              <  2:5261072/101‑1 (MQ=255)
                                                    cTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAggg                                             <  2:5260806/101‑1 (MQ=255)
                                                    cTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAggg                                             <  2:5260664/101‑1 (MQ=255)
                                                     tCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGt                                            >  1:5261302/1‑101 (MQ=255)
                                                     tCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGt                                            <  1:5260975/101‑1 (MQ=255)
                                                     tCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGt                                            >  2:5261301/1‑101 (MQ=255)
                                                      cgcgCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTg                                           >  2:5261305/1‑101 (MQ=255)
                                                       gcgcACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCa                                                >  1:5261309/1‑95 (MQ=255)
                                                       gcgcACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGa                                          <  1:5260545/101‑1 (MQ=255)
                                                        cgcACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGaa                                         <  1:5260495/101‑1 (MQ=255)
                                                          cACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAAcc                                       >  2:5261311/1‑101 (MQ=255)
                                                           aCCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACca                                      >  1:5261317/1‑101 (MQ=255)
                                                           aCCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACca                                      <  2:5260862/101‑1 (MQ=255)
                                                           aCCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACca                                      <  1:5260591/101‑1 (MQ=255)
                                                           aCCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACca                                      <  2:5260740/101‑1 (MQ=255)
                                                           aCCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACca                                      >  2:5261314/1‑101 (MQ=255)
                                                              aTCGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACat                                   >  1:5261321/1‑101 (MQ=255)
                                                                cGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACcac                                     <  1:5260904/97‑1 (MQ=255)
                                                                cGGTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACatat                                 >  2:5261324/1‑101 (MQ=255)
                                                                  gTACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATcc                               <  2:5260817/101‑1 (MQ=255)
                                                                   tACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCg                              <  2:5260409/101‑1 (MQ=255)
                                                                   tACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCg                              >  2:5261328/1‑101 (MQ=255)
                                                                   tACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCg                              <  2:5260765/101‑1 (MQ=255)
                                                                   tACCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCg                              <  2:5260848/101‑1 (MQ=255)
                                                                    aCCGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGc                             <  1:5260794/101‑1 (MQ=255)
                                                                     ccGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAAc                                        <  2:5260728/89‑1 (MQ=255)
                                                                     ccGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCa                            >  1:5261332/1‑101 (MQ=255)
                                                                      cGCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCaa                           >  2:5261334/1‑101 (MQ=255)
                                                                       gCATTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCaaa                          <  2:5260837/101‑1 (MQ=255)
                                                                         aTTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAAcg                        >  1:5261340/1‑101 (MQ=255)
                                                                         aTTGGCCCGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAAcg                        >  2:5261338/1‑101 (MQ=255)
                                                                              cccGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCt                   >  1:5261347/1‑101 (MQ=255)
                                                                              cccGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCt                   >  2:5261346/1‑101 (MQ=255)
                                                                                cGCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCttt                 <  1:5260769/101‑1 (MQ=255)
                                                                                 gCACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCtttt                >  2:5261350/1‑101 (MQ=255)
                                                                                  cACCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTg               >  1:5261355/1‑101 (MQ=255)
                                                                                   aCCCAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTgg              >  2:5261357/1‑101 (MQ=255)
                                                                                     ccAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCCCTCCAGGGCTCCCCtttttt                                                  <  2:5260870/63‑1 (MQ=255)
                                                                                      cAATAATGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCat           >  2:5261361/1‑101 (MQ=255)
                                                                                       aataatGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCata          <  2:5260714/101‑1 (MQ=255)
                                                                                        ataatGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCatat         >  2:5261366/1‑101 (MQ=255)
                                                                                         taatGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATc        >  1:5261367/1‑101 (MQ=255)
                                                                                         taatGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATc        >  2:5261370/1‑101 (MQ=255)
                                                                                           atGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCg      >  2:5261374/1‑101 (MQ=255)
                                                                                            tGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCgg     >  1:5261380/1‑101 (MQ=255)
                                                                                            tGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCgg     <  1:5260636/101‑1 (MQ=255)
                                                                                            tGGTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCgg     >  2:5261378/1‑101 (MQ=255)
                                                                                              gTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGt   <  2:5260896/101‑1 (MQ=255)
                                                                                              gTCATCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGt   <  1:5261134/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCACTGATTTTTCCCTCGCGCACCATCGGTACCGCATTGGCCCGCACCCAATAATGGTCACCATTTTTGCGGCGATTTTTCACGATGCCGCTCCAGGGCTCCCCTTTTTTCAGGGTGAACCACATATCCGCAAACGCCGCTTTTGGCATATCCGGGTG  >  NC_000913/3218710‑3218905

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: