Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,333,582 (A)8→9 intergenic (‑131/‑128) rplU ← / → ispB 50S ribosomal subunit protein L21/octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,333,5741.A95.3% 490.1 / 15.5 149intergenic (‑123/‑136)rplU/ispB50S ribosomal subunit protein L21/octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (2/5);  new base A (84/58);  total (86/63)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 169 total aligned reads displayed.

CCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTG  >  NC_000913/3333475‑3333668
                                                                                                      |                                                                                                   
ccTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag                                                                                                      >  1:712045/1‑100 (MQ=255)
ccTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag                                                                                                      >  1:1784000/1‑100 (MQ=255)
 cTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaag‑‑a                                                                                                   <  2:2636399/100‑1 (MQ=255)
 cTTTTCAATGACTCATAGTGCGCCATAAATACTCACAACACGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCCATACACgaag‑‑a                                                                                                   <  2:2495998/100‑1 (MQ=255)
      cAATGATTCAGCGAGCGCTATAAATATTCACAATATGGCGCGAATATTATCCAAAACGAACACCTTTGACAAGTGCTACCGTCCGTAAACGAAG‑‑aaaa                                                                                                >  2:1317660/1‑98 (MQ=255)
      cAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaa‑‑‑a                                                                                           >  2:2736002/1‑100 (MQ=255)
       aaTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaa‑‑‑aa                                                                                         >  1:3007315/1‑100 (MQ=255)
        aTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaa‑‑‑aaa                                                                                        <  2:1756551/100‑1 (MQ=255)
        aTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaa‑‑‑aaa                                                                                        <  1:3436354/100‑1 (MQ=255)
          gATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAaa‑‑aaaaa‑‑‑aaaac                                                                                      <  2:218710/100‑3 (MQ=255)
          gATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑aaaaa‑‑‑aaaac                                                                                      <  1:1902208/100‑3 (MQ=255)
          gATTCAGAGTGAGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGGGCTACAGTCAATACACGAAG‑‑aaaaa‑‑‑aaaac                                                                                      <  1:2280596/100‑3 (MQ=255)
            ttCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTAAATACACGAAG‑‑aaaaa‑‑‑aaaacac                                                                                    <  1:1499260/100‑5 (MQ=255)
              cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACaa                                                                                   >  2:260895/1‑100 (MQ=255)
              cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACaa                                                                                   >  1:323226/1‑100 (MQ=255)
              cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACaa                                                                                   >  1:785228/1‑100 (MQ=255)
                gagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAA‑‑‑AAAaacacaac                                                                                >  2:3152461/1‑96 (MQ=255)
                  gTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACt                                                                                 >  1:3431225/1‑98 (MQ=255)
                  gTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTg                                                                               >  2:2401965/1‑100 (MQ=255)
                    gggcTATAAATATGTACAGTAGGGGGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTg                                                                               <  2:1178077/96‑1 (MQ=255)
                    gcgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTa                                                                             >  2:2745534/1‑100 (MQ=255)
                    gcgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTa                                                                             >  1:1620754/1‑100 (MQ=255)
                    gcgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTa                                                                             <  2:1783987/100‑1 (MQ=255)
                      gcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACg                                                                           >  2:3605952/1‑100 (MQ=255)
                      gcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACg                                                                           >  1:3218690/1‑100 (MQ=255)
                         ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTa                                                                        >  1:3381589/1‑100 (MQ=255)
                         ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAAAAGAAAAAAAAAA‑‑‑AAACACAACTTGTACgg                                                                          <  1:841347/99‑1 (MQ=255)
                          taAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTAACAAGTGCTACAGTCAATACACGAAAAAAAAAA‑‑‑AAACACAACTTGTACGGTa                                                                        <  2:2839348/100‑1 (MQ=255)
                           aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAA‑‑‑AAACACAACTTGTACGGTaa                                                                       >  2:2315835/1‑100 (MQ=255)
                           aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAAc                                                                      >  2:1561799/1‑100 (MQ=255)
                           aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACCTGTACGGGaa                                                                       >  2:1871695/1‑99 (MQ=255)
                             atatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACACTCGATACACTAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGt                                                                    >  2:2929773/1‑100 (MQ=255)
                             aaatgcacaAAAGGGCGCGAGTAGTACGCAAAACGCACGCCTTTGATGAGTGCTACAGTCAATAAAAAAAAA‑‑AAAAA‑‑‑AAACACAACTTGTACGGTAACGt                                                                    <  2:3562580/95‑1 (MQ=255)
                              tatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGtt                                                                   <  1:1149946/100‑1 (MQ=255)
                                ttCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTa                                                                 >  2:2234495/1‑100 (MQ=255)
                                ttCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTa                                                                 >  1:2214696/1‑100 (MQ=255)
                                 tCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTa                                                                 >  2:1045709/1‑99 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGTTAACGTTTAt                                                                >  1:1703547/1‑99 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATc                                                               >  1:409752/1‑100 (MQ=255)
                                   acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATc                                                               <  2:897013/99‑1 (MQ=255)
                                    caATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACTGAAACGTTTATc                                                               >  2:1236672/1‑98 (MQ=255)
                                     aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCtgt                                                            <  1:3645852/100‑1 (MQ=255)
                                     aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACCTTTATc                                                               >  1:2620455/1‑97 (MQ=255)
                                     aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACAAAAAAAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCtgt                                                            <  2:1620742/100‑1 (MQ=255)
                                         gggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCa                                                        >  1:2671611/1‑100 (MQ=255)
                                         gggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCa                                                        >  1:2619215/1‑100 (MQ=255)
                                          ttccgatctatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                                       <  2:4128848/92‑1 (MQ=255)
                                          ggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                                       >  1:3753653/1‑100 (MQ=255)
                                           gcgcGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCAtt                                                      >  2:2699808/1‑100 (MQ=255)
                                            cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCAtt                                                      >  2:2032792/1‑99 (MQ=255)
                                              cgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCAttttt                                                   <  1:411309/100‑1 (MQ=255)
                                                aaTATTACGCAAAACGCACGCCTTTGACAAGGGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  1:3845914/100‑1 (MQ=255)
                                                aaTATTACGCAAAACGAACGCATTTGACAAGTGCTACAGTCAATAAAAAAAAA‑‑AAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  2:1492453/100‑1 (MQ=255)
                                                aaTATTACGCAAAACGAACACCTTTGACAAGTGCTACAGTCAATACATGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  2:215895/100‑1 (MQ=255)
                                                 atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                <  1:1988958/100‑1 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCACTACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               <  1:1881760/100‑1 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               >  2:3232377/1‑100 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               <  2:4080432/100‑1 (MQ=255)
                                                   atTACGCAAAACGAACGCCTTTAACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                <  2:3027672/98‑1 (MQ=255)
                                                    ttACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTa                                             <  1:3152123/100‑1 (MQ=255)
                                                       cGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAAACTAAACACAACTTCTACGGTAACGTTTTTCTGCGCCATTTTTTCAg                                               >  2:2383933/1‑97 (MQ=255)
                                                            aaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAcc                                     <  2:50170/100‑1 (MQ=255)
                                                             acgcacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAccc                                    <  1:2906428/100‑1 (MQ=255)
                                                              cgcacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCt                                   >  2:2130538/1‑100 (MQ=255)
                                                                cacgcCTTTGAGAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGCGCCATTTTTTTAGTATACTCCCCatat                                 >  2:1055071/1‑100 (MQ=255)
                                                                cacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtat                                 >  1:1910617/1‑100 (MQ=255)
                                                                 acgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtata                                <  1:2107816/100‑1 (MQ=255)
                                                                  cgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTatat                               <  1:1212448/100‑1 (MQ=255)
                                                                   gcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGATTATCTGTGCCATTTTTTCAGTACAATCACCCTATAtt                              >  2:312906/1‑100 (MQ=255)
                                                                    ccTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATAtt                              <  1:333170/100‑1 (MQ=255)
                                                                    ccTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTc                             >  1:632257/1‑100 (MQ=255)
                                                                     cTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATTACCCTATATTcc                            >  1:2034143/1‑100 (MQ=255)
                                                                       ttGACAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTa                          <  1:1416161/100‑1 (MQ=255)
                                                                          aCAAGTGCTTCATTCAATAAAAAAAAA‑‑AAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc                        <  2:3627789/99‑1 (MQ=255)
                                                                          aCAAGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                            <  2:3857050/95‑1 (MQ=255)
                                                                            aaGTGCTACAGTCAATACACAAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt                           <  2:3381550/94‑1 (MQ=255)
                                                                             aGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATa                    <  1:132097/100‑1 (MQ=255)
                                                                             aGTGCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATa                    <  2:127783/100‑1 (MQ=255)
                                                                              gggCTACAGGCAAAAAAAGAAAAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc                        <  2:2656692/93‑1 (MQ=255)
                                                                                gCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTTCCATTTTTTCAATACAATCACCCTATATTTCTAACCATAAAc                 >  1:1288957/1‑100 (MQ=255)
                                                                                gCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAc                 <  2:291148/100‑1 (MQ=255)
                                                                                gCTACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTAGTACGGTAACGTTCATCTGTGCCATTTCTTCAGCACAAACACCCTATATTCCTAACCATaa                   >  1:2952822/1‑98 (MQ=255)
                                                                                gCTACAGTCAATACACAAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAc                 <  2:3824598/100‑1 (MQ=255)
                                                                                  tACAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAATACAATCACCTTATATTCCTAACCATaatctc               >  2:80672/1‑96 (MQ=255)
                                                                                   aCAGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCt              >  1:615780/1‑100 (MQ=255)
                                                                                    cggTCAATACACGAAGAAAAAA‑‑‑AAACACAACATGTACGGTAAGGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTa             >  1:2821268/3‑100 (MQ=255)
                                                                                     aGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTaa            >  2:628404/1‑100 (MQ=255)
                                                                                     aGTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTATGCCATTTTATCAGTTCAATCACCCTATATTCCTAACCATAAACCCTa             >  2:3696780/1‑99 (MQ=255)
                                                                                      gTCAATACAGGAAG‑‑AAAAA‑‑‑ACAATACACCTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATCTTCCTAACCATAAACCCt              <  2:3334003/97‑1 (MQ=255)
                                                                                      gTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAg           <  1:2721075/100‑1 (MQ=255)
                                                                                      gTCAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAg           >  1:911756/1‑100 (MQ=255)
                                                                                        cAATACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGtt         >  2:3908646/1‑100 (MQ=255)
                                                                                         aaTACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTTTTTACAATAAACCTATATTCCTCATCATAAACTCTAAttt         >  2:824742/1‑99 (MQ=255)
                                                                                           tACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGcc      >  1:935967/1‑100 (MQ=255)
                                                                                           tACACGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGcc      <  1:436006/100‑1 (MQ=255)
                                                                                                cGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCATTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  >  2:2225229/1‑100 (MQ=255)
                                                                                                cGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCtt    >  1:3906187/1‑98 (MQ=255)
                                                                                                cGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  <  1:1058590/100‑1 (MQ=255)
                                                                                                cGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  <  1:1045714/100‑1 (MQ=255)
                                                                                                cGAAGAAAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACACTATATTCCTAACCATAAACCCTAAggttccgt    >  2:343101/1‑91 (MQ=255)
                                                                                                      |                                                                                                   
CCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAG‑‑AAAAA‑‑‑AAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTG  >  NC_000913/3333475‑3333668

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