Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,193,289 T→A 24.0% 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,2890TA24.0% 246.5 / 76.6 125Q162L (CAG→CTG) mnmAtRNA(Gln,Lys,Glu) U34 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base T (47/48);  new base A (13/17);  total (60/65)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACGT  >  NC_000913/1193192‑1193388
                                                                                                 |                                                                                                   
gTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGc                                                                                                   >  2:1184993/1‑100 (MQ=255)
  ctcttttCTGCGCGGTGACCACCCACAGATGCTCCGCAATCTTACGCACCCGCGGTTCTTCCACTTCGCCGACCGGGAACAGGCTTTGTTCAGTCTGCTc                                                                                                 <  2:1803451/97‑1 (MQ=255)
  cTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTc                                                                                                 >  1:2957903/1‑100 (MQ=255)
   ttttttCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCa                                                                                                <  1:3247262/100‑1 (MQ=255)
    tttttCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCAt                                                                                               >  1:2658483/1‑100 (MQ=255)
       ttcttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGc                                                                                            <  2:2735989/100‑1 (MQ=255)
       ttcttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGc                                                                                            <  2:1784886/100‑1 (MQ=255)
        tcttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCt                                                                                           >  2:49089/1‑100 (MQ=255)
         cttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTCCGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTg                                                                                          <  1:2932745/100‑1 (MQ=255)
         cttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCt                                                                                           >  2:2983280/1‑99 (MQ=255)
         cttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTACAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTg                                                                                          <  2:3055736/100‑1 (MQ=255)
          ttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGa                                                                                         <  1:2172097/100‑1 (MQ=255)
             gcgGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCg                                                                                      >  1:2372426/1‑100 (MQ=255)
             gcgGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCg                                                                                      <  1:2171490/100‑1 (MQ=255)
               ggTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGta                                                                                    <  1:4087826/100‑1 (MQ=255)
                gTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGtat                                                                                   >  1:2365441/1‑100 (MQ=255)
                  gaccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGtatca                                                                                 >  2:2531624/1‑98 (MQ=255)
                    ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAg                                                                               <  1:3817824/100‑1 (MQ=255)
                    ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAg                                                                               <  2:1977565/100‑1 (MQ=255)
                    ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAg                                                                               <  2:2047555/100‑1 (MQ=255)
                    ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAg                                                                               <  1:760026/100‑1 (MQ=255)
                      agaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGa                                                                             <  1:3464767/100‑1 (MQ=255)
                        accaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAg                                                                           <  1:166596/100‑1 (MQ=255)
                         ccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGt                                                                          <  1:200913/100‑1 (MQ=255)
                         ccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGt                                                                          <  2:2230233/100‑1 (MQ=255)
                           aaGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAg                                                                        <  1:325536/100‑1 (MQ=255)
                            aGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGc                                                                       <  1:2113895/100‑1 (MQ=255)
                              aTCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATCACTGAGCGTATAAAGGCAGTAGctc                                                                     >  2:3920095/1‑100 (MQ=255)
                                 ttCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTgg                                                                  <  1:1804738/100‑1 (MQ=255)
                                 ttCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTgg                                                                  <  2:378299/100‑1 (MQ=255)
                                  tCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGt                                                                 >  2:2066264/1‑100 (MQ=255)
                                   cAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTc                                                                <  2:1055694/100‑1 (MQ=255)
                                   cAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTggcc                                                                >  1:3098997/1‑98 (MQ=255)
                                   cAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTc                                                                <  2:2869663/100‑1 (MQ=255)
                                    aGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCt                                                               >  2:3581757/1‑100 (MQ=255)
                                    aGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCt                                                               >  2:368749/1‑100 (MQ=255)
                                    aGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTc                                                                >  2:3713110/1‑99 (MQ=255)
                                     gCAATCTTACGCACCTGCGGTTTTTTTTTTTTTTTTACAGGGAACAGGCGTTGCGCAATTTTCTCGTGGCTGGGCGTATAAAGGAAGTAGCTCTGGTCtt                                                              >  2:1213790/1‑100 (MQ=255)
                                     gCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCtt                                                              <  2:2891948/100‑1 (MQ=255)
                                       aaTCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTa                                                            <  1:466850/100‑1 (MQ=255)
                                         tCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTAtt                                                          <  1:946328/100‑1 (MQ=255)
                                          cTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTg                                                         >  1:3714019/1‑100 (MQ=255)
                                          cTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTg                                                         <  1:1484715/100‑1 (MQ=255)
                                          cTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTg                                                         <  2:842665/100‑1 (MQ=255)
                                           ttACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATTAAGGAAGTAGCTCTGGTCTTTATTGc                                                        <  2:3769515/100‑1 (MQ=255)
                                           ttACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGGAGTAGCTCTGGTCTTTATTGc                                                        >  2:1362631/1‑100 (MQ=255)
                                           ttACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGc                                                        <  1:612658/100‑1 (MQ=255)
                                           ttACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGc                                                        <  1:525067/100‑1 (MQ=255)
                                             aCGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTg                                                      >  1:1255930/1‑100 (MQ=255)
                                              cGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGt                                                     <  1:368752/100‑1 (MQ=255)
                                              cGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGt                                                     >  1:4092326/1‑100 (MQ=255)
                                              cGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGt                                                     >  2:396428/1‑100 (MQ=255)
                                               gCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTc                                                    >  2:3073961/1‑100 (MQ=255)
                                                cACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTcc                                                   <  1:1017674/100‑1 (MQ=255)
                                                cACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTcc                                                   <  1:3847544/100‑1 (MQ=255)
                                                   cTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGa                                                <  2:1648414/100‑1 (MQ=255)
                                                    tGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGAc                                               >  2:3317279/1‑100 (MQ=255)
                                                     gCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGAcc                                              <  1:2234217/100‑1 (MQ=255)
                                                       ggTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCAc                                            <  1:1334074/100‑1 (MQ=255)
                                                         tttttCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgc                                          <  2:1659225/100‑1 (MQ=255)
                                                             tCCAGTTCGCCGACCGGGAACAGGCTTTGCGCACTCTGCTCATGGCAGAGCGTAGGAAAGAAGTAGCTCCGGTCTGTATTGCTGTCCAGACCACgc                                          >  2:2572172/1‑96 (MQ=255)
                                                             tCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcagca                                      >  1:965716/1‑100 (MQ=255)
                                                              ccAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCCTGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGcagcag                                     >  2:457463/1‑100 (MQ=255)
                                                               cAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAgg                                    >  1:2247306/1‑100 (MQ=255)
                                                               cAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAgg                                    >  1:3017681/1‑100 (MQ=255)
                                                               cAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAgg                                    >  2:4157341/1‑100 (MQ=255)
                                                               cAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAgg                                    >  1:1745862/1‑100 (MQ=255)
                                                                 gTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcg                                  >  1:1028044/1‑100 (MQ=255)
                                                                 gTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcg                                  >  1:265864/1‑100 (MQ=255)
                                                                 gTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcg                                  <  2:3082918/100‑1 (MQ=255)
                                                                  ttCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcgg                                 >  2:3200995/1‑100 (MQ=255)
                                                                   tCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcggc                                <  2:1186682/100‑1 (MQ=255)
                                                                   tCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcggc                                <  2:1967050/100‑1 (MQ=255)
                                                                   tCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcggc                                <  2:1575914/100‑1 (MQ=255)
                                                                   tCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcggc                                <  1:4157375/100‑1 (MQ=255)
                                                                   tCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcgg                                 >  2:893048/1‑99 (MQ=255)
                                                                    caccGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGct                               >  1:2819950/3‑100 (MQ=255)
                                                                    cGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGct                               >  1:1981196/1‑100 (MQ=255)
                                                                    cGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGct                               >  2:1168696/1‑100 (MQ=255)
                                                                    cGCCGACCGGGAACAGGCTTTGCGCAATAAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGct                               >  1:2069200/1‑100 (MQ=255)
                                                                     gCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGctc                              <  2:4092294/100‑1 (MQ=255)
                                                                      ccGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGctct                             >  1:4082702/1‑100 (MQ=255)
                                                                        gACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTg                           <  2:1524293/100‑1 (MQ=255)
                                                                         aCCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGc                          <  2:3943009/100‑1 (MQ=255)
                                                                          ccGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGcc                         >  1:3714630/1‑100 (MQ=255)
                                                                          ccGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGcc                         >  2:1941492/1‑100 (MQ=255)
                                                                              gAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCg                     <  1:1048797/100‑1 (MQ=255)
                                                                              gAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCg                     >  1:3922782/1‑100 (MQ=255)
                                                                                aCAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGAc                   <  2:856244/100‑1 (MQ=255)
                                                                                aCAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGAc                   <  1:3364312/100‑1 (MQ=255)
                                                                                  aGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACAt                 >  2:3653414/1‑100 (MQ=255)
                                                                                   ggCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATc                >  2:2322619/1‑100 (MQ=255)
                                                                                   ggCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATc                >  2:247762/1‑100 (MQ=255)
                                                                                   gccTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATc                <  1:2147749/98‑1 (MQ=255)
                                                                                       ttGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGcc            <  1:439122/100‑1 (MQ=255)
                                                                                        tGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGccc           <  1:920084/100‑1 (MQ=255)
                                                                                        tGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGcc            <  1:3908084/99‑1 (MQ=255)
                                                                                         gcgcAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCg          >  1:649652/1‑100 (MQ=255)
                                                                                             aaTCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGacgn      <  1:1533413/100‑2 (MQ=255)
                                                                                                 tGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGacgtacgt  >  1:1064028/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACGT  >  NC_000913/1193192‑1193388

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