Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,143,641 G→A 100% H243Y (CAT→TAT)  rne ← endoribonuclease; RNA‑binding protein;RNA degradosome binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,143,6410GA98.9% 310.0 / ‑4.5 89H243Y (CAT→TAT) rneendoribonuclease; RNA‑binding protein;RNA degradosome binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (41/47);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

ATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGC  >  NC_000913/1143544‑1143736
                                                                                                 |                                                                                               
aTCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACt                                                                                               <  1:361050/100‑1 (MQ=255)
  cTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCTGGGCGACCTAATGCAGCGATATACTGa                                                                                             >  2:988694/1‑100 (MQ=255)
    gtgACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACg                                                                                           <  1:60925/100‑1 (MQ=255)
    gtgACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACg                                                                                           <  2:121178/100‑1 (MQ=255)
    gtgACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACg                                                                                           >  1:347138/1‑100 (MQ=255)
    gtgACTCGATCTGGTAGTGGCTGAACAGCAGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGACATCCGGGCGACCTAATGCAGCGATATACTGACg                                                                                           <  1:1832594/100‑1 (MQ=255)
     tgACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACg                                                                                           >  1:1792577/1‑99 (MQ=255)
     tgACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGt                                                                                          >  2:57647/1‑100 (MQ=255)
      gACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTg                                                                                         >  1:1469955/1‑100 (MQ=255)
      gACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTg                                                                                         >  1:1942459/1‑100 (MQ=255)
          cGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAg                                                                                     <  1:89292/100‑1 (MQ=255)
             tCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTc                                                                                  <  2:860581/100‑1 (MQ=255)
              cTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCg                                                                                 <  1:169271/100‑1 (MQ=255)
                ggTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAg                                                                               <  2:1538085/100‑1 (MQ=255)
                ggTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAg                                                                               >  2:2307440/1‑100 (MQ=255)
                     tGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACtt                                                                          <  1:366803/100‑1 (MQ=255)
                      ggCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACttt                                                                         >  1:344271/1‑100 (MQ=255)
                      ggCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACttt                                                                         >  2:2430216/1‑100 (MQ=255)
                      ggCTGAACAGCGGGAGCTCGCCGGTGTACAGTGTGATTTGGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTGCGGGCACttt                                                                         >  1:1825225/1‑100 (MQ=255)
                         tGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGCCCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCgg                                                                      >  2:1721467/1‑100 (MQ=255)
                         tGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCgg                                                                      <  1:1154458/100‑1 (MQ=255)
                           aaCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGt                                                                    <  2:1176630/100‑1 (MQ=255)
                           aaCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGt                                                                    <  1:1028755/100‑1 (MQ=255)
                            aCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGtt                                                                   <  2:2619058/100‑1 (MQ=255)
                            aCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGtt                                                                   >  1:963606/1‑100 (MQ=255)
                             cAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCCGTTCGAGCACTTTCGGGTTa                                                                  <  1:1872324/100‑1 (MQ=255)
                             cAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTa                                                                  >  2:387255/1‑100 (MQ=255)
                             cAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTa                                                                  >  2:556934/1‑100 (MQ=255)
                              aGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTAt                                                                 <  2:1629699/100‑1 (MQ=255)
                               gCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATc                                                                >  2:379941/1‑100 (MQ=255)
                               gCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATc                                                                <  2:2275033/100‑1 (MQ=255)
                               gCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATc                                                                <  1:2036904/100‑1 (MQ=255)
                                cGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCg                                                               <  1:2350378/100‑1 (MQ=255)
                                cGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCg                                                               <  2:350719/100‑1 (MQ=255)
                                cGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCg                                                               <  1:57646/100‑1 (MQ=255)
                                cGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCg                                                               <  1:2180787/100‑1 (MQ=255)
                                 gggATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGa                                                              <  2:1173075/100‑1 (MQ=255)
                                   gATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATa                                                            <  1:826575/100‑1 (MQ=255)
                                     tctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAg                                                          >  2:665327/1‑100 (MQ=255)
                                      ctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAgg                                                         <  2:2360260/100‑1 (MQ=255)
                                      ctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAgg                                                         <  2:2742753/100‑1 (MQ=255)
                                         gCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGAtt                                                      <  1:37377/100‑1 (MQ=255)
                                         gCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGAtt                                                      >  1:283948/1‑100 (MQ=255)
                                            ggTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCg                                                   <  1:1725391/100‑1 (MQ=255)
                                               gtACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGcc                                                 >  1:1441266/1‑99 (MQ=255)
                                               gtACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCg                                                <  1:12611/100‑1 (MQ=255)
                                               gtACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCg                                                <  2:2972298/100‑1 (MQ=255)
                                               gtACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCg                                                <  1:2440957/100‑1 (MQ=255)
                                               gtACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCg                                                <  1:1411239/100‑1 (MQ=255)
                                                tACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTCTCGATCAGGATTTCGCCg                                                >  1:2412494/1‑99 (MQ=255)
                                                tACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGa                                               <  1:1112971/100‑1 (MQ=255)
                                                tACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGa                                               >  1:2665728/1‑100 (MQ=255)
                                                tACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGa                                               <  1:1390568/100‑1 (MQ=255)
                                                tACAGTTTGAGTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGa                                               <  1:84829/100‑1 (MQ=255)
                                                  cAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATg                                             <  1:2212566/100‑1 (MQ=255)
                                                   aGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGt                                            >  2:1731547/1‑100 (MQ=255)
                                                    gTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTc                                           >  2:2820493/1‑100 (MQ=255)
                                                    gTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTc                                           <  2:499779/100‑1 (MQ=255)
                                                       tGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTCCGAGCACTTTCCGGGTATCGATAAGGATTTCGCCGATGTCCt                                         >  1:2434489/1‑99 (MQ=255)
                                                         aTTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGAc                                      <  2:2263586/100‑1 (MQ=255)
                                                              gctgctGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaag                                 <  1:2511844/100‑1 (MQ=255)
                                                               ctgctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaagt                                >  1:500673/1‑100 (MQ=255)
                                                                 gctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaagtaa                              >  2:2135102/1‑100 (MQ=255)
                                                                 gctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaagtaa                              >  2:2025512/1‑100 (MQ=255)
                                                                 gctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaagtaa                              <  1:430894/100‑1 (MQ=255)
                                                                 gctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaagtaa                              >  2:3057810/1‑100 (MQ=255)
                                                                   tgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATc                            <  1:2269391/100‑1 (MQ=255)
                                                                   tgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATc                            >  1:2854527/1‑100 (MQ=255)
                                                                   tgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATc                            >  2:1585743/1‑100 (MQ=255)
                                                                      aaTCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATcgcg                         >  2:1217436/1‑100 (MQ=255)
                                                                        tCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGa                       >  2:672641/1‑100 (MQ=255)
                                                                             gCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATgcg                  <  1:379937/100‑1 (MQ=255)
                                                                              cGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATgcgc                 <  1:2134765/100‑1 (MQ=255)
                                                                               gACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGcgcg                <  1:2800770/100‑1 (MQ=255)
                                                                               gACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGcgcg                >  1:2856077/1‑100 (MQ=255)
                                                                               gACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGACTTCGCCGATGTCCTGACGTAAGTAATCGCGGAAGGcgcg                <  2:1261230/100‑1 (MQ=255)
                                                                                  cTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAAc             >  1:1850465/1‑100 (MQ=255)
                                                                                   tAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACg            >  2:490027/1‑100 (MQ=255)
                                                                                       gCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTa        <  2:2240649/100‑1 (MQ=255)
                                                                                       gCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTa        >  1:517321/1‑100 (MQ=255)
                                                                                       gCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTa        >  2:1425785/1‑100 (MQ=255)
                                                                                        cAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTAc       >  2:2953067/1‑100 (MQ=255)
                                                                                        cAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTAc       >  1:2791716/1‑100 (MQ=255)
                                                                                        cAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTGCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTAc       >  2:468171/1‑100 (MQ=255)
                                                                                        cAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGATTATCGATAAGGGTTTCGCGGATGTCCTGGCGTAAGTAATCGCGGAATGCGCGAACGGtt         >  1:2709149/1‑98 (MQ=255)
                                                                                         aGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACg      >  1:649156/1‑100 (MQ=255)
                                                                                           cGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGtt    <  2:542215/100‑1 (MQ=255)
                                                                                             atatACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc  <  1:2342953/100‑1 (MQ=255)
                                                                                             atatACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc  <  2:1469961/100‑1 (MQ=255)
                                                                                                 |                                                                                               
ATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGC  >  NC_000913/1143544‑1143736

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