Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,191,997 G→C P213P (CCC→CCG thiG ← thiamine biosynthesis ThiGH complex subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,191,9970GC98.9% 339.6 / ‑4.0 92P213P (CCC→CCGthiGthiamine biosynthesis ThiGH complex subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (43/48);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCA  >  NC_000913/4191900‑4192093
                                                                                                 |                                                                                                
tggtggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACcgg                                                                                                <  2:1113424/100‑1 (MQ=255)
 ggtggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGa                                                                                               >  2:2662575/1‑100 (MQ=255)
 ggtggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGa                                                                                               >  1:1562030/1‑100 (MQ=255)
  gtggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGAt                                                                                              <  1:2020440/100‑1 (MQ=255)
      cATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTc                                                                                          >  2:359687/1‑100 (MQ=255)
      cATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTc                                                                                          >  2:1068295/1‑100 (MQ=255)
       aTGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTcc                                                                                         <  2:1986274/100‑1 (MQ=255)
        tGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcg                                                                                        >  1:2506167/1‑100 (MQ=255)
         gAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcgc                                                                                       <  2:2311976/100‑1 (MQ=255)
         gAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcgc                                                                                       <  1:592637/100‑1 (MQ=255)
           gCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcgcg                                                                                      >  2:659750/1‑99 (MQ=255)
           gCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGa                                                                                     >  2:1759570/1‑100 (MQ=255)
            cAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGAc                                                                                    >  1:2211082/1‑100 (MQ=255)
             aaaaTAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACg                                                                                   <  1:792578/100‑1 (MQ=255)
              aaaTAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACgg                                                                                  <  1:2413140/100‑1 (MQ=255)
              aaaTAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACgg                                                                                  <  2:1111250/100‑1 (MQ=255)
                      gCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGcc                                                                          <  2:240937/100‑1 (MQ=255)
                         gCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTa                                                                       >  1:1505332/1‑100 (MQ=255)
                          cTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTAt                                                                      >  1:593197/1‑100 (MQ=255)
                           tGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTAtt                                                                     >  1:946813/1‑100 (MQ=255)
                             cccGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCa                                                                   >  1:727516/1‑100 (MQ=255)
                             cccGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCa                                                                   >  1:532469/1‑100 (MQ=255)
                              ccGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCAc                                                                  <  2:1831456/100‑1 (MQ=255)
                              ccGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTCGACCGGATCGTCCGCGACGGCAATCGCCGTATTCAc                                                                  >  1:1980123/1‑100 (MQ=255)
                               cGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACt                                                                 <  1:2646478/100‑1 (MQ=255)
                               cGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACt                                                                 <  2:550597/100‑1 (MQ=255)
                               cGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACt                                                                 <  1:138134/100‑1 (MQ=255)
                                     ggACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACAcc                                                           >  2:1218824/1‑100 (MQ=255)
                                     ggACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACAcc                                                           >  1:1755744/1‑100 (MQ=255)
                                        cTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACcgcg                                                        >  2:885652/1‑100 (MQ=255)
                                        cTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACcgcg                                                        <  2:88532/100‑1 (MQ=255)
                                         tGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCCTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGt                                                       <  2:1505326/100‑1 (MQ=255)
                                          gACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTc                                                      >  2:475050/1‑100 (MQ=255)
                                          gACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTc                                                      <  1:1771626/100‑1 (MQ=255)
                                             gTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGc                                                   <  2:1705541/100‑1 (MQ=255)
                                              tGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGcc                                                  >  1:2788989/1‑100 (MQ=255)
                                              tGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGcc                                                  >  1:1265642/1‑100 (MQ=255)
                                                 cAGTATGCCTGCTTCTACCGCCAGTCGAAATGCTTTCACCATGTTGACCGGCTCGTCCGCGACGGCAATCTTCGTATTCACTAACACCGCGTCGGcccc                                                <  2:568134/99‑1 (MQ=255)
                                                 cAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGccccc                                               >  2:2080686/1‑100 (MQ=255)
                                                     gggCCTGCTTCTCCCGCCAGACGAAATGCCTTTGCCATGTTGGCCGGATCGGCCGCGACGGCAATCTCCGTATTCACTAACACCGCGTCGGCCCCCAttt                                           <  2:1682432/99‑1 (MQ=255)
                                                     aGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCAttt                                           >  2:1074048/1‑100 (MQ=255)
                                                       gCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTcc                                         >  1:2332317/1‑100 (MQ=255)
                                                       gCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTcc                                         <  2:2250660/100‑1 (MQ=255)
                                                        ccTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCa                                        >  2:1325722/1‑100 (MQ=255)
                                                        ccTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCa                                        <  2:1322014/100‑1 (MQ=255)
                                                         cTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAg                                       <  2:1254079/100‑1 (MQ=255)
                                                             ttCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGcc                                   <  1:1218829/100‑1 (MQ=255)
                                                             ttCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGcc                                   >  1:2031686/1‑100 (MQ=255)
                                                               cTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTg                                 <  1:885655/100‑1 (MQ=255)
                                                                tACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTgc                                >  2:597403/1‑100 (MQ=255)
                                                                tACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTgc                                >  1:338838/1‑100 (MQ=255)
                                                                 aCCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTgcg                               <  1:1958592/100‑1 (MQ=255)
                                                                   cGCCAGACGAAATGCCTTCGCCATGTTGACCGGCTCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGcgcg                             >  1:1277921/1‑100 (MQ=255)
                                                                   cGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTgcgc                              >  1:1343889/1‑99 (MQ=255)
                                                                   cGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGcgcg                             >  1:2160323/1‑100 (MQ=255)
                                                                    gCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCgg                            >  2:503510/1‑100 (MQ=255)
                                                                    gCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCgg                            >  2:1532080/1‑100 (MQ=255)
                                                                      cAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCa                          >  1:2090713/1‑100 (MQ=255)
                                                                      cAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCa                          >  2:389593/1‑100 (MQ=255)
                                                                       aGACGACATGCCTTCGCCATGTTGACCGGATCTTCCGCGACGGCAATCACCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGcgcg                             <  2:1562169/96‑1 (MQ=255)
                                                                       aGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCAt                         >  2:55221/1‑100 (MQ=255)
                                                                       aGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCAt                         <  1:2412996/100‑1 (MQ=255)
                                                                       aGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCAt                         <  1:1728677/100‑1 (MQ=255)
                                                                        gACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATg                        >  1:1818616/1‑100 (MQ=255)
                                                                          cGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGc                      <  1:1588990/100‑1 (MQ=255)
                                                                            aaaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTg                    <  2:2211071/100‑1 (MQ=255)
                                                                            aaaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTg                    >  2:2755039/1‑100 (MQ=255)
                                                                             aaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTgg                   <  1:760178/100‑1 (MQ=255)
                                                                             aaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTgg                   <  1:756837/100‑1 (MQ=255)
                                                                              cgatcTTCGCCATGTTGACCGGATCGTCCGCGACGGCAAGCGCCGTATTCACTAACCCCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTggg                  <  1:2523220/96‑1 (MQ=255)
                                                                              aTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTggg                  <  1:2319645/100‑1 (MQ=255)
                                                                              aTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTgg                   <  1:2202289/99‑1 (MQ=255)
                                                                              aTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGACATGGCTggg                  <  2:2506155/100‑1 (MQ=255)
                                                                               tGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGa                 <  2:1552271/100‑1 (MQ=255)
                                                                               tGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGa                 <  2:204876/100‑1 (MQ=255)
                                                                                gCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGaa                >  1:365540/1‑100 (MQ=255)
                                                                                  cTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAAc               >  1:2352413/1‑99 (MQ=255)
                                                                                  cTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACg              >  1:5552/1‑100 (MQ=255)
                                                                                  cTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACg              <  1:1590127/100‑1 (MQ=255)
                                                                                  cTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACg              <  1:2484402/100‑1 (MQ=255)
                                                                                  cTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGGCGGCCCCCATTTCCAGCGCCTGCGCGTCATGGCTGGGAACg              <  2:1277916/100‑1 (MQ=255)
                                                                                   ttCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGc             <  1:1797239/100‑1 (MQ=255)
                                                                                    tCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGagatcg            >  2:2523208/1‑95 (MQ=255)
                                                                                    tCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGcc            <  2:310225/100‑1 (MQ=255)
                                                                                    tCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGcc            <  1:1633710/100‑1 (MQ=255)
                                                                                     cGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCg           <  1:2755051/100‑1 (MQ=255)
                                                                                         aTGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGc       <  1:1759579/100‑1 (MQ=255)
                                                                                          tGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGcc      <  2:1994121/100‑1 (MQ=255)
                                                                                          tGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGcc      <  1:12887/100‑1 (MQ=255)
                                                                                          tGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGcc      <  2:1343883/100‑1 (MQ=255)
                                                                                           gTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCa     >  1:271810/1‑100 (MQ=255)
                                                                                            ttGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAg    >  1:263592/1‑100 (MQ=255)
                                                                                              gACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCa  >  2:99923/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TGGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCA  >  NC_000913/4191900‑4192093

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