Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,632,141 T→C S166P (TCA→CCA)  b1545 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,632,1410TC100.0% 390.8 / NA 137S166P (TCA→CCA) b1545hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (66/71);  total (66/71)

Only 100 of 140 total aligned reads displayed.

ATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGC  >  NC_000913/1632051‑1632211
                                                                                          |                                                                      
aTAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCa                                                                      >  2:1336848/1‑93 (MQ=37)
 tAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAg                                                                     <  2:2660011/93‑1 (MQ=37)
 tAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAg                                                                     <  2:114365/93‑1 (MQ=37)
 tAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAg                                                                     <  1:5356189/93‑1 (MQ=255)
  aGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAgg                                                                    <  2:5957997/93‑1 (MQ=255)
   gTAAGAGCCCGCGGCGCGGGGAAACGTTTTGGCCGACCACCTGGGTTCAATGAAGACCAGAAAAAGGCGGAATTCGACCGAATTAAGCCAggg                                                                   >  2:193861/1‑92 (MQ=11)
   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                   <  1:1269535/93‑1 (MQ=255)
   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                   <  2:6189687/93‑1 (MQ=255)
   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                   <  2:5500451/93‑1 (MQ=37)
   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                   >  2:1096686/1‑93 (MQ=37)
   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                   >  1:338011/1‑93 (MQ=37)
   gTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgt                                                                   <  1:2325627/93‑1 (MQ=255)
    tAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtg                                                                  <  1:4089929/93‑1 (MQ=255)
    tAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtg                                                                  >  2:6055850/1‑93 (MQ=37)
     aaGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGgtgt                                                                 >  2:1549644/1‑93 (MQ=37)
      agagCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTa                                                                <  2:1266449/93‑1 (MQ=255)
       gagCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTaa                                                               >  1:1711884/1‑93 (MQ=37)
        agCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAg                                                              <  1:4005338/93‑1 (MQ=255)
           ccGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGtat                                                           >  2:2862/1‑93 (MQ=255)
              cggcgCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAg                                                        <  2:4723298/93‑1 (MQ=255)
                gcgcAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTg                                                      >  2:1085418/1‑93 (MQ=255)
                  gcATGGAAACGTTTTGGTCTACCACCTGTGTTAAATGAAGAACAGACACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGcc                                                    <  2:3983176/93‑1 (MQ=11)
                   cAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCa                                                   <  2:2428637/93‑1 (MQ=255)
                    aGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAt                                                  <  2:173040/93‑1 (MQ=255)
                     gggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                 >  1:4770185/1‑93 (MQ=255)
                     aggAAACGTGTTGGTCGACCACGTGTGTTAAAGGAAGAACAGAGACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtt                                                 <  2:2607515/92‑1 (MQ=11)
                      ggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCAtcg                                                <  2:6482602/93‑3 (MQ=255)
                      ggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTg                                                >  1:3183895/1‑93 (MQ=255)
                      ggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTg                                                <  1:2551557/93‑1 (MQ=255)
                      ggAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTg                                                <  2:5004283/93‑1 (MQ=255)
                          aCGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCg                                            <  2:5610146/93‑1 (MQ=255)
                          aCGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCg                                            <  2:1090601/93‑1 (MQ=255)
                            gTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCggg                                          >  1:5425438/1‑93 (MQ=255)
                            gTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCggg                                          >  2:836261/1‑93 (MQ=255)
                            gTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCggg                                          >  1:4520930/1‑93 (MQ=255)
                               ttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAt                                       >  1:6343127/1‑93 (MQ=255)
                               ttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAt                                       >  1:1774547/1‑93 (MQ=37)
                               ttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAt                                       <  1:4497138/93‑1 (MQ=255)
                               ttGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAt                                       <  1:1065443/93‑1 (MQ=255)
                                tGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAAtt                                      <  2:1780496/93‑1 (MQ=37)
                                 ggTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTc                                     <  2:4561667/93‑1 (MQ=255)
                                    cGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaa                                  <  1:2782275/93‑1 (MQ=255)
                                     gACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCaaaa                                 >  2:5641715/1‑93 (MQ=255)
                                      accaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                >  2:3644672/1‑93 (MQ=37)
                                      accaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAc                                >  2:1981285/1‑93 (MQ=255)
                                       ccaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAcc                               <  1:3419370/93‑1 (MQ=37)
                                       ccaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAcc                               >  2:2419646/1‑93 (MQ=255)
                                       ccaccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAAcc                               <  1:2737126/93‑1 (MQ=255)
                                        caccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCt                              <  1:6484903/93‑1 (MQ=25)
                                        caccTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCt                              <  2:5357638/93‑1 (MQ=25)
                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTc                             >  1:222618/1‑93 (MQ=25)
                                         accTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTc                             >  1:4054130/1‑93 (MQ=25)
                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgc                           >  1:3732426/1‑93 (MQ=18)
                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgc                           >  2:4673693/1‑93 (MQ=18)
                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgc                           <  1:3085092/93‑1 (MQ=18)
                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgc                           <  2:5630801/93‑1 (MQ=18)
                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgc                           <  2:5887692/93‑1 (MQ=18)
                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgc                           <  2:3809777/93‑1 (MQ=18)
                                           cTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTcgc                           >  2:2647917/1‑93 (MQ=18)
                                            tgtgTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGAAAGAATAATTGCCATTGCCCGGGAATTCAAAACCTggcg                          >  1:6171586/1‑93 (MQ=11)
                                              tgtTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGc                        >  1:5003848/1‑93 (MQ=21)
                                                ttAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCaa                      >  2:1762603/1‑93 (MQ=21)
                                                 tAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCaaa                     >  1:1921229/1‑93 (MQ=21)
                                                  aaaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                    >  1:6180907/1‑93 (MQ=18)
                                                  aaaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                    <  1:6272022/93‑1 (MQ=18)
                                                  aaaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAc                    >  1:5384290/1‑93 (MQ=21)
                                                   aaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAcc                   >  1:5452714/1‑93 (MQ=21)
                                                   aaTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAAcc                   <  2:6612942/93‑1 (MQ=18)
                                                    aTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCa                  >  1:1947298/1‑93 (MQ=21)
                                                    aTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCa                  >  2:3505788/1‑93 (MQ=18)
                                                    aTGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCa                  >  2:6500098/1‑93 (MQ=18)
                                                      gaagaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAtt                <  1:3093252/93‑1 (MQ=18)
                                                       aagaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAttt               <  1:2620805/93‑1 (MQ=18)
                                                        agaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCAtttt              >  2:2275508/1‑93 (MQ=21)
                                                         gaaCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTa             >  1:6311926/1‑93 (MQ=21)
                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag           >  1:1385464/1‑93 (MQ=21)
                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag           <  2:4049078/93‑1 (MQ=18)
                                                           aCAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAag           >  1:4457468/1‑93 (MQ=18)
                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga          <  2:5839244/93‑1 (MQ=18)
                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga          <  2:6396493/93‑1 (MQ=18)
                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga          <  1:6524908/93‑1 (MQ=18)
                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga          >  2:1380704/1‑93 (MQ=18)
                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga          >  2:1780913/1‑93 (MQ=18)
                                                            cAGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAaga          <  1:2604061/93‑1 (MQ=18)
                                                             aGAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAagag         <  1:5765651/93‑1 (MQ=18)
                                                              gAAACAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGc        >  1:5513491/1‑93 (MQ=18)
                                                               aaaCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGcc       <  2:44490/93‑1 (MQ=18)
                                                                aaCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCa      >  1:940816/1‑93 (MQ=18)
                                                                 aCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaa     >  1:933111/1‑93 (MQ=21)
                                                                 aCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaa     <  2:4198767/93‑1 (MQ=18)
                                                                 aCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaa     >  2:3754835/1‑93 (MQ=18)
                                                                 aCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaa     >  2:6379743/1‑93 (MQ=21)
                                                                 aCAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCCGCAAACCATTTTAAGAGCCaa     <  2:5509483/93‑1 (MQ=16)
                                                                  cAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaaa    <  2:2131319/93‑1 (MQ=18)
                                                                  cAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaaa    <  2:6702150/93‑1 (MQ=18)
                                                                  cAGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCaaa    >  2:4528054/1‑93 (MQ=18)
                                                                   aGGCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAg   <  1:1705493/93‑1 (MQ=18)
                                                                    ggCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGc  >  1:6694154/1‑93 (MQ=18)
                                                                    ggCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGc  <  2:301410/93‑1 (MQ=18)
                                                                    ggCGGTATTCGAACGAATTAAGCCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGc  <  2:5438357/93‑1 (MQ=18)
                                                                                          |                                                                      
ATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAATGAAGAACAGAAACAGGCGGTATTCGAACGAATTAAGTCAGGTGTAAGTATAAGTGCCATTGCCCGGGAATTCAAAACCTCGCGGCAAACCATTTTAAGAGCCAAAGC  >  NC_000913/1632051‑1632211

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: