Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,712,402 A→G 42.6% S121S (AGT→AGC trmN ← tRNA1(Val) (adenine(37)‑N6)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,712,4020AG42.6% 34.9 / 109.7 92S121S (AGT→AGCtrmNtRNA1(Val) (adenine(37)‑N6)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (30/22);  new base G (18/21);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

AGGTCAGCAACGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGACTGGTTGAT  >  NC_000913/2712305‑2712500
                                                                                                 |                                                                                                  
aGGTCAGCAACGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCt                                                                                                  <  2:264586/100‑1 (MQ=255)
 ggTCAGCAACGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCtg                                                                                                 <  1:1328415/100‑1 (MQ=255)
 ggTCAGCAACGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACtg                                                                                                 <  1:1979869/100‑1 (MQ=255)
  gTCAGCAACGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACtga                                                                                                <  1:1085265/100‑1 (MQ=255)
    cAGCAACGATGGGTGATCAAGCGTAGTGGTATAGTGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTCGGTTGCtgatg                                                                                              <  2:1261151/100‑1 (MQ=255)
    cAGCAACGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACtgatg                                                                                              >  2:1485479/1‑100 (MQ=255)
     aGCAACGATGGGTGATCAAGCGTAGTGGTGTAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTTCTCCCTGCTGATAGTAAGGTGGGTTGCtgatga                                                                                             <  2:1433861/100‑1 (MQ=255)
     aGCAACGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACtgatga                                                                                             <  2:2341742/100‑1 (MQ=255)
      gCAACGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTgatgat                                                                                            >  2:2366227/1‑100 (MQ=255)
        aaCGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTa                                                                                          >  2:140599/1‑100 (MQ=255)
            aTGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATc                                                                                      <  1:1730910/100‑1 (MQ=255)
             tGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCg                                                                                     <  2:1009479/100‑1 (MQ=255)
                gTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAg                                                                                  <  2:506089/100‑1 (MQ=255)
                  gATAAAGCGTAGTAATATAGCGAGCCTTTTCCCGTTGAGGTGTCGAGCACGCTACTCCTTGCTGATAGTAAGGGGTGTTGCGGATGAGGAAATCGAAGCg                                                                                >  2:2003926/1‑100 (MQ=255)
                   aTCAAGCGTAGTGGTATAGCGAGCATGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGt                                                                               >  2:384544/1‑100 (MQ=255)
                    tCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTa                                                                              >  2:419750/1‑100 (MQ=255)
                     cAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTAc                                                                             >  1:1163652/1‑100 (MQ=255)
                     cAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTAc                                                                             <  1:1858622/100‑1 (MQ=255)
                     cAAGCGTAGGGGTATAGCGAGCCTGTGCCCGTTGAGGTGTTGAGCAATCTACTCCCTGCGGATAGTAAGGTGTGTTGTTGATGATTAAATCGAAGCGTAc                                                                             <  2:1389234/100‑1 (MQ=255)
                      aaGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACt                                                                            >  2:1071038/1‑100 (MQ=255)
                       aGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTg                                                                           >  2:1200552/1‑100 (MQ=255)
                       aGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTg                                                                           >  2:1631192/1‑100 (MQ=255)
                       aGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTg                                                                           >  1:235863/1‑100 (MQ=255)
                       aGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTg                                                                           >  1:236602/1‑100 (MQ=255)
                       aGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTg                                                                           >  2:2379824/1‑100 (MQ=255)
                         cGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTc                                                                         <  2:1429618/100‑1 (MQ=255)
                         cGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTc                                                                         <  1:1705044/100‑1 (MQ=255)
                            aGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAATGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTctgc                                                                      <  2:1472534/100‑1 (MQ=255)
                            aGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTctgc                                                                      <  1:1200563/100‑1 (MQ=255)
                            aGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTctgc                                                                      <  1:1377094/100‑1 (MQ=255)
                               ggTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCtgt                                                                   >  1:161118/1‑100 (MQ=255)
                                 tataGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCtgtgt                                                                 <  2:1814998/100‑1 (MQ=255)
                                    aGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCtgtgt                                                                 <  2:170520/97‑1 (MQ=255)
                                    aGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGAt                                                              >  1:1042300/1‑100 (MQ=255)
                                    aGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCCGTGTGAt                                                              >  2:268749/1‑100 (MQ=255)
                                    aGCGAGCCTGTTCCCGTTGAGGTGGTGAGCACTCTACTCCCTTCTGATAGGAAGGTGGGTTCCTGATGATTAAATCGAAGCGTACTATCTGCTGTATGAt                                                              >  2:1165379/1‑100 (MQ=255)
                                     gCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATc                                                             <  2:37842/100‑1 (MQ=255)
                                      cGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACAGTCTGCTATGTGAt                                                              >  1:81838/1‑98 (MQ=255)
                                      cGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATcc                                                            <  1:1662440/100‑1 (MQ=255)
                                      cGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATcc                                                            >  1:1566802/1‑100 (MQ=255)
                                      cGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATcc                                                            <  1:1492484/100‑1 (MQ=255)
                                      cGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATcc                                                            <  2:550678/100‑1 (MQ=255)
                                        aGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTTCTGTGTGATcccc                                                          <  1:1626319/100‑3 (MQ=255)
                                           cTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCActgc                                                       <  1:533430/100‑1 (MQ=255)
                                           cTGTTCCCGTGGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCActgc                                                       <  1:1071048/100‑1 (MQ=255)
                                            tGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGGTGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCActgc                                                       <  1:373674/99‑1 (MQ=255)
                                            tGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCActgct                                                      >  2:1112984/1‑100 (MQ=255)
                                             gTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCActgctg                                                     >  2:1820527/1‑100 (MQ=255)
                                              ttCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGa                                                    <  1:303777/100‑1 (MQ=255)
                                              ttCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGa                                                    <  2:531158/100‑1 (MQ=255)
                                                cccGTTGATGTGTCGAGCACTCTACGCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGaa                                                   <  2:830789/99‑1 (MQ=255)
                                                  cGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAatat                                                >  2:371495/1‑100 (MQ=255)
                                                    ttGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATcc                                              <  1:1594223/100‑1 (MQ=255)
                                                        ggTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATccgccg                                          >  2:2103705/1‑100 (MQ=255)
                                                        ggTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATccgccg                                          >  2:1803273/1‑100 (MQ=255)
                                                        ggTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATccgccg                                          >  2:1377134/1‑100 (MQ=255)
                                                        ggTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATccgccg                                          >  1:1235392/1‑100 (MQ=255)
                                                         gtgtCGAGCACTCTACTCTCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGt                                         <  2:2355819/100‑1 (MQ=255)
                                                         gtgtCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGt                                         <  1:1936742/100‑1 (MQ=255)
                                                          tgtCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTa                                        >  1:2173451/1‑100 (MQ=255)
                                                          tgtCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTa                                        >  1:570824/1‑100 (MQ=255)
                                                              gAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGa                                    >  2:2058581/1‑100 (MQ=255)
                                                                gCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACg                                  >  1:629501/1‑100 (MQ=255)
                                                                 cACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGt                                 >  2:1869527/1‑100 (MQ=255)
                                                                   ctctACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGtt                                <  1:1369156/99‑1 (MQ=255)
                                                                     ctACTCCCTGCTGATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAAt                             >  2:816419/1‑100 (MQ=255)
                                                                     ctACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAAt                             >  1:997956/1‑100 (MQ=255)
                                                                       aCTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATcc                           >  2:609935/1‑100 (MQ=255)
                                                                        cTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCg                          >  2:513302/1‑100 (MQ=255)
                                                                          cccTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGct                        <  1:223737/100‑1 (MQ=255)
                                                                           ccTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGctc                       >  2:857875/1‑100 (MQ=255)
                                                                           ccTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGctc                       >  1:2192849/1‑100 (MQ=255)
                                                                            ctgctgATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGctct                      <  1:2392224/100‑1 (MQ=255)
                                                                            ctgctgAGAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGctct                      <  2:1934619/100‑1 (MQ=255)
                                                                              gctgATGGTAGGGTGGGTGGCGGACGATGAATTCGATGCCGTCTGTCGCCTGGGTGATCCACTGCTGAATATCAGCCGTATGGCCGTTATTCCGctc                       <  2:1493479/97‑1 (MQ=255)
                                                                              gctgATAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGc                    <  1:1003284/100‑1 (MQ=255)
                                                                              gctgATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGc                    <  1:1112994/100‑1 (MQ=255)
                                                                              gctgATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGc                    >  2:1327875/1‑100 (MQ=255)
                                                                                  aTAGTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCAc                >  1:2298651/1‑100 (MQ=255)
                                                                                     gTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACggg             >  1:1597607/1‑100 (MQ=255)
                                                                                     gTAAGGTGGGTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACggg             >  1:394260/1‑100 (MQ=255)
                                                                                       aaGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGa           >  2:1704036/1‑100 (MQ=255)
                                                                                        aGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGAc          <  2:2010621/100‑1 (MQ=255)
                                                                                         ggTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGAc          <  1:1485495/99‑1 (MQ=255)
                                                                                            gggTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGACTGGt      >  1:1391099/1‑100 (MQ=255)
                                                                                             ggTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGACTGGtt     >  1:1448115/1‑100 (MQ=255)
                                                                                             ggTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGACTGGtt     >  2:133677/1‑100 (MQ=255)
                                                                                              gTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGACTGGTTg    >  1:1145869/1‑100 (MQ=255)
                                                                                              gTTGCTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGACTGGTTg    <  2:81836/100‑1 (MQ=255)
                                                                                              gTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCCAggggtgaggtt     >  1:2246201/1‑92 (MQ=255)
                                                                                              gTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGACTGGTTg    >  1:2139757/1‑100 (MQ=255)
                                                                                                tACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGACTGGTTGat  >  2:2118785/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
AGGTCAGCAACGATGGGTGATCAAGCGTAGTGGTATAGCGAGCCTGTTCCCGTTGAGGTGTCGAGCACTCTACTCCCTGCTGATAGTAAGGTGGGTTACTGATGATTAAATCGAAGCGTACTGTCTGCTGTGTGATCCACTGCTGAATATCCGCCGTATGGACGTTAATCCGCTCTGCCCACGGGGACTGGTTGAT  >  NC_000913/2712305‑2712500

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: