Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,799,488 C→T 18.1% V88I (GTT→ATT)  rplT ← 50S ribosomal subunit protein L20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,799,4880CT18.1% 202.4 / 34.8 90V88I (GTT→ATT) rplT50S ribosomal subunit protein L20
Reads supporting (aligned to +/- strand):  ref base C (48/25);  new base T (10/6);  total (59/31)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

TGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGAACTGAC  >  NC_000913/1799389‑1799586
                                                                                                   |                                                                                                  
tGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAAc                                                                                                    <  2:1475816/100‑1 (MQ=255)
   cTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACaga                                                                                                 >  2:420439/1‑100 (MQ=255)
    ttATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATagag                                                                                                <  1:547178/100‑1 (MQ=255)
     tATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAgg                                                                                               <  1:376508/100‑1 (MQ=255)
      aTGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGACTCCTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGc                                                                                              >  2:664050/1‑100 (MQ=255)
       tGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGc                                                                                              >  2:1504852/1‑99 (MQ=255)
       tGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCt                                                                                             >  1:697264/1‑100 (MQ=255)
       tGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCt                                                                                             >  2:53733/1‑100 (MQ=255)
         ccAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTCCTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGGCGATTTCAAGAGCGGCttt                                                                                           >  1:1527447/1‑100 (MQ=255)
          cAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCtttt                                                                                          >  2:62689/1‑100 (MQ=255)
          cAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCtttt                                                                                          >  2:970717/1‑100 (MQ=255)
            gagCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCtttttt                                                                                        >  1:1452591/1‑100 (MQ=255)
              gctgctTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCa                                                                                      <  1:1254070/100‑1 (MQ=255)
               ctgctTTCGCTTTTTCAACCAGAGCGGTGAACGCTAATTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg                                                                                     >  1:1658842/1‑100 (MQ=255)
                 gctTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGc                                                                                   >  1:104842/1‑100 (MQ=255)
                 gctTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGc                                                                                   >  1:1809362/1‑100 (MQ=255)
                   tttCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCa                                                                                 >  2:983715/1‑100 (MQ=255)
                   tttCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCa                                                                                 >  1:1257457/1‑100 (MQ=255)
                     tCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCAtt                                                                               >  2:1817282/1‑100 (MQ=255)
                      cGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg                                                                              >  2:1265292/1‑100 (MQ=255)
                      cGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg                                                                              >  2:2244689/1‑100 (MQ=255)
                      cGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATAAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTACCGGCCATtg                                                                              >  2:354904/1‑100 (MQ=255)
                       gCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtga                                                                             >  1:333223/1‑100 (MQ=255)
                       gCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg                                                                              <  2:2258439/99‑1 (MQ=255)
                        cTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtgat                                                                            <  2:1679254/100‑1 (MQ=255)
                             tCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAAtt                                                                       >  1:292129/1‑100 (MQ=255)
                              cAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAAttt                                                                      >  1:1632941/1‑100 (MQ=255)
                                aCCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGc                                                                    >  2:1211373/1‑100 (MQ=255)
                                aCCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGc                                                                    <  2:2081666/100‑1 (MQ=255)
                                aCCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGc                                                                    <  1:402210/100‑1 (MQ=255)
                                aCCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGc                                                                    <  2:1090221/100‑1 (MQ=255)
                                 ccAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCt                                                                   >  1:1971565/1‑100 (MQ=255)
                                 ccAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCt                                                                   <  1:2329940/100‑1 (MQ=255)
                                  cAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTg                                                                  <  1:399462/100‑1 (MQ=255)
                                     agCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTaa                                                               <  2:665004/100‑1 (MQ=255)
                                        ggTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGaa                                                            >  2:1116651/1‑100 (MQ=255)
                                          tGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAAt                                                          <  1:1343105/100‑1 (MQ=255)
                                           gAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATa                                                         >  2:825313/1‑100 (MQ=255)
                                           gAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAAt                                                          <  2:2148111/99‑1 (MQ=255)
                                           cttccgatcTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATa                                                         <  1:2075188/92‑1 (MQ=255)
                                             aCGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATAcc                                                       >  1:652049/1‑100 (MQ=255)
                                                cTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGtt                                                    >  2:830571/1‑100 (MQ=255)
                                                 tACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTc                                                   >  2:1224400/1‑100 (MQ=255)
                                                   cTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTg                                                 <  1:1265297/100‑1 (MQ=255)
                                                     ttGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATaagatcggaag                                               >  2:2075171/1‑90 (MQ=255)
                                                      tGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACg                                              <  1:1100963/100‑1 (MQ=255)
                                                      tGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACg                                              <  1:1393792/100‑1 (MQ=255)
                                                      tGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACg                                              <  2:997892/100‑1 (MQ=255)
                                                      tGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACg                                              <  1:2084320/100‑1 (MQ=255)
                                                       gTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGt                                             >  1:1810737/1‑100 (MQ=255)
                                                          gAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGtgct                                          >  1:295372/1‑100 (MQ=255)
                                                          gAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGtgct                                          >  1:2019031/1‑100 (MQ=255)
                                                          gAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGtgct                                          >  1:1849367/1‑100 (MQ=255)
                                                            aTACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGtgctgc                                        <  1:808559/100‑1 (MQ=255)
                                                                tGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCg                                    >  1:467432/1‑100 (MQ=255)
                                                                  cGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGtt                                  <  1:2119131/100‑1 (MQ=255)
                                                                   gATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTg                                 <  2:742835/100‑1 (MQ=255)
                                                                   gATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTg                                 <  1:535519/100‑1 (MQ=255)
                                                                     tatCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGAt                               >  1:185481/1‑100 (MQ=255)
                                                                     tatCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTATAAGAAATACCGTTCTGACGTGCTGCTGCGTTGAt                               >  2:770296/1‑100 (MQ=255)
                                                                      atCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATa                              >  2:1195520/1‑100 (MQ=255)
                                                                      atCAGCCAGGAGCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATa                              <  2:2223375/100‑1 (MQ=255)
                                                                        cAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATAcg                            >  1:951228/1‑100 (MQ=255)
                                                                          gCCAGGATCTTACGGTCGATTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATAcgcg                          <  2:43583/100‑1 (MQ=255)
                                                                          gCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATAcgcg                          >  2:2268388/1‑100 (MQ=255)
                                                                          gCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATAcgcg                          >  1:1064872/1‑100 (MQ=255)
                                                                          gCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATAcgcg                          >  1:2089732/1‑100 (MQ=255)
                                                                            cAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTGTGACGTGCTGCTACGTTGATAcgcgaa                        >  2:326909/1‑98 (MQ=255)
                                                                            cAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCa                        >  1:2240469/1‑100 (MQ=255)
                                                                               gATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATc                     >  1:1936858/1‑100 (MQ=255)
                                                                                aTCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATcc                    >  1:1275376/1‑100 (MQ=255)
                                                                                  cTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCcac                  >  1:447909/1‑100 (MQ=255)
                                                                                  cTTACGGTCGATTTCAACAGAGGCTTGTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCcac                  <  2:1880231/100‑1 (MQ=255)
                                                                                   ttACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCcaca                 >  1:524918/1‑100 (MQ=255)
                                                                                    tACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAg                >  2:1756941/1‑100 (MQ=255)
                                                                                      cGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGtt              >  1:737827/1‑100 (MQ=255)
                                                                                       ggTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTg             >  2:2035854/1‑100 (MQ=255)
                                                                                          cGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACg          >  1:166303/1‑100 (MQ=255)
                                                                                          cGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACg          >  1:2298721/1‑100 (MQ=255)
                                                                                          cGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACg          >  1:2329972/1‑100 (MQ=255)
                                                                                           gATTTCAACAGAGGCTTTTTTTTTTTTTTTTATGAATTTGCTGTAAGAAATACCATTCTGACGTGCTGCTGCGTTTATACGCGCAATCCACAGTTGAc           >  2:1667856/1‑98 (MQ=255)
                                                                                           gATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACgg         <  2:107560/100‑1 (MQ=255)
                                                                                            aTTTCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGa        >  1:1337866/1‑100 (MQ=255)
                                                                                            aTTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGa        >  2:1539012/1‑100 (MQ=255)
                                                                                            aTTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGa        <  1:1299575/100‑1 (MQ=255)
                                                                                              ttCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGAAc      >  2:358549/1‑100 (MQ=255)
                                                                                               tCAATAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGAAc      >  2:1500321/1‑99 (MQ=255)
                                                                                               tCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGAACt     <  2:2190003/100‑1 (MQ=255)
                                                                                                cAACAGAGGCTTTTTCCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGAACTg    <  1:2317752/100‑1 (MQ=255)
                                                                                                 aaCAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATATCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGAACTGa   <  1:920301/100‑1 (MQ=255)
                                                                                                  aCAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGAACTGAc  >  2:1410901/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTTGCTGTAAGAAATACCGTTCTGACGTGCTGCTGCGTTGATACGCGCAATCCACAGTTGACGGAACTGAC  >  NC_000913/1799389‑1799586

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: