Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,171,570:1 +T 29.2% intergenic (+26/‑67) 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,5691.T29.2% 219.8 / 209.3 154intergenic (+25/‑68)rrlB/rrfB23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base . (27/82);  new base T (2/43);  total (29/125)
Fisher's exact test for biased strand distribution p-value = 2.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 155 total aligned reads displayed.

CGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCT  >  NC_000913/4171503‑4171668
                                                                    |                                                                                                   
cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAA                                                                    >  2:3699533‑M1/1‑100 (MQ=255)
cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAA                                                                    <  1:1673055‑M1/100‑1 (MQ=255)
cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAACTGTTTTGGCGG‑A‑TGAAAGAAGATTTTCAGCCTGATACAGATTAAA                                                                    >  2:425019‑M1/1‑100 (MQ=255)
 gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAAT                                                                   >  2:2363510‑M1/1‑100 (MQ=255)
 gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAAT                                                                   >  2:2135617‑M1/1‑100 (MQ=255)
 gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAAT                                                                   >  1:1484505‑M1/1‑100 (MQ=255)
  ttGAGCTAACCGGTAGTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                                  <  1:3887084‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGTTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                                  <  1:3827557‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                                  <  1:846854‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                                  <  1:1028382‑M1/100‑1 (MQ=255)
  ttGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATC                                                                  <  1:227036‑M1/100‑1 (MQ=255)
    gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAATCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGACACAGATTAAATCAG                                                                >  2:2301235‑M1/1‑100 (MQ=255)
    gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAG                                                                >  2:3637157‑M1/1‑100 (MQ=255)
    gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAG                                                                >  1:1496126‑M1/1‑100 (MQ=255)
    gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAG                                                                <  2:1922328‑M1/100‑1 (MQ=255)
      gCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAc                                                              <  1:2874891‑M1/100‑2 (MQ=255)
                 cTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgt                                                   <  1:703062‑M1/100‑13 (MQ=255)
                  tAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgtt                                                  >  2:3184019‑M1/1‑87 (MQ=255)
                   aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgtt                                                  >  1:2787051‑M1/1‑87 (MQ=255)
                    aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttc                                                 >  1:2604649‑M1/1‑86 (MQ=255)
                    aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttct                                                <  1:3483814‑M1/100‑16 (MQ=255)
                      gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctga                                              <  1:1729486‑M1/100‑18 (MQ=255)
                      gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctga                                              <  2:415235‑M1/100‑18 (MQ=255)
                      gAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctga                                              <  1:3718273‑M1/100‑18 (MQ=255)
                      gAACCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctga                                              <  1:2583609‑M1/100‑18 (MQ=255)
                       aaCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                             <  2:2807966‑M1/100‑19 (MQ=255)
                       aaCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                             <  2:496899‑M1/100‑19 (MQ=255)
                       aaCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                             <  2:3609277‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                             <  2:774936‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                             <  1:3464813‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTGTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                             <  1:2736571‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGTT‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                             <  2:2764023‑M1/100‑19 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  1:3295208‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  1:322734‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  1:711567‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  2:3047966‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  2:3066462‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  1:568618‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  2:1962096‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATAAAATCAGAccgcttctgcgttctgatt                                            <  2:1823511‑M1/100‑20 (MQ=255)
                        aCCGTGAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  1:2536083‑M1/100‑20 (MQ=255)
                         ccGTGAGTCTTAACCTTACAACGCCTAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  1:1348272‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCTAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  2:3496351‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  2:141903‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  2:2604632‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  2:2536937‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGTGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  1:313367‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatt                                            <  1:146010‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGACTAAATCAGAccgcttctgcgttctgatt                                            <  2:3318596‑M1/100‑20 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  2:3254676‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  1:1036832‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  2:2305577‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  2:885925‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  2:65588‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  1:2326573‑M1/100‑21 (MQ=255)
                         ccGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  1:3322905‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  2:570371‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  2:3538682‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  2:876550‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  1:408824‑M1/100‑21 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgat                                             <  2:3726969‑M1/98‑19 (MQ=255)
                          cGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTGTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgattt                                           <  2:619512‑M1/99‑21 (MQ=255)
                           gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGTGG‑A‑TGAGGGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaa                                         >  2:1313704‑M1/1‑78 (MQ=255)
                           gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaa                                         <  2:1951627‑M1/100‑23 (MQ=255)
                           gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaa                                         <  1:2998316‑M1/100‑23 (MQ=255)
                           gTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaa                                         <  2:3878473‑M1/100‑23 (MQ=255)
                            tGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaat                                        <  1:427081‑M1/100‑24 (MQ=255)
                            tGAGGCTTAACCTTACAACGCCGAAGCTGTGTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaat                                        <  2:2367557‑M1/100‑24 (MQ=255)
                             gAGGCTTAACCTTACAACGCCTAATCTGTTTTGGCGG‑ATTGAGCGAAGACTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaat                                        <  1:1189916‑M1/100‑24 (MQ=255)
                             gAGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                                       <  1:1969659‑M1/100‑25 (MQ=255)
                              agtcttaACCTTACAACGCCGAAGCTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                                       <  2:2723294‑M1/97‑25 (MQ=255)
                              aGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                                       <  2:910492‑M1/100‑25 (MQ=255)
                              aGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatc                                       <  2:553020‑M1/100‑25 (MQ=255)
                               ggCTTAACCTTACAACGCCGAATCTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatct                                      <  2:3672629‑M1/100‑26 (MQ=255)
                               ggCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctg                                     <  1:2210014‑M1/100‑27 (MQ=255)
                               ggCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctg                                     <  1:2085429‑M1/100‑27 (MQ=255)
                                gCTTAACCTTACAACGCCGAAGGTGGTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctg                                     <  2:1938441‑M1/100‑27 (MQ=255)
                                gCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctg                                     <  2:315343‑M1/100‑27 (MQ=255)
                                        ttACAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctgtatcaggc                             <  1:2842541‑M1/100‑35 (MQ=255)
                                          aCAACGCCGAAGGTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctgtatcaggctg                           <  1:527868‑M1/100‑37 (MQ=255)
                                            aaCGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctgtatcaggctgaaa                        <  1:3045758‑M1/100‑40 (MQ=255)
                                              cGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGGTTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctgtatcaggctgaaaat                      <  1:3798530‑M1/100‑42 (MQ=255)
                                              cGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCa                      >  2:2949471/1‑100 (MQ=11)
                                              cGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGAGGTCTGATAAAACAGAATTTGCCTGGCGGCa                      >  1:2224438/1‑100 (MQ=11)
                                               gCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCGGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                     <  2:1985665/100‑1 (MQ=14)
                                               gCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                     >  1:2649922/1‑100 (MQ=31)
                                               gCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAg                     <  1:820501/100‑1 (MQ=31)
                                                 cGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGCTACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa                   >  2:1086027/1‑100 (MQ=17)
                                                 cGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTa                   >  2:2994643/1‑100 (MQ=32)
                                                    agatgttTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctgtatcaggctgaaaatcttctc                <  2:3920378‑M1/97‑48 (MQ=255)
                                                     gCTGTTTTGGCGG‑ATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAccgcttctgcgttctgatttaatctgtatcaggctgaaaatcttctc                <  1:3230090‑M1/100‑48 (MQ=255)
                                                     gCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAgcgc               <  1:3675267/100‑1 (MQ=33)
                                                        gTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGGGGCAGTAgcgc               >  2:956241/1‑97 (MQ=17)
                                                        gTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggt            >  2:2025549/1‑100 (MQ=33)
                                                         ttttGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtg           >  1:2523326/1‑100 (MQ=33)
                                                         ttttGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtg           <  1:1252495/100‑1 (MQ=33)
                                                            tGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc        <  1:3185668/100‑1 (MQ=33)
                                                            tGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTc        <  1:164754/100‑1 (MQ=33)
                                                             ggCGG‑A‑TGAGAGAAGATTTACAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTcc       >  2:1784714/1‑100 (MQ=17)
                                                               cGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCa     >  1:3068934/1‑100 (MQ=33)
                                                                    |                                                                                                   
CGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGG‑A‑TGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCT  >  NC_000913/4171503‑4171668

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