Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,171,557 C→G 23.7% intergenic (+13/‑80) rrlB → / → rrfB 23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,171,5570CG23.7% 252.4 / 89.4 142intergenic (+13/‑80)rrlB/rrfB23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base C (19/83);  new base G (2/27);  total (22/120)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-02

Only 100 of 142 total aligned reads displayed.

CGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGC  >  NC_000913/4171503‑4171655
                                                      |                                                                                                    
cGTTGAGCTAATCGGTACTAATTAACCCTGAGGCTTAACCTTACAACGTCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAA                                                       <  2:486253‑M1/100‑1 (MQ=255)
cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAA                                                       <  1:1673055‑M1/100‑1 (MQ=255)
cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAA                                                       <  1:1678990‑M1/100‑1 (MQ=255)
cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAA                                                       >  2:3699533‑M1/1‑100 (MQ=255)
cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAACTGTTTTGGCGGATGAAAGAAGATTTTCAGCCTGATACAGATTAAA                                                       >  2:425019‑M1/1‑100 (MQ=255)
 gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAAT                                                      >  1:3233641‑M1/1‑100 (MQ=255)
 gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAAT                                                      >  1:1484505‑M1/1‑100 (MQ=255)
 gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAAT                                                      >  2:2363510‑M1/1‑100 (MQ=255)
 gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAAT                                                      >  2:2135617‑M1/1‑100 (MQ=255)
 gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCGGATACAGAATAAAT                                                      >  2:937003‑M1/1‑100 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGTTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                     <  1:3827557‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                     <  1:1028382‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                     <  2:3108395‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                     <  1:846854‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                     <  1:227036‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                     <  2:536375‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                     <  1:1165692‑M1/100‑1 (MQ=255)
   tGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCA                                                    <  2:3630693‑M1/100‑1 (MQ=255)
    gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAATCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGACACAGATTAAATCAG                                                   >  2:2301235‑M1/1‑100 (MQ=255)
    gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAG                                                   <  1:1945634‑M1/100‑1 (MQ=255)
    gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAG                                                   <  2:1922328‑M1/100‑1 (MQ=255)
    gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAG                                                   >  1:1496126‑M1/1‑100 (MQ=255)
    gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAG                                                   >  2:3637157‑M1/1‑100 (MQ=255)
      gCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAc                                                 <  1:2874891‑M1/100‑2 (MQ=255)
                 cTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgt                                      <  1:703062‑M1/100‑13 (MQ=255)
                  tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgtt                                     >  2:3184019‑M1/1‑87 (MQ=255)
                   aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgtt                                     >  1:2787051‑M1/1‑87 (MQ=255)
                    aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttc                                    >  1:2604649‑M1/1‑86 (MQ=255)
                    aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttc                                    <  2:2787034‑M1/100‑15 (MQ=255)
                    aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttct                                   <  1:3483814‑M1/100‑16 (MQ=255)
                      gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctga                                 <  1:3718273‑M1/100‑18 (MQ=255)
                      gAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctga                                 <  1:2583609‑M1/100‑18 (MQ=255)
                       aaCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctga                                 <  1:3733472‑M1/100‑18 (MQ=255)
                       aaCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                <  2:3609277‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                <  2:774936‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                <  1:3464813‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTGTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                <  1:2736571‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGTTTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                <  2:2764023‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  1:568618‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  2:1962096‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  1:711567‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  2:1986979‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  2:3047966‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  1:3295208‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  1:322734‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATAAAATCAGAccgcttctgcgttctgatt                               <  2:1823511‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  1:2536083‑M1/100‑20 (MQ=255)
                        aCCGAGAGCCTTAACCGTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGCTTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctga                                 <  2:1945543‑M1/98‑18 (MQ=255)
                         ccGTGAGTCTTAACCTTACAACGCCTAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  1:1348272‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCTAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:3496351‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  2:2604632‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  2:2536937‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  2:3531475‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  1:2537540‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGTGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                               <  1:313367‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGACTAAATCAGAccgcttctgcgttctgatt                               <  2:3318596‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  1:1036832‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:65588‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:2305577‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:3254676‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  1:3322905‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  1:2326573‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:369369‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTTATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:3933264‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:197403‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  1:3062365‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  1:380644‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:3034382‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  1:503548‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  1:2085443‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:876550‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  1:408824‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTGTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                              <  2:619512‑M1/99‑21 (MQ=255)
                           gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGTGGATGAGGGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaa                            >  2:1313704‑M1/1‑78 (MQ=255)
                           gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaa                            <  2:3878473‑M1/100‑23 (MQ=255)
                           gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaa                            <  1:2998316‑M1/100‑23 (MQ=255)
                           gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGAATAATACAGATTAAATCAGAccgcttctgcgttctgatttaa                            <  2:283804‑M1/100‑23 (MQ=255)
                             gAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                          <  1:3191567‑M1/100‑25 (MQ=255)
                             gAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                          <  1:1969659‑M1/100‑25 (MQ=255)
                              agtcttaACCTTACAACGCCGAAGCTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                          <  2:2723294‑M1/97‑25 (MQ=255)
                              aGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                          <  2:910492‑M1/100‑25 (MQ=255)
                              aGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                          <  1:1325308‑M1/100‑25 (MQ=255)
                              aGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                          <  2:553020‑M1/100‑25 (MQ=255)
                              aGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatct                         <  2:3786500‑M1/100‑26 (MQ=255)
                              aGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatct                         <  1:391258‑M1/100‑26 (MQ=255)
                              aGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatct                         <  1:1481101‑M1/100‑26 (MQ=255)
                               ggCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctg                        <  1:2210014‑M1/100‑27 (MQ=255)
                                gCTTAACCTTACAACGCCGAAGGTGGTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctg                        <  2:1938441‑M1/100‑27 (MQ=255)
                                gCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctg                        <  2:315343‑M1/100‑27 (MQ=255)
                                        ttACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctgtatcaggc                <  1:2842541‑M1/100‑35 (MQ=255)
                                          aCAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctgtatcaggctg              <  1:527868‑M1/100‑37 (MQ=255)
                                              cGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa         >  2:2949471/1‑100 (MQ=11)
                                              cGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGAGGTCTGATAAAACAGAATTTGCCTGGCGGCa         >  1:2224438/1‑100 (MQ=11)
                                               gCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg        <  1:820501/100‑1 (MQ=31)
                                               gCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg        >  1:2649922/1‑100 (MQ=31)
                                                 cGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGCTACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa      >  2:1086027/1‑100 (MQ=17)
                                                 cGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGCCTGATAAAACAGAATTTGCCTGGCGGAAGTa      >  2:583324/1‑100 (MQ=9)
                                                    aGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTTATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgcg   <  2:1656504/100‑1 (MQ=17)
                                                     gCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgcgc  <  1:3675267/100‑1 (MQ=33)
                                                     gCTGTTTTGGCGGATGAGAGAAGATTTTCAGCATGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgcgc  <  2:1310924/100‑1 (MQ=18)
                                                      |                                                                                                    
CGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGC  >  NC_000913/4171503‑4171655

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