Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,593,917 G→A 27.8% A648V (GCG→GTG)  tmcA ← elongator methionine tRNA (ac4C34) acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,593,9170GA27.8% 150.1 / 60.2 97A648V (GCG→GTG) tmcAelongator methionine tRNA (ac4C34) acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base A (13/14);  total (46/51)
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.84e-01

ACCATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTGTAAA  >  NC_000913/2593821‑2594015
                                                                                                |                                                                                                  
aCCATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCt                                                                                                 >  1:826763/1‑100 (MQ=255)
 ccATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTg                                                                                                >  2:230326/1‑100 (MQ=255)
 ccATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTg                                                                                                >  2:1873077/1‑100 (MQ=255)
 ccATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTg                                                                                                >  2:1784169/1‑100 (MQ=255)
 ccATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTg                                                                                                >  2:1329928/1‑100 (MQ=255)
  cATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTgc                                                                                               >  2:496372/1‑100 (MQ=255)
  cATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTgc                                                                                               >  1:1948851/1‑100 (MQ=255)
   aTGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTgcg                                                                                              >  1:1001070/1‑100 (MQ=255)
    tGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGGGCGGACGCTCAGTGCGAACATCATTAAGTGGGAACAACGCCTgcg                                                                                              >  1:1914698/1‑99 (MQ=255)
     gagaTAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGcgcg                                                                                            <  1:2135109/100‑1 (MQ=255)
      agaTAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGAGCgg                                                                                           >  2:925473/1‑100 (MQ=255)
        aTAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTgcgagccc                                                                                         >  2:1415351/1‑95 (MQ=255)
        aTAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGcc                                                                                         >  1:571554/1‑100 (MQ=255)
         tAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGcc                                                                                         <  1:392200/99‑1 (MQ=255)
          aaaaTTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCctc                                                                                       >  1:1509240/1‑100 (MQ=255)
           aaaTTGAAGGAGTCAGTGAAATAATTGCCTTTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACCTCATTAAGTGCGAACAGCGCCTGCGCGGCctct                                                                                      <  2:1441123/100‑1 (MQ=255)
            aaTTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCtt                                                                                     >  2:288194/1‑100 (MQ=255)
             aTTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTc                                                                                    >  2:1710548/1‑100 (MQ=255)
              ttGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTcc                                                                                   <  1:1644590/100‑1 (MQ=255)
                gAAGTAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCg                                                                                 >  1:49286/1‑100 (MQ=255)
                gAAGGAGTAAGTGAAATAATTGCCATTGTGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCg                                                                                 <  1:199040/100‑1 (MQ=255)
                 aaGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGc                                                                                >  1:1659879/1‑100 (MQ=255)
                 aaGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGc                                                                                <  1:1992734/100‑1 (MQ=255)
                 aaGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGc                                                                                <  1:2175166/100‑1 (MQ=255)
                  aGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCt                                                                               >  1:207414/1‑100 (MQ=255)
                   ggAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTg                                                                              <  1:783699/100‑1 (MQ=255)
                     aGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGAc                                                                            >  1:186311/1‑100 (MQ=255)
                        cAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTACGCGGCCTCTTCCCGCTGACGGa                                                                         >  1:355113/1‑100 (MQ=255)
                          gTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGAcc                                                                       <  1:1263748/100‑1 (MQ=255)
                          gTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGAcc                                                                       >  2:728572/1‑100 (MQ=255)
                          gTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGAcc                                                                       <  1:151670/100‑1 (MQ=255)
                          gTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGAcc                                                                       <  1:1166830/100‑1 (MQ=255)
                            gAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAg                                                                     >  1:2400391/1‑100 (MQ=255)
                            gAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAg                                                                     <  2:661205/100‑1 (MQ=255)
                            gAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAg                                                                     >  1:1943039/1‑100 (MQ=255)
                            gAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAg                                                                     >  1:1106910/1‑100 (MQ=255)
                            gAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAg                                                                     <  1:1507318/100‑1 (MQ=255)
                            gAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAg                                                                     >  1:1749575/1‑100 (MQ=255)
                            gAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAg                                                                     <  1:1695302/100‑1 (MQ=255)
                             aaaTAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGt                                                                    <  2:841370/100‑1 (MQ=255)
                             aaaTAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGt                                                                    <  2:2298888/100‑1 (MQ=255)
                              aataatTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTa                                                                   >  2:2082132/1‑100 (MQ=255)
                               ataatTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTaa                                                                  <  2:1806863/100‑1 (MQ=255)
                                  atTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACAt                                                               >  1:1663109/1‑100 (MQ=255)
                                   ttGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATc                                                              <  1:2136277/100‑1 (MQ=255)
                                       cATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGc                                                          <  2:468741/100‑1 (MQ=255)
                                         ttGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCgg                                                        <  2:1068240/100‑1 (MQ=255)
                                           gCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGcc                                                      <  2:480230/100‑1 (MQ=255)
                                            cGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCt                                                     <  2:2143018/100‑1 (MQ=255)
                                             gTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTg                                                    <  2:1749554/100‑1 (MQ=255)
                                             gTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTg                                                    >  2:1989628/1‑100 (MQ=255)
                                             gTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTg                                                    <  2:469975/100‑1 (MQ=255)
                                                aTGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCAGCCTGaaa                                                 >  2:518496/1‑100 (MQ=255)
                                                 tGCGATCGCTCAGACGCTCAGTGCACATAGCATTAAGTGCGAATAGCACTTGCGCTGCCTCTTTCCGGTGACTCACCAGTAACATCTTGCGGCCTGAAAg                                                <  2:1529042/100‑1 (MQ=255)
                                                 tGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGaaa                                                 >  2:380537/1‑99 (MQ=255)
                                                 tGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAg                                                <  1:2029806/100‑1 (MQ=255)
                                                  gCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGt                                               <  1:1813490/100‑1 (MQ=255)
                                                   cGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGtt                                              <  2:407053/100‑1 (MQ=255)
                                                    gATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGttt                                             >  1:810334/1‑100 (MQ=255)
                                                     aTCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTa                                            <  1:720462/100‑1 (MQ=255)
                                                     aTCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTa                                            >  1:2270234/1‑100 (MQ=255)
                                                      tCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTaa                                           >  2:501301/1‑100 (MQ=255)
                                                         cgcAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAgtg                                        >  1:64801/1‑100 (MQ=255)
                                                            aGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTgg                                     <  1:1630448/100‑1 (MQ=255)
                                                             gACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGt                                    <  2:403819/100‑1 (MQ=255)
                                                             gACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGt                                    <  2:403547/100‑1 (MQ=255)
                                                             gACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTATAAGTGTGGt                                    >  1:940186/1‑100 (MQ=255)
                                                              aCGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTa                                   <  2:1273426/100‑1 (MQ=255)
                                                              aCGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTa                                   <  1:817529/100‑1 (MQ=255)
                                                                gCTCAGTGCGAACATCTTTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACa                                 <  2:1400646/100‑1 (MQ=255)
                                                                   cAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTTCGGCCTGAAAGTTTAAGTGTGGTACATaa                              <  2:2180000/100‑1 (MQ=255)
                                                                    aGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAAc                             <  2:476937/100‑1 (MQ=255)
                                                                       gCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTgc                          <  1:1781422/100‑1 (MQ=255)
                                                                        cGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTgcg                         <  1:901416/100‑1 (MQ=255)
                                                                         gAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTgcgc                        >  2:1187862/1‑100 (MQ=255)
                                                                         gAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTgcgc                        >  1:1025912/1‑100 (MQ=255)
                                                                            catcatTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATc                     <  2:2028747/100‑1 (MQ=255)
                                                                            catcatTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATc                     <  1:518944/100‑1 (MQ=255)
                                                                             atcatTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCa                    <  2:635811/100‑1 (MQ=255)
                                                                                 ttAAGGGCGAACAGCTCTTTCGGGTCCTCATACCGCGGACTGACAAGTACCATGTTGCGTTGTTAAAGTTTATGTGTTGTACATAACTGCGCATCACTgg                <  2:1250203/100‑1 (MQ=255)
                                                                                   aaGTGCGAACAGCACCTGCGCGGCCTCTTCTCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGGTTAAGTGTGGTACATAACTGCGCATCACTGGCg              <  2:157123/100‑1 (MQ=255)
                                                                                    aGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCg              >  1:46333/1‑99 (MQ=255)
                                                                                    aGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGt             >  1:850292/1‑100 (MQ=255)
                                                                                     gTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGt             >  1:668402/1‑99 (MQ=255)
                                                                                     gTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGATTAAGTGTGGTACATAACTGCGCATCACTGGCGt             >  2:681913/1‑99 (MQ=255)
                                                                                      tGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGttt           <  1:2357523/100‑1 (MQ=255)
                                                                                       gCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGtttt          >  2:2401967/1‑100 (MQ=255)
                                                                                        cGAACAGCGCCTGCGCGTCCTCTTCCCGCTGACGCACAAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGtttt          >  2:1358025/1‑99 (MQ=255)
                                                                                        cGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGttttt         <  1:496375/100‑1 (MQ=255)
                                                                                        cGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGttttt         >  2:1358043/1‑100 (MQ=255)
                                                                                        cGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGttttt         >  2:1531585/1‑100 (MQ=255)
                                                                                         gAACAGCGCCTGCGCCGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTTGTACATAACTGCGCATCACTGGCGTTTTTc        <  2:858846/100‑1 (MQ=255)
                                                                                           aCAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTg      <  1:1948628/100‑1 (MQ=255)
                                                                                           aCAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATATCTGCGCATCACTGGCGTTTTTCTg      <  2:955769/100‑1 (MQ=255)
                                                                                            cAGCGCCTGCGCGTCTTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAATGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTGt     <  2:451970/100‑1 (MQ=255)
                                                                                             aGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTGTa    <  2:958909/100‑1 (MQ=255)
                                                                                             aGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGCGTGGTGCATAACTGAGCATCACTGGCGTTTTTCTGTa    >  2:1614974/1‑100 (MQ=255)
                                                                                              gcgcCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTGTaa   <  2:742380/100‑1 (MQ=255)
                                                                                               cgcCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTGTaaa  >  1:1804853/1‑100 (MQ=255)
                                                                                                |                                                                                                  
ACCATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTGTAAA  >  NC_000913/2593821‑2594015

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