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.

CCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTG  >  NC_000913/3333475‑3333668
                                                                                                      |                                                                                                   
ccTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACg‑aag                                                                                                      >  1:712045/1‑100 (MQ=255)
ccTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACg‑aag                                                                                                      >  1:1784000/1‑100 (MQ=255)
 cTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACg‑aag‑‑a                                                                                                   <  2:2636399/100‑1 (MQ=255)
 cTTTTCAATGACTCATAGTGCGCCATAAATACTCACAACACGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCCATAC‑ACg‑aag‑‑a                                                                                                   <  2:2495998/100‑1 (MQ=255)
      cAATGATTCAGCGAGCGCTATAAATATTCACAATATGGCGCGAATATTATCCAAAACGAACACCTTTGACAAGTGCTACCGTCCGTAA‑ACG‑AAG‑‑aaaa                                                                                                >  2:1317660/1‑98 (MQ=255)
      cAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑‑aaaaa‑‑‑a                                                                                           >  2:2736002/1‑100 (MQ=255)
       aaTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑‑aaaaa‑‑‑a‑a                                                                                         >  1:3007315/1‑100 (MQ=255)
        aTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑‑aaaaa‑‑‑a‑aa                                                                                        <  2:1756551/100‑1 (MQ=255)
        aTGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑‑aaaaa‑‑‑a‑aa                                                                                        <  1:3436354/100‑1 (MQ=255)
          gATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑Aaa‑‑aaaaa‑‑‑a‑aaac                                                                                      <  2:218710/100‑3 (MQ=255)
          gATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑‑aaaaa‑‑‑a‑aaac                                                                                      <  1:1902208/100‑3 (MQ=255)
          gATTCAGAGTGAGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGGGCTACAGTCAATAC‑ACG‑AAG‑‑aaaaa‑‑‑a‑aaac                                                                                      <  1:2280596/100‑3 (MQ=255)
            ttCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTAAATAC‑ACG‑AAG‑‑aaaaa‑‑‑a‑aaacac                                                                                    <  1:1499260/100‑5 (MQ=255)
              cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACaa                                                                                   >  2:260895/1‑100 (MQ=255)
              cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACaa                                                                                   >  1:323226/1‑100 (MQ=255)
              cAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACaa                                                                                   >  1:785228/1‑100 (MQ=255)
                gagTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑‑AAAAA‑‑‑A‑AAaacacaac                                                                                >  2:3152461/1‑96 (MQ=255)
                  gTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACt                                                                                 >  1:3431225/1‑98 (MQ=255)
                  gTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTg                                                                               >  2:2401965/1‑100 (MQ=255)
                    gggcTATAAATATGTACAGTAGGGGGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTg                                                                               <  2:1178077/96‑1 (MQ=255)
                    gcgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTa                                                                             >  2:2745534/1‑100 (MQ=255)
                    gcgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTa                                                                             >  1:1620754/1‑100 (MQ=255)
                    gcgcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTa                                                                             <  2:1783987/100‑1 (MQ=255)
                      gcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACg                                                                           >  2:3605952/1‑100 (MQ=255)
                      gcTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACg                                                                           >  1:3218690/1‑100 (MQ=255)
                         ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTa                                                                        >  1:3381589/1‑100 (MQ=255)
                         ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAA‑AAG‑AAAAAAAAAA‑‑‑A‑AACACAACTTGTACgg                                                                          <  1:841347/99‑1 (MQ=255)
                          taAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTAACAAGTGCTACAGTCAATAC‑ACG‑AAAAAAAAAA‑‑‑A‑AACACAACTTGTACGGTa                                                                        <  2:2839348/100‑1 (MQ=255)
                           aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAA‑‑‑A‑AACACAACTTGTACGGTaa                                                                       >  2:2315835/1‑100 (MQ=255)
                           aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAAc                                                                      >  2:1561799/1‑100 (MQ=255)
                           aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACCTGTACGGGaa                                                                       >  2:1871695/1‑99 (MQ=255)
                             atatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACACTCGATAC‑ACT‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGt                                                                    >  2:2929773/1‑100 (MQ=255)
                             aaatgcacaAAAGGGCGCGAGTAGTACGCAAAACGCACGCCTTTGATGAGTGCTACAGTCAATAAAAAA‑AAA‑‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGt                                                                    <  2:3562580/95‑1 (MQ=255)
                              tatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGtt                                                                   <  1:1149946/100‑1 (MQ=255)
                                ttCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTa                                                                 >  2:2234495/1‑100 (MQ=255)
                                ttCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTa                                                                 >  1:2214696/1‑100 (MQ=255)
                                 tCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTa                                                                 >  2:1045709/1‑99 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGTTAACGTTTAt                                                                >  1:1703547/1‑99 (MQ=255)
                                  cacaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATc                                                               >  1:409752/1‑100 (MQ=255)
                                   acaATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATc                                                               <  2:897013/99‑1 (MQ=255)
                                    caATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACTGAAACGTTTATc                                                               >  2:1236672/1‑98 (MQ=255)
                                     aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCtgt                                                            <  1:3645852/100‑1 (MQ=255)
                                     aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACCTTTATc                                                               >  1:2620455/1‑97 (MQ=255)
                                     aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑AAA‑AAAA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCtgt                                                            <  2:1620742/100‑1 (MQ=255)
                                         gggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCa                                                        >  1:2671611/1‑100 (MQ=255)
                                         gggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCa                                                        >  1:2619215/1‑100 (MQ=255)
                                          ttccgatctatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                                       <  2:4128848/92‑1 (MQ=255)
                                          ggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                                       >  1:3753653/1‑100 (MQ=255)
                                           gcgcGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCAtt                                                      >  2:2699808/1‑100 (MQ=255)
                                            cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCAtt                                                      >  2:2032792/1‑99 (MQ=255)
                                              cgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCAttttt                                                   <  1:411309/100‑1 (MQ=255)
                                                aaTATTACGCAAAACGCACGCCTTTGACAAGGGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  1:3845914/100‑1 (MQ=255)
                                                aaTATTACGCAAAACGAACGCATTTGACAAGTGCTACAGTCAATAA‑AAAAAAA‑‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  2:1492453/100‑1 (MQ=255)
                                                aaTATTACGCAAAACGAACACCTTTGACAAGTGCTACAGTCAATAC‑ATG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTc                                                 <  2:215895/100‑1 (MQ=255)
                                                 atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                <  1:1988958/100‑1 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCACTAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               <  1:1881760/100‑1 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               >  2:3232377/1‑100 (MQ=255)
                                                  tatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAg                                               <  2:4080432/100‑1 (MQ=255)
                                                   atTACGCAAAACGAACGCCTTTAACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCa                                                <  2:3027672/98‑1 (MQ=255)
                                                    ttACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTa                                             <  1:3152123/100‑1 (MQ=255)
                                                       cGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑‑AAAAAACTA‑AACACAACTTCTACGGTAACGTTTTTCTGCGCCATTTTTTCAg                                               >  2:2383933/1‑97 (MQ=255)
                                                            aaCGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAcc                                     <  2:50170/100‑1 (MQ=255)
                                                             acgcacgcCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAccc                                    <  1:2906428/100‑1 (MQ=255)
                                                              cgcacgcCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCt                                   >  2:2130538/1‑100 (MQ=255)
                                                                cacgcCTTTGAGAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGCGCCATTTTTTTAGTATACTCCCCatat                                 >  2:1055071/1‑100 (MQ=255)
                                                                cacgcCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtat                                 >  1:1910617/1‑100 (MQ=255)
                                                                 acgcCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtata                                <  1:2107816/100‑1 (MQ=255)
                                                                  cgcCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTatat                               <  1:1212448/100‑1 (MQ=255)
                                                                   gcCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGATTATCTGTGCCATTTTTTCAGTACAATCACCCTATAtt                              >  2:312906/1‑100 (MQ=255)
                                                                    ccTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGAAAAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATAtt                              <  1:333170/100‑1 (MQ=255)
                                                                    ccTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTc                             >  1:632257/1‑100 (MQ=255)
                                                                     cTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATTACCCTATATTcc                            >  1:2034143/1‑100 (MQ=255)
                                                                       ttGACAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTa                          <  1:1416161/100‑1 (MQ=255)
                                                                          aCAAGTGCTTCATTCAATAA‑AAAAAAA‑‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc                        <  2:3627789/99‑1 (MQ=255)
                                                                          aCAAGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTcc                            <  2:3857050/95‑1 (MQ=255)
                                                                            aaGTGCTACAGTCAATAC‑ACA‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt                           <  2:3381550/94‑1 (MQ=255)
                                                                             aGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATa                    <  1:132097/100‑1 (MQ=255)
                                                                             aGTGCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATa                    <  2:127783/100‑1 (MQ=255)
                                                                              gggCTACAGGCAAAAA‑AAG‑AAAA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAAc                        <  2:2656692/93‑1 (MQ=255)
                                                                                gCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTTCCATTTTTTCAATACAATCACCCTATATTTCTAACCATAAAc                 >  1:1288957/1‑100 (MQ=255)
                                                                                gCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAc                 <  2:291148/100‑1 (MQ=255)
                                                                                gCTACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTAGTACGGTAACGTTCATCTGTGCCATTTCTTCAGCACAAACACCCTATATTCCTAACCATaa                   >  1:2952822/1‑98 (MQ=255)
                                                                                gCTACAGTCAATAC‑ACA‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAAc                 <  2:3824598/100‑1 (MQ=255)
                                                                                  tACAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAATACAATCACCTTATATTCCTAACCATaatctc               >  2:80672/1‑96 (MQ=255)
                                                                                   aCAGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCt              >  1:615780/1‑100 (MQ=255)
                                                                                    cggTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACATGTACGGTAAGGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTa             >  1:2821268/3‑100 (MQ=255)
                                                                                     aGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTaa            >  2:628404/1‑100 (MQ=255)
                                                                                     aGTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTATGCCATTTTATCAGTTCAATCACCCTATATTCCTAACCATAAACCCTa             >  2:3696780/1‑99 (MQ=255)
                                                                                      gTCAATAC‑AGG‑AAG‑‑AAAAA‑‑‑ACAATACACCTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATCTTCCTAACCATAAACCCt              <  2:3334003/97‑1 (MQ=255)
                                                                                      gTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAg           <  1:2721075/100‑1 (MQ=255)
                                                                                      gTCAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAg           >  1:911756/1‑100 (MQ=255)
                                                                                        cAATAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGtt         >  2:3908646/1‑100 (MQ=255)
                                                                                         aaTAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTTTTTACAATAAACCTATATTCCTCATCATAAACTCTAAttt         >  2:824742/1‑99 (MQ=255)
                                                                                           tAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGcc      >  1:935967/1‑100 (MQ=255)
                                                                                           tAC‑ACG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGcc      <  1:436006/100‑1 (MQ=255)
                                                                                                cG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCATTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  >  2:2225229/1‑100 (MQ=255)
                                                                                                cG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCtt    >  1:3906187/1‑98 (MQ=255)
                                                                                                cG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  <  1:1058590/100‑1 (MQ=255)
                                                                                                cG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTg  <  1:1045714/100‑1 (MQ=255)
                                                                                                cG‑AAGA‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACACTATATTCCTAACCATAAACCCTAAggttccgt    >  2:343101/1‑91 (MQ=255)
                                                                                                      |                                                                                                   
CCTTTTCAATGATTCAGAGTGCGCTATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATAC‑ACG‑AAG‑‑AAAAA‑‑‑A‑AACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTG  >  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: ‑