Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,602,033 A→G V22A (GTC→GCC)  yjjB ← putative succinate exporter YjjB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,602,0330AG99.7% 1379.8 / NA 336V22A (GTC→GCC) yjjBputative succinate exporter YjjB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (144/191);  minor base C (0/1);  total (144/192)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

TGGTCATCAAGATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGACCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCAGTGCCATCACT  >  NC_000913/4601935‑4602130
                                                                                                  |                                                                                                 
  gTCATCAAGATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGc                                                                                               >  2:8949376/1‑101 (MQ=255)
    catcaAGATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAg                                                                                             >  2:8949381/1‑101 (MQ=255)
    catcaAGATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAg                                                                                             >  1:8949379/1‑101 (MQ=255)
      tcaAGATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGa                                                                                           >  1:8949384/1‑101 (MQ=255)
        aaGATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAAt                                                                                         <  1:8948785/101‑1 (MQ=255)
         aGATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATg                                                                                        >  1:8949389/1‑101 (MQ=255)
          gATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATgg                                                                                       >  2:8949391/1‑101 (MQ=255)
          gATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATgg                                                                                       <  2:8948756/101‑1 (MQ=255)
             cATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATggcgg                                                                                    <  2:8949039/101‑1 (MQ=255)
             cATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATggcgg                                                                                    <  2:8948627/101‑1 (MQ=255)
               ttCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGgcggcg                                                                                  >  1:8949397/1‑101 (MQ=255)
                 cGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAg                                                                                <  1:8948777/101‑1 (MQ=255)
                  gggAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAgg                                                                               >  2:8949401/1‑101 (MQ=255)
                  gggAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAgg                                                                               >  1:8949402/1‑101 (MQ=255)
                   ggAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATggcggc                                                                                   <  2:8948780/96‑1 (MQ=255)
                   ggAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGa                                                                              <  2:8948833/101‑1 (MQ=255)
                   ggAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGCCCGCAGGAATGGCGGCGAGGa                                                                              <  1:8948649/101‑1 (MQ=255)
                    gAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGAt                                                                             >  1:8949406/1‑101 (MQ=255)
                       ccATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCat                                                                          <  2:8948852/101‑1 (MQ=255)
                       ccATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCat                                                                          <  1:8948797/101‑1 (MQ=255)
                        cATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCata                                                                         >  2:8949412/1‑101 (MQ=255)
                         aTGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCatat                                                                        <  1:8948939/101‑1 (MQ=255)
                         aTGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCatat                                                                        >  1:8949414/1‑101 (MQ=255)
                         aTGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCatat                                                                        <  1:8948877/101‑1 (MQ=255)
                           gACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATcc                                                                      <  1:8948927/101‑1 (MQ=255)
                            aCCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCt                                                                     <  2:8949090/101‑1 (MQ=255)
                            aCCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCt                                                                     >  2:8949422/1‑101 (MQ=255)
                            aCCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCt                                                                     <  2:8948996/101‑1 (MQ=255)
                            aCCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCt                                                                     <  2:8948554/101‑1 (MQ=255)
                            aCCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCt                                                                     >  1:8949421/1‑101 (MQ=255)
                              cTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTgc                                                                   >  2:8949426/1‑101 (MQ=255)
                                aTCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATggcg                                                                                     <  2:8948685/81‑1 (MQ=255)
                                 tCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGcc                                                                >  1:8949429/1‑101 (MQ=255)
                                 tCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGcc                                                                <  1:8948704/101‑1 (MQ=255)
                                 tCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGcc                                                                <  2:8949057/101‑1 (MQ=255)
                                      ccAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGc                                                           >  1:8949436/1‑101 (MQ=255)
                                      ccAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGc                                                           <  2:8948460/101‑1 (MQ=255)
                                      ccAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGc                                                           >  2:8949435/1‑101 (MQ=255)
                                        aaGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTa                                                         >  1:8949440/1‑101 (MQ=255)
                                        aaGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTa                                                         <  2:8948844/101‑1 (MQ=255)
                                        aaGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTa                                                         <  2:8948909/101‑1 (MQ=255)
                                         agcagcGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGAt                                                                             >  1:8949443/1‑80 (MQ=255)
                                         agcagcGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCt                                                          <  1:8948689/99‑1 (MQ=255)
                                           cagcGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACa                                                      <  2:8948940/101‑1 (MQ=255)
                                           cagcGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACa                                                      >  2:8949445/1‑101 (MQ=255)
                                               gcACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGaaa                                                  >  1:8949449/1‑101 (MQ=255)
                                               gcACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGaaa                                                  <  1:8948986/101‑1 (MQ=255)
                                                cacaCCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAAt                                                 <  1:8949055/101‑1 (MQ=255)
                                                 acacCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAAtt                                                <  2:8948881/101‑1 (MQ=255)
                                                  cacCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTaa                                                        <  2:8948824/92‑1 (MQ=255)
                                                  cacCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTc                                               <  2:8948931/101‑1 (MQ=255)
                                                  cacCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTc                                               <  1:8949016/101‑1 (MQ=255)
                                                  cacCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTc                                               <  1:8949072/101‑1 (MQ=255)
                                                   acCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCg                                              >  2:8949454/1‑101 (MQ=255)
                                                     cAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGAt                                            <  1:8949181/101‑1 (MQ=255)
                                                     cAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGAt                                            <  1:8948825/101‑1 (MQ=255)
                                                      aGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATc                                           <  1:8948855/101‑1 (MQ=255)
                                                      aGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATc                                           <  2:8948946/101‑1 (MQ=255)
                                                      aGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATc                                           <  1:8949151/101‑1 (MQ=255)
                                                          tAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCacac                                       >  1:8949458/1‑101 (MQ=255)
                                                           aaCGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACAcc                                      <  1:8949028/101‑1 (MQ=255)
                                                            aCGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACAccc                                     >  2:8949461/1‑101 (MQ=255)
                                                            aCGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACAccc                                     >  1:8949463/1‑101 (MQ=255)
                                                            aCGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACAccc                                     <  2:8948884/101‑1 (MQ=255)
                                                             cGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCa                                    <  2:8948808/101‑1 (MQ=255)
                                                              gCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCAt                                   >  1:8949467/1‑101 (MQ=255)
                                                              gCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCAt                                   <  2:8948889/101‑1 (MQ=255)
                                                                ccGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATcc                                 >  1:8949471/1‑101 (MQ=255)
                                                                ccGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATcc                                 >  2:8949473/1‑101 (MQ=255)
                                                                 cGCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCcg                                <  1:8948733/101‑1 (MQ=255)
                                                                  gCACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCcgc                               >  1:8949477/1‑101 (MQ=255)
                                                                   cACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCcgcg                              >  2:8949480/1‑101 (MQ=255)
                                                                   cACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCcgcg                              <  2:8948744/101‑1 (MQ=255)
                                                                   cACGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCcgcg                              <  2:8949019/101‑1 (MQ=255)
                                                                    aCGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGt                             >  2:8949483/1‑101 (MQ=255)
                                                                    aCGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGt                             <  2:8948667/101‑1 (MQ=255)
                                                                     cGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTa                            <  2:8949113/101‑1 (MQ=255)
                                                                     cGGGAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTa                            >  2:8949486/1‑101 (MQ=255)
                                                                      gggAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAg                           <  1:8948552/101‑1 (MQ=255)
                                                                        gAACGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGcc                         >  1:8949491/1‑101 (MQ=255)
                                                                         aaCGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACCCCCATCCGCGt                             <  1:8948947/96‑1 (MQ=255)
                                                                          aCGTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGcccc                       >  2:8949493/1‑101 (MQ=255)
                                                                            gTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCa                     >  2:8949500/1‑101 (MQ=255)
                                                                            gTTGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCa                     >  1:8949498/1‑101 (MQ=255)
                                                                             ttGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCaa                    <  1:8948907/101‑1 (MQ=255)
                                                                              tGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCaaa                   <  2:8948863/101‑1 (MQ=255)
                                                                              tGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCaaa                   <  1:8949120/101‑1 (MQ=255)
                                                                              tGAACACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCaaa                   <  1:8948866/101‑1 (MQ=255)
                                                                                aaCACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATc                 <  2:8948706/101‑1 (MQ=255)
                                                                                aaCACCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATc                 >  1:8949503/1‑101 (MQ=255)
                                                                                   acCATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTc              >  1:8949508/1‑101 (MQ=255)
                                                                                    ccATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGAt                                            <  2:8948905/70‑1 (MQ=255)
                                                                                    ccATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCa             >  2:8949511/1‑101 (MQ=255)
                                                                                     cATCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCAg            <  2:8948972/101‑1 (MQ=255)
                                                                                       tCGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCAGTg          <  1:8948598/101‑1 (MQ=255)
                                                                                        cGCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCAGTGc         >  2:8949516/1‑101 (MQ=255)
                                                                                         gCAAAGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCAGTGcc        <  1:8949103/101‑1 (MQ=255)
                                                                                           aaaGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCAGTGCCAt      >  1:8949519/1‑101 (MQ=255)
                                                                                           aaaGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCAGTGCCAt      <  1:8949086/101‑1 (MQ=255)
                                                                                            aaGCCGGCCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCAGTGCCATc     >  1:8949523/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
TGGTCATCAAGATCATTCGGGAACCATGACCTATCGAGCCAAGCAGCGCACACCAGCGTAACGCCCGCACGGGAACGTTGAACACCATCGCAAAGCCGACCGCAGGAATGGCGGCGAGGATCATATCCTGCGCCAACGCTAACAGAAATTCGATCACACCCATCCGCGTAGCCCCCAAATCGTCAGTGCCATCACT  >  NC_000913/4601935‑4602130

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