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,9970GC99.3% 503.9 / ‑3.6 136P213P (CCC→CCGthiGthiamine biosynthesis ThiGH complex subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (74/61);  total (74/62)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 136 total aligned reads displayed.

GGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCA  >  NC_000913/4191901‑4192093
                                                                                                |                                                                                                
 gtggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGAt                                                                                              >  2:497227/1‑100 (MQ=255)
 gtggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGAt                                                                                              <  2:1880074/100‑1 (MQ=255)
 gtggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGAt                                                                                              <  1:216199/100‑1 (MQ=255)
 ctggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGAt                                                                                              <  1:1697220/99‑1 (MQ=255)
   ggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGAta                                                                                             >  2:1697226/1‑98 (MQ=255)
   ggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCg                                                                                            >  2:1001744/1‑100 (MQ=255)
   ggCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCg                                                                                            >  2:1788978/1‑100 (MQ=255)
    gCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGt                                                                                           >  2:188329/1‑100 (MQ=255)
       tGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcg                                                                                        <  2:1497645/100‑1 (MQ=255)
       tGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcg                                                                                        <  2:2873717/100‑1 (MQ=255)
         aGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcgcg                                                                                      <  1:2943207/100‑1 (MQ=255)
          gCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCcgcgc                                                                                     >  2:3249720/1‑99 (MQ=255)
           cAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGAc                                                                                    <  1:3371927/100‑1 (MQ=255)
           cAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGAc                                                                                    >  1:3027074/1‑100 (MQ=255)
           cAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGAc                                                                                    >  1:1342272/1‑100 (MQ=255)
             aaaTAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCACGACgg                                                                                  <  2:367711/100‑1 (MQ=255)
              aaTAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGc                                                                                 >  2:3098947/1‑100 (MQ=255)
                 aaCTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAAt                                                                              <  2:3123075/100‑1 (MQ=255)
                 aaCTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAAt                                                                              >  2:3436742/1‑100 (MQ=255)
                   cTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCg                                                                            >  2:3381695/1‑100 (MQ=255)
                   cTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCg                                                                            >  2:1429634/1‑100 (MQ=255)
                   cTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCg                                                                            >  2:2635186/1‑100 (MQ=255)
                    tGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGc                                                                           <  2:2879898/100‑1 (MQ=255)
                      cGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCTCCATGTTGACCGGATCGTCCGCGACGGCAATCGccn                                                                         <  1:1294600/100‑2 (MQ=255)
                      cGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGCCCGGATCGTCCGCGACGGCAATCGCCg                                                                         <  2:2247915/100‑1 (MQ=255)
                      cGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCg                                                                         <  2:2509758/100‑1 (MQ=255)
                       ggCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGt                                                                        >  2:389460/1‑100 (MQ=255)
                        gCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGCCCGGATCGTCCGCGACGGCAATCGCCGTa                                                                       >  2:2095147/1‑100 (MQ=255)
                        gCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTa                                                                       >  2:1785316/1‑100 (MQ=255)
                           gCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTc                                                                    <  2:3381332/100‑1 (MQ=255)
                           gCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTc                                                                    >  1:2238652/1‑100 (MQ=255)
                           gCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCCTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTc                                                                    >  2:695781/1‑100 (MQ=255)
                             ccGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCAc                                                                  <  2:156592/100‑1 (MQ=255)
                              cGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACt                                                                 <  2:1625253/100‑1 (MQ=255)
                              cGGTCCGGACTGACGTGCCAGTAGGACTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACt                                                                 <  1:291215/100‑1 (MQ=255)
                                  ccGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAaca                                                             >  1:3072518/1‑100 (MQ=255)
                                  ccGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAaca                                                             <  2:459546/100‑1 (MQ=255)
                                  ccGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAaca                                                             <  1:3098927/100‑1 (MQ=255)
                                   cGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAacac                                                            >  1:474525/1‑100 (MQ=255)
                                    ggACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACAcc                                                           >  1:1220388/1‑100 (MQ=255)
                                    ggACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACAcc                                                           <  1:1815686/100‑1 (MQ=255)
                                     gACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACcg                                                          >  1:2299725/1‑100 (MQ=255)
                                       cTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGATATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACcgcg                                                        >  1:1895304/1‑100 (MQ=255)
                                        tGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTAGTCACTAACACCGCGt                                                       >  2:3489900/1‑100 (MQ=255)
                                         gACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTc                                                      >  1:6300/1‑100 (MQ=255)
                                         gACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTc                                                      <  1:650025/100‑1 (MQ=255)
                                          aCGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCg                                                     <  2:162003/100‑1 (MQ=255)
                                           cGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCgg                                                    >  2:1126915/1‑100 (MQ=255)
                                           cGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCg                                                     >  1:557161/1‑99 (MQ=255)
                                             tGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGcc                                                  >  1:921703/1‑100 (MQ=255)
                                               ccAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGcccc                                                >  2:1464367/1‑100 (MQ=255)
                                                cAGTCGGCCTGCTTCCACCGCCCGCCGAAACGCCTTCGCCTTGTTGACCGGAGCGTCCGCTACGGCCATCTCCGTAGTCACTTACACCGCGTCGGcccct                                               <  1:2950025/100‑2 (MQ=255)
                                                cAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGccccc                                               >  1:2912108/1‑100 (MQ=255)
                                                cAGTAGGCCTGCTTCTACCGACAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGccccc                                               >  1:2162807/1‑100 (MQ=255)
                                                 aGTAGGCCTGCTTCTACCGCCAGACGAAATGCCATCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCa                                              >  1:3287082/1‑100 (MQ=255)
                                                    aGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCAttt                                           <  1:3204709/100‑1 (MQ=255)
                                                     ggCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTc                                          >  2:1535632/1‑100 (MQ=255)
                                                      gCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTcc                                         <  1:2247031/100‑1 (MQ=255)
                                                      gCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTCGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTcc                                         >  1:2172840/1‑100 (MQ=255)
                                                       ccTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCa                                        >  1:2004401/1‑100 (MQ=255)
                                                       ccTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCATCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCa                                        <  1:451395/100‑1 (MQ=255)
                                                            ntCTACCGCCAGACGAAATGCCTTCGCCATGTTGACCGTATCGTCCGCGACGGCAAACGCCGTATTCACTAACCCCGCGTCGGCCCCCATTTCCAGCGcc                                   >  1:3047493/2‑100 (MQ=255)
                                                               tACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTgc                                >  2:1744474/1‑100 (MQ=255)
                                                               tACCGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTgc                                >  2:1744595/1‑100 (MQ=255)
                                                                  cGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGcgcg                             >  1:2648072/1‑100 (MQ=255)
                                                                  cGCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGcgcg                             <  1:1811956/100‑1 (MQ=255)
                                                                  cGCCAGACGAAAAGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGcgcg                             >  2:2320997/1‑100 (MQ=255)
                                                                   gCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCgg                            >  2:3031200/1‑100 (MQ=255)
                                                                   gCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCgg                            >  1:2017403/1‑100 (MQ=255)
                                                                   gCCAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCgg                            >  1:709024/1‑100 (MQ=255)
                                                                     cAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCa                          >  2:1999592/1‑100 (MQ=255)
                                                                     cAGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCa                          >  2:236270/1‑100 (MQ=255)
                                                                      aGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCAt                         >  1:1836343/1‑100 (MQ=255)
                                                                      aGACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCAt                         >  2:2337473/1‑100 (MQ=255)
                                                                       gACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATg                        >  1:112107/1‑100 (MQ=255)
                                                                       gACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATg                        <  2:1342278/100‑1 (MQ=255)
                                                                       gACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATg                        >  1:3359837/1‑100 (MQ=255)
                                                                       gACGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATg                        <  2:1304357/100‑1 (MQ=255)
                                                                        aCGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACCAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATgg                       >  1:1348014/1‑100 (MQ=255)
                                                                         cGAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGc                      >  2:3498747/1‑100 (MQ=255)
                                                                          gAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGc                      <  1:2516989/99‑1 (MQ=255)
                                                                          gAAATGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCt                     >  1:2403182/1‑100 (MQ=255)
                                                                           aaaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTg                    >  2:1168205/1‑100 (MQ=255)
                                                                            aaTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTgg                   <  2:129273/100‑1 (MQ=255)
                                                                             aTGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTggg                  <  2:474527/100‑1 (MQ=255)
                                                                             aTGCCTTCGCCATGTTGACCGGATCGGCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTggg                  <  1:2196184/100‑1 (MQ=255)
                                                                              tGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTgggn                 <  1:885387/100‑2 (MQ=255)
                                                                              tGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGa                 <  1:2592092/100‑1 (MQ=255)
                                                                              tGCCTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGa                 <  1:1535626/100‑1 (MQ=255)
                                                                                 cTTCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAAcn              <  2:1801251/100‑2 (MQ=255)
                                                                                   tCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGcc            <  2:569215/100‑1 (MQ=255)
                                                                                   tCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGcc            <  1:345272/100‑1 (MQ=255)
                                                                                   tCGCCATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGcc            >  1:3422090/1‑100 (MQ=255)
                                                                                      ccATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGAt         >  1:404239/1‑100 (MQ=255)
                                                                                       cATGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATg        <  2:1020603/100‑1 (MQ=255)
                                                                                        aTGTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGc       >  1:387412/1‑100 (MQ=255)
                                                                                          gTTGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCa     <  1:497225/100‑1 (MQ=255)
                                                                                           ttGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAg    <  1:1402995/100‑1 (MQ=255)
                                                                                           ttGACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAg    <  2:207096/100‑1 (MQ=255)
                                                                                             gACCGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCa  >  2:399527/1‑100 (MQ=255)
                                                                                                |                                                                                                
GGTGGCATGAGCAAAATAACTGCGGCTGCCCGGTCCGGACTGACGTGCCAGTAGGCCTGCTTCTACCGCCAGACGAAATGCCTTCGCCATGTTGACGGGATCGTCCGCGACGGCAATCGCCGTATTCACTAACACCGCGTCGGCCCCCATTTCCAGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCA  >  NC_000913/4191901‑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: