Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,193,289 T→A Q162L (CAG→CTG)  mnmA ← tRNA(Gln,Lys,Glu) U34 2‑thiouridylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,193,2890TA96.6% 300.2 / 1.6 91Q162L (CAG→CTG) mnmAtRNA(Gln,Lys,Glu) U34 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base A (40/46);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

TCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTA  >  NC_000913/1193193‑1193385
                                                                                                |                                                                                                
tCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCt                                                                                               <  2:159714/100‑1 (MQ=255)
tCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCt                                                                                               <  1:1570174/100‑1 (MQ=255)
  ttttttCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCa                                                                                             >  1:280471/1‑100 (MQ=255)
    ttttCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATg                                                                                           >  1:1117527/1‑100 (MQ=255)
     tttCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCGGCTCATgg                                                                                          >  2:677285/1‑100 (MQ=255)
      ttcttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGc                                                                                         <  1:1416650/100‑1 (MQ=255)
      ttcttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGc                                                                                         >  1:1775855/1‑100 (MQ=255)
        cttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTg                                                                                       <  1:1919150/100‑1 (MQ=255)
          tcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAg                                                                                     >  1:1900213/1‑100 (MQ=255)
          tcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAg                                                                                     <  2:1859422/100‑1 (MQ=255)
            gcgGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCg                                                                                   >  2:908506/1‑100 (MQ=255)
             cgGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGt                                                                                  >  1:674754/1‑100 (MQ=255)
              ggTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGta                                                                                 >  2:408646/1‑100 (MQ=255)
               gTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGtat                                                                                <  2:674753/100‑1 (MQ=255)
                tGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGtata                                                                               <  1:1639153/100‑1 (MQ=255)
                tGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGtata                                                                               >  1:1694009/1‑100 (MQ=255)
                 gaccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATaa                                                                              >  2:1103442/1‑100 (MQ=255)
                 gaccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGAAACAGGATTTGCGCAATCAGCTCATGGCTGAGCGTATaa                                                                              >  2:1639148/1‑100 (MQ=255)
                  accagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATaaa                                                                             >  1:1153776/1‑100 (MQ=255)
                   ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTTTTGTTTTCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAg                                                                            >  2:124139/1‑100 (MQ=255)
                   ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCGGCTCATGGCTGAGCGTATAAAg                                                                            <  1:677286/100‑1 (MQ=255)
                   ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAg                                                                            <  2:1024049/100‑1 (MQ=255)
                   ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAg                                                                            >  2:1914279/1‑100 (MQ=255)
                     agaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGGGCGTATAAAGGa                                                                          >  1:1177323/1‑100 (MQ=255)
                     agaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGa                                                                          <  2:1882260/100‑1 (MQ=255)
                     agaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGa                                                                          <  1:1260342/100‑1 (MQ=255)
                     agaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGa                                                                          <  2:435351/100‑1 (MQ=255)
                      gaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGaa                                                                         <  1:2082865/100‑1 (MQ=255)
                       tccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAg                                                                        <  1:326877/99‑1 (MQ=255)
                        ccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAga                                                                       >  2:326875/1‑99 (MQ=255)
                        ccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGt                                                                       <  2:1085444/100‑1 (MQ=255)
                        ccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGt                                                                       >  2:1293349/1‑100 (MQ=255)
                         caAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTa                                                                      <  2:1900206/100‑1 (MQ=255)
                         caAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTa                                                                      <  2:1784548/100‑1 (MQ=255)
                           aGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAgn                                                                    <  1:1930906/100‑2 (MQ=255)
                           aGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGc                                                                    <  1:219773/100‑1 (MQ=255)
                            gATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGct                                                                   <  2:311256/100‑1 (MQ=255)
                               gatctgCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTg                                                                <  2:298950/95‑1 (MQ=255)
                               cTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTg                                                                <  2:2196469/100‑1 (MQ=255)
                               cTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGATCATGGCTGAGCGTATAAAGGAAGTAGCTCTg                                                                >  2:1644710/1‑100 (MQ=255)
                                ttCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTgg                                                               <  1:1978638/100‑1 (MQ=255)
                                  cAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTc                                                             >  2:1203417/1‑100 (MQ=255)
                                   aGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCt                                                            <  2:1954050/100‑1 (MQ=255)
                                    gCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTgagatc                                                           >  1:298952/1‑95 (MQ=255)
                                    gCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCtt                                                           <  1:124141/100‑1 (MQ=255)
                                    gCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCACAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCtt                                                           >  2:1665548/1‑100 (MQ=255)
                                     cAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCttt                                                          >  1:157999/1‑100 (MQ=255)
                                        tCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTAtt                                                       <  2:821939/100‑1 (MQ=255)
                                          ttACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGGGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGc                                                     <  2:1177320/100‑1 (MQ=255)
                                          ttACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGc                                                     <  1:2096289/100‑1 (MQ=255)
                                           tACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGc                                                     <  2:273256/99‑1 (MQ=255)
                                           tACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCt                                                    >  2:2239754/1‑100 (MQ=255)
                                            aCGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTg                                                   <  1:374068/100‑1 (MQ=255)
                                            aCGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTg                                                   >  1:273843/1‑100 (MQ=255)
                                            aCGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTg                                                   <  1:609279/100‑1 (MQ=255)
                                               cACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTcc                                                >  2:1861409/1‑100 (MQ=255)
                                                 ccTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAg                                              <  1:34892/100‑1 (MQ=255)
                                                   tGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGAc                                            >  1:1685215/1‑100 (MQ=255)
                                                    gCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGAcc                                           >  2:11450/1‑100 (MQ=255)
                                                        gatctCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgc                                       <  2:1806663/96‑1 (MQ=255)
                                                         ttttCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgca                                      >  1:1129839/1‑100 (MQ=255)
                                                         ttttCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgca                                      >  1:1960917/1‑100 (MQ=255)
                                                         ttttCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgca                                      <  2:1775850/100‑1 (MQ=255)
                                                         ttttCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgca                                      <  2:208461/100‑1 (MQ=255)
                                                            tCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcagca                                   >  1:664899/1‑100 (MQ=255)
                                                            tCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcagca                                   >  2:2022829/1‑100 (MQ=255)
                                                             ccAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcagatc                                  >  1:1806668/1‑97 (MQ=255)
                                                             ccAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGcagcag                                  >  2:1150896/1‑100 (MQ=255)
                                                              cAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAgg                                 >  2:423971/1‑100 (MQ=255)
                                                               aGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggc                                <  1:1480670/100‑1 (MQ=255)
                                                               aGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggc                                >  1:2059565/1‑100 (MQ=255)
                                                                gTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcg                               >  2:796818/1‑100 (MQ=255)
                                                                  tCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcggc                             <  1:1914286/100‑1 (MQ=255)
                                                                  tCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcggc                             <  2:1624241/100‑1 (MQ=255)
                                                                     ccGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGctct                          <  2:1153774/100‑1 (MQ=255)
                                                                     ccGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGctct                          >  2:1763783/1‑100 (MQ=255)
                                                                      cGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCtt                         >  1:500988/1‑100 (MQ=255)
                                                                        aCCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGAAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGc                       <  1:394105/100‑1 (MQ=255)
                                                                           gggAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCAt                    <  1:908508/100‑1 (MQ=255)
                                                                                  ggCTTTGCGCAATCAGCTCATGTCTGAGCGTGTAAAGTAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATc             <  2:1579869/100‑1 (MQ=255)
                                                                                   gCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCg            >  1:1544611/1‑100 (MQ=255)
                                                                                   gCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCg            >  1:1752204/1‑100 (MQ=255)
                                                                                    cTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCgg           <  1:1103444/100‑1 (MQ=255)
                                                                                     tttGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGc          <  1:408648/100‑1 (MQ=255)
                                                                                      ttGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGcc         <  1:2093131/100‑1 (MQ=255)
                                                                                       tGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGccc        >  2:533153/1‑100 (MQ=255)
                                                                                        gcgcAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCg       <  2:172216/100‑1 (MQ=255)
                                                                                          gcAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGac     >  1:2119449/1‑100 (MQ=255)
                                                                                          gcAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGac     <  2:273841/100‑1 (MQ=255)
                                                                                           cAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGacg    <  2:280469/100‑1 (MQ=255)
                                                                                             aTCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGacgta  >  1:781195/1‑100 (MQ=255)
                                                                                                |                                                                                                
TCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTA  >  NC_000913/1193193‑1193385

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