Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,022,987 G→T A84A (GCG→GCT tatC → TatABCE protein translocation system subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,022,9870GT84.8% 238.6 / 37.7 92A84A (GCG→GCTtatCTatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base T (36/42);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

AGTTGCCGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTGGTTTCCAGCTCTCTGCT  >  NC_000913/4022890‑4023085
                                                                                                 |                                                                                                  
aGTTGCCGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTcc                                                                                                  <  2:1586829/100‑1 (MQ=255)
 gTTGCCGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCg                                                                                                 <  1:1148781/100‑1 (MQ=255)
  ttGCCGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCgg                                                                                                <  2:1067784/100‑1 (MQ=255)
   tGCCGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGt                                                                                               <  1:686459/100‑1 (MQ=255)
   tGCCGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGt                                                                                               <  2:1797982/100‑1 (MQ=255)
   tGCCGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGt                                                                                               <  2:445718/100‑1 (MQ=255)
   tGCCGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGt                                                                                               <  2:1672543/100‑1 (MQ=255)
      cGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGAt                                                                                            >  1:100445/1‑100 (MQ=255)
      cGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGAt                                                                                            <  2:594813/100‑1 (MQ=255)
       gCAAGGTTCAACGATGATCGCCAGCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTAAGCTCCGGTGAtt                                                                                           <  1:178849/100‑1 (MQ=255)
         aaGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATtct                                                                                         >  2:1424020/1‑100 (MQ=255)
          aGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATtctc                                                                                        <  2:300271/100‑1 (MQ=255)
           ggTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTGTGGTGTCGCTGATTCTGTCAGCTCCGGTGATTctct                                                                                       >  2:1852397/1‑100 (MQ=255)
            gTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTa                                                                                      >  1:1462531/1‑100 (MQ=255)
            gTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTa                                                                                      >  1:706624/1‑100 (MQ=255)
            gTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTa                                                                                      >  2:275099/1‑100 (MQ=255)
            gTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTa                                                                                      >  1:1375083/1‑100 (MQ=255)
              tCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATc                                                                                    >  1:698156/1‑100 (MQ=255)
              tCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATc                                                                                    >  1:634258/1‑100 (MQ=255)
               cAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCa                                                                                   <  2:396810/100‑1 (MQ=255)
                 aCGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAgg                                                                                 >  2:10289/1‑100 (MQ=255)
                     tgatCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGtgtg                                                                             <  1:260907/100‑1 (MQ=255)
                      tccgatctccGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTgg                                                                            <  1:1702487/92‑1 (MQ=255)
                      gatCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTgg                                                                            <  2:31808/100‑1 (MQ=255)
                      gatCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTgg                                                                            >  1:53387/1‑100 (MQ=255)
                        tCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGc                                                                          <  1:1856182/100‑1 (MQ=255)
                         cGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCa                                                                         <  2:1793246/100‑1 (MQ=255)
                         cGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCa                                                                         <  2:1730506/100‑1 (MQ=255)
                          gCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCAt                                                                        >  1:850869/1‑100 (MQ=255)
                          gCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCAt                                                                        >  1:850885/1‑100 (MQ=255)
                           ccACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCAtt                                                                       >  1:1539739/1‑100 (MQ=255)
                             aCCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTa                                                                     <  2:1080052/100‑1 (MQ=255)
                             aCCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCCTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTa                                                                     >  1:1752638/1‑100 (MQ=255)
                              ccGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTggagatcgga                                                                    >  2:1702498/1‑92 (MQ=255)
                               cGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTAt                                                                    <  1:932960/99‑1 (MQ=255)
                               cGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATc                                                                   <  1:1191397/100‑1 (MQ=255)
                                   gTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGccc                                                               <  1:1417641/100‑1 (MQ=255)
                                    tGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGcccc                                                              <  2:312519/100‑1 (MQ=255)
                                    tGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGcccc                                                              <  1:1727242/100‑1 (MQ=255)
                                       ccTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAgc                                                           >  2:649172/1‑100 (MQ=255)
                                        cTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAgcg                                                          >  1:274207/1‑100 (MQ=255)
                                          cGCCGTGCTTTACGCCGATCAAGCTGACCTTTATTGTGTCGCTGATTCGGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAgcgc                                                         <  2:354274/99‑1 (MQ=255)
                                             cGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGta                                                     >  1:1308432/1‑100 (MQ=255)
                                                 cTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAg                                                 <  2:952233/100‑1 (MQ=255)
                                                 cTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAg                                                 <  2:1347612/100‑1 (MQ=255)
                                                  tttACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGc                                                >  2:1582331/1‑100 (MQ=255)
                                                   ttACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCa                                               <  2:169179/100‑1 (MQ=255)
                                                     ncGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATg                                             >  1:89094/2‑100 (MQ=255)
                                                     aCGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATg                                             <  2:1901104/100‑1 (MQ=255)
                                                      cGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGa                                            >  1:792976/1‑100 (MQ=255)
                                                      cGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGa                                            >  1:1356395/1‑100 (MQ=255)
                                                      cGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGa                                            <  1:1814598/100‑1 (MQ=255)
                                                      cGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGa                                            <  1:1738552/100‑1 (MQ=255)
                                                       gCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGaa                                           >  1:565199/1‑100 (MQ=255)
                                                         cGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACg                                         >  1:403647/1‑100 (MQ=255)
                                                         cGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACg                                         <  2:1164009/100‑1 (MQ=255)
                                                         cGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACg                                         >  1:1004955/1‑100 (MQ=255)
                                                         cGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACg                                         >  2:1401224/1‑100 (MQ=255)
                                                           aTCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTc                                       >  2:221941/1‑100 (MQ=255)
                                                           aTCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTc                                       <  1:963843/100‑1 (MQ=255)
                                                             cAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGc                                     >  1:1596458/1‑100 (MQ=255)
                                                              aaGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGGATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGcc                                    <  1:1195387/100‑1 (MQ=255)
                                                               aGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCt                                   <  2:1422532/100‑1 (MQ=255)
                                                               aGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCt                                   >  1:1401353/1‑100 (MQ=255)
                                                               aGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCt                                   <  2:1462539/100‑1 (MQ=255)
                                                                 cTGACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCtgg                                 >  1:1886873/1‑100 (MQ=255)
                                                                  ngACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATACGCATGAACGTCGCCtggt                                >  2:734970/2‑100 (MQ=255)
                                                                   gACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCtggtg                               <  1:1148074/100‑1 (MQ=255)
                                                                   gACCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCtggtg                               >  2:1832413/1‑100 (MQ=255)
                                                                    aCCTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCtggtgg                              <  1:1794321/100‑1 (MQ=255)
                                                                     ccTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTCTCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTggtggt                             <  2:479/100‑1 (MQ=255)
                                                                     ccTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTggtggt                             <  2:1191409/100‑1 (MQ=255)
                                                                     ccTTTATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTggtggt                             >  1:330196/1‑100 (MQ=255)
                                                                       tttATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGc                           >  1:504959/1‑100 (MQ=255)
                                                                        ttATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGcc                          <  2:739331/100‑1 (MQ=255)
                                                                        ttATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGcc                          <  2:672285/100‑1 (MQ=255)
                                                                        ttATGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGcc                          <  1:1557494/100‑1 (MQ=255)
                                                                        ttATGGTGTCGCTGATGCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATGTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGcc                          <  2:1541026/100‑1 (MQ=255)
                                                                           tGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCgct                       <  2:1226222/100‑1 (MQ=255)
                                                                           tGGTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCgct                       >  2:1284815/1‑100 (MQ=255)
                                                                            ggTGTCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCgctg                      >  2:365300/1‑100 (MQ=255)
                                                                            ggTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCgctg                      >  1:675169/1‑100 (MQ=255)
                                                                             gtgtCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCgctgc                     <  1:1209017/100‑1 (MQ=255)
                                                                                tCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTgg                  <  1:1521596/100‑1 (MQ=255)
                                                                                tCGCTGATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTgg                  <  1:1424012/100‑1 (MQ=255)
                                                                                     gATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTGGTTTcc             <  2:330198/100‑1 (MQ=255)
                                                                                     gATTCTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTGGTTTcc             >  2:350454/1‑100 (MQ=255)
                                                                                         cTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTGGTTTCCAGct         >  2:1489614/1‑100 (MQ=255)
                                                                                         cTGTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTGGTTTCCAGc          >  1:1692358/1‑99 (MQ=255)
                                                                                           gTCAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGCCGCCTGGTGGTGCCGCTGCTGGTTTCCAGctc        >  2:1936389/1‑99 (MQ=255)
                                                                                            tCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTGGTTTCCAGctctc      >  1:917655/1‑100 (MQ=255)
                                                                                             cAGCTCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTGGTTTCCAGctctct     <  1:365298/100‑1 (MQ=255)
                                                                                                cgcCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTGGTTTCCAGCTCTctgct  >  1:19277/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
AGTTGCCGCAAGGTTCAACGATGATCGCCACCGACGTGGCCTCGCCGTTCTTTACGCCGATCAAGCTGACCTTTATGGTGTCGCTGATTCTGTCAGCGCCGGTGATTCTCTATCAGGTGTGGGCATTTATCGCCCCAGCGCTGTATAAGCATGAACGTCGCCTGGTGGTGCCGCTGCTGGTTTCCAGCTCTCTGCT  >  NC_000913/4022890‑4023085

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: