Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,193,289 T→A 54.8% 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,2890TA54.8% 33.9 / 215.8 155Q162L (CAG→CTG) mnmAtRNA(Gln,Lys,Glu) U34 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base T (29/41);  new base A (40/45);  total (69/86)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 155 total aligned reads displayed.

GTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  NC_000913/1193192‑1193387
                                                                                                 |                                                                                                  
gTCTTTTTTCTTCGCGGTGACCAGCCCAAGATCTTCAGCACTCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTg                                                                                                   <  1:424098/99‑1 (MQ=255)
gTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGc                                                                                                  <  1:3803569/100‑1 (MQ=255)
 tCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCt                                                                                                 <  1:2721489/100‑1 (MQ=255)
  cTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTc                                                                                                >  1:1605146/1‑100 (MQ=255)
       ttcttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGc                                                                                           <  1:2425408/100‑1 (MQ=255)
       ttcttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGc                                                                                           <  1:3253053/100‑1 (MQ=255)
       ttcttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGc                                                                                           <  2:734210/100‑1 (MQ=255)
        tcttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCt                                                                                          >  2:259816/1‑100 (MQ=255)
          ttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTGCAGTTCGCCGACCGGGAACAGGCTTTGAGCAATCAGCTCATGGCTGa                                                                                        >  1:3878664/1‑100 (MQ=255)
          ttcGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGa                                                                                        >  2:679797/1‑100 (MQ=255)
             gcgGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTGCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCg                                                                                     <  1:665937/100‑1 (MQ=255)
                 tGACCAGACCAAGATCTTCAGCAATCTTACGCTCCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGTAATCAGCTCATGGCTGAGCGtata                                                                                 <  1:438200/100‑1 (MQ=255)
                 tGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGtata                                                                                 <  2:2319200/100‑1 (MQ=255)
                 tGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGtata                                                                                 <  2:3892812/100‑1 (MQ=255)
                  gaccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATaa                                                                                >  2:4042146/1‑100 (MQ=255)
                   accagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGGGTATaaa                                                                               <  2:2240133/100‑1 (MQ=255)
                    ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAg                                                                              <  1:2969013/100‑1 (MQ=255)
                    ccagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTAGAAAg                                                                              >  1:2664332/1‑100 (MQ=255)
                     cagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATaaa                                                                               <  1:745060/98‑1 (MQ=255)
                     cagaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAgg                                                                             >  1:1703318/1‑100 (MQ=255)
                      agaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGa                                                                            <  2:1353069/100‑1 (MQ=255)
                      agaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGta                                                                                   <  1:3654332/93‑1 (MQ=255)
                      agaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGa                                                                            <  2:3553013/100‑1 (MQ=255)
                      agaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGa                                                                            <  1:397457/100‑1 (MQ=255)
                       gaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTTTTTTTTTTTTTCCGGGAACAGGCTTTGCGCAATCATCTTATGGCTGAGCGTATAAAGGa                                                                            >  1:3201831/1‑99 (MQ=255)
                       gaccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGaa                                                                           >  2:3980500/1‑100 (MQ=255)
                        accaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAg                                                                          <  2:2817538/100‑1 (MQ=255)
                        accaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAg                                                                          <  1:1884832/100‑1 (MQ=255)
                         ccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGt                                                                         >  2:3385823/1‑100 (MQ=255)
                         ccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGt                                                                         <  2:2421787/100‑1 (MQ=255)
                         ccaAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCCACCGGGACCCGGCTTTGCGCAACCTGCTCCTCGCTGAGCGTATCAAGGaact                                                                         >  1:3781023/1‑98 (MQ=255)
                          caAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTa                                                                        <  1:184447/100‑1 (MQ=255)
                           aaGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAg                                                                       >  2:505718/1‑100 (MQ=255)
                            aGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGc                                                                      <  1:909433/100‑1 (MQ=255)
                            aGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGc                                                                      <  1:1487910/100‑1 (MQ=255)
                            aGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGAATAAAGGAAGTAGc                                                                      >  2:3351160/1‑100 (MQ=255)
                             gATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGct                                                                     <  1:2358710/100‑1 (MQ=255)
                              aTCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGctc                                                                    <  1:1523024/100‑1 (MQ=255)
                              aTCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGctc                                                                    >  1:308232/1‑100 (MQ=255)
                               tCTTTAGCAAGTTTACGCACGTGCGGTTTTTCGAGTTCGCCGACCGGGATCAGGCGTTGCGCTATCTGCTCATGGCTGAGCGCCTCAAGGAAGTa                                                                        >  1:2916770/1‑95 (MQ=255)
                                 ttCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTgg                                                                 <  2:1039434/100‑1 (MQ=255)
                                    cgCAATCTTACGCATCTGCGGCATTTCCAGTTCGCCGACCAGGAACAGACTTTGCGCAATCTGCTCATGACTGAGCGTATAAAGAAAGTAACTATGGTCt                                                              >  1:703416/2‑100 (MQ=255)
                                     gCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCtt                                                             >  2:3350855/1‑100 (MQ=255)
                                     gCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCtt                                                             >  2:2434296/1‑100 (MQ=255)
                                       aaTCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTa                                                           <  1:1104147/100‑1 (MQ=255)
                                        aTCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTAt                                                          <  2:1079205/100‑1 (MQ=255)
                                        aTCTTACACACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTAt                                                          >  1:3955890/1‑100 (MQ=255)
                                           ttACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGc                                                       <  1:2160103/100‑1 (MQ=255)
                                            tACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTgat                                                      >  2:2207885/1‑98 (MQ=255)
                                              cGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGt                                                    >  2:17039/1‑100 (MQ=255)
                                                cACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTcc                                                  <  2:633527/100‑1 (MQ=255)
                                                cACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTcc                                                  <  1:1100513/100‑1 (MQ=255)
                                                  ccTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAg                                                <  2:3791998/100‑1 (MQ=255)
                                                    tGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGATCTGGTCTTTATTGCTGTCCAGAc                                              >  1:2103152/1‑100 (MQ=255)
                                                          ttttCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgca                                        >  2:2088137/1‑100 (MQ=255)
                                                           tttCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcag                                       >  2:214614/1‑100 (MQ=255)
                                                            ttCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcagc                                      <  1:2394221/100‑1 (MQ=255)
                                                            ttCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcagc                                      <  2:2083975/100‑1 (MQ=255)
                                                             tCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcagca                                     >  2:3275604/1‑100 (MQ=255)
                                                             tCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcagca                                     >  2:2974397/1‑100 (MQ=255)
                                                             tCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCGCATGGCTGAGCGTATAAAGGAAGTAGCGCTGGTCTTTATTGCTGTCCAGACCACgcag                                       <  1:3982696/98‑1 (MQ=255)
                                                              ccAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGcagcag                                    <  1:3322515/100‑1 (MQ=255)
                                                                aagtCGCCTACCGGGAACAGGCTTGGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggc                                  <  2:785932/97‑1 (MQ=255)
                                                                aGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggc                                  <  1:2678587/100‑1 (MQ=255)
                                                                aGTTCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggc                                  <  1:434356/100‑1 (MQ=255)
                                                                 gTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcg                                 <  2:3514138/100‑1 (MQ=255)
                                                                   tCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcggc                               >  2:700809/1‑100 (MQ=255)
                                                                   tCGCCGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcgac                               <  2:2886532/100‑3 (MQ=255)
                                                                    ngCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGCGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcggc                               >  1:2522459/2‑99 (MQ=255)
                                                                    cGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGct                              >  2:2150396/1‑100 (MQ=255)
                                                                    cGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGct                              >  2:1763597/1‑100 (MQ=255)
                                                                      ccGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGctct                            >  2:585811/1‑100 (MQ=255)
                                                                      ccGACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGctct                            <  2:293567/100‑1 (MQ=255)
                                                                        gACCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTg                          >  1:1306120/1‑100 (MQ=255)
                                                                         aCCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGc                         <  2:3619745/100‑1 (MQ=255)
                                                                         aCCGGGAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCAGCTCTTGc                         <  1:1779148/100‑1 (MQ=255)
                                                                            gggAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCttagatc                      >  2:3634464/1‑95 (MQ=255)
                                                                            gggAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCAt                      >  2:4166048/1‑100 (MQ=255)
                                                                             ggAACAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATc                     <  1:3053546/100‑1 (MQ=255)
                                                                              gAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCg                    <  1:3655171/100‑1 (MQ=255)
                                                                                aCAGGCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGAc                  >  1:3299679/1‑100 (MQ=255)
                                                                                aCAGGCTTTGCGCAATCAGCCCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGAc                  <  2:2103156/100‑1 (MQ=255)
                                                                                 cAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACa                 <  1:4204152/100‑1 (MQ=255)
                                                                                  aGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGTAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACAt                >  1:989167/1‑100 (MQ=255)
                                                                                   ggCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATc               >  1:2195042/1‑100 (MQ=255)
                                                                                   ggCTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATc               >  1:2853322/1‑100 (MQ=255)
                                                                                    gCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCg              >  1:1312452/1‑100 (MQ=255)
                                                                                     cTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCgg             >  2:3500936/1‑100 (MQ=255)
                                                                                     cTTTGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCgg             >  1:358820/1‑100 (MQ=255)
                                                                                      tttGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGc            <  2:237473/100‑1 (MQ=255)
                                                                                       ttGCGCCATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGcc           <  2:2421394/100‑1 (MQ=255)
                                                                                       ttGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGcc           <  2:4141554/100‑1 (MQ=255)
                                                                                       ttGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGcc           <  1:3885575/100‑1 (MQ=255)
                                                                                       ttGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGcc           >  2:1381995/1‑100 (MQ=255)
                                                                                        tGCGCAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGccc          >  2:390266/1‑100 (MQ=255)
                                                                                          cgcAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGa        >  2:3074055/1‑100 (MQ=255)
                                                                                           gcAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGac       >  2:2693002/1‑100 (MQ=255)
                                                                                           gcAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGac       >  1:4047909/1‑100 (MQ=255)
                                                                                            cAATCAGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGcc           >  1:4022884/1‑95 (MQ=255)
                                                                                               tCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCAAGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGacgtac   <  2:3863755/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
GTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  NC_000913/1193192‑1193387

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: