Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 668,691 C→A G9* (GGA→TGA)  rlmH ← 23S rRNA m(3)Psi1915 pseudouridine methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913668,6910CA98.3% 227.0 / ‑2.8 60G9* (GGA→TGA) rlmH23S rRNA m(3)Psi1915 pseudouridine methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (28/31);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGCCG  >  NC_000913/668544‑668836
                                                                                                                                                   |                                                                                                                                                 
cACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCacac                                                                                                                                                <  1:239550/151‑1 (MQ=255)
   ccTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGc                                                                                                                                             <  2:952669/151‑1 (MQ=255)
     tctTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCg                                                                                                                                            >  1:520419/1‑150 (MQ=255)
     tctTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCCCGAAGGGCTTATATTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGCTTGTACCCAGTCCGGCATTTTAGTTCACATgg                                                                                                                                              >  2:458355/1‑148 (MQ=255)
      ctTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAGAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGAc                                                                                                                                          >  2:1211012/1‑151 (MQ=255)
        ttGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGa                                                                                                                                           <  2:954127/148‑1 (MQ=255)
        ttGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACa                                                                                                                                         >  1:386183/1‑150 (MQ=255)
          gTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAg                                                                                                                                       <  2:75422/150‑1 (MQ=255)
           tCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGtt                                                                                                                                     >  1:1107705/1‑151 (MQ=255)
                 aTGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCt                                                                                                                               <  1:968803/151‑1 (MQ=255)
                            gTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTa                                                                                                                     >  1:1212415/1‑150 (MQ=255)
                            gTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTa                                                                                                                     <  2:462698/150‑1 (MQ=255)
                              ccGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACt                                                                                                                  <  2:756151/151‑1 (MQ=255)
                                    tCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCag                                                                                                              >  2:370330/1‑149 (MQ=255)
                                       tGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCagag                                                                                                            <  2:53827/148‑1 (MQ=255)
                                        gCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGtttt                                                                                                        <  2:213938/151‑1 (MQ=255)
                                        gCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGtttt                                                                                                        >  2:370359/1‑151 (MQ=255)
                                         ccGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGttttt                                                                                                       >  1:518089/1‑151 (MQ=255)
                                         ccGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCAACAAGTTGCAGCTTCACGCATTAACTCCAGAGttttt                                                                                                       >  2:733944/1‑151 (MQ=255)
                                          cgcgTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTc                                                                                                      >  2:488552/1‑151 (MQ=255)
                                           gcgTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTcc                                                                                                     <  2:1107815/151‑1 (MQ=255)
                                           gcgTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTcc                                                                                                     <  2:1124478/151‑1 (MQ=255)
                                             gTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTcc                                                                                                     >  2:365548/1‑149 (MQ=255)
                                             gTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAg                                                                                                   <  2:791253/151‑1 (MQ=255)
                                             gTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAg                                                                                                   <  2:796619/151‑1 (MQ=255)
                                                   cggccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATa                                                                                             <  1:661306/151‑1 (MQ=255)
                                                    ggccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATAc                                                                                            >  2:419732/1‑151 (MQ=255)
                                                      ccggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATAc                                                                                            >  2:291331/1‑149 (MQ=255)
                                                       cggAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAgg                                                                                         >  2:1066520/1‑151 (MQ=255)
                                                         gAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCg                                                                                       >  1:1023202/1‑151 (MQ=255)
                                                          aaTTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCg                                                                                       >  1:13726/1‑150 (MQ=255)
                                                          aaTTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGa                                                                                      <  1:617852/151‑1 (MQ=255)
                                                           aTTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGAc                                                                                     >  2:12233/1‑151 (MQ=255)
                                                             ttCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACg                                                                                    <  2:693787/150‑1 (MQ=255)
                                                                aaTCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGctc                                                                                >  1:879865/1‑151 (MQ=255)
                                                                            gggCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATg                                                                    <  1:953848/151‑1 (MQ=255)
                                                                              gCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGAc                                                                  <  2:399666/151‑1 (MQ=255)
                                                                               cATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGAc                                                                  >  2:221077/1‑150 (MQ=255)
                                                                                  atCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATgg                                                              >  1:383510/1‑151 (MQ=255)
                                                                                      ttCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACa                                                           >  1:1231644/1‑150 (MQ=255)
                                                                                       tCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAAt                                                         >  1:609315/1‑151 (MQ=255)
                                                                                              acgacgCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCg                                                  <  1:877272/151‑1 (MQ=255)
                                                                                                gacgCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGcc                                                >  2:560911/1‑151 (MQ=255)
                                                                                                    cAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCaaa                                            <  2:581913/151‑1 (MQ=255)
                                                                                                      ggTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCaaa                                            <  2:402888/149‑1 (MQ=255)
                                                                                                         aCTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCAc                                       <  1:613260/151‑1 (MQ=255)
                                                                                                                 aaaCCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATcc                               <  2:945593/151‑1 (MQ=255)
                                                                                                                 aaaCCGGTTTATACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCATCTTCACGCATTAACTCCAGAGTTTTACCAGTTCATACAGGAGACGGCTCTCTTCCTGCATGACATAGACAATCACATCGCCCAAATCCACGACAATcc                               <  1:458308/151‑1 (MQ=255)
                                                                                                                   aCCGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAg                             <  1:370325/151‑1 (MQ=255)
                                                                                                                    ccGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGt                            <  2:423318/151‑1 (MQ=255)
                                                                                                                     cGGTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTc                           <  2:386220/151‑1 (MQ=255)
                                                                                                                      ggTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTc                           >  1:734400/1‑150 (MQ=255)
                                                                                                                       gTTTGTACCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCgg                         <  1:1226328/151‑1 (MQ=255)
                                                                                                                             aCCCAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGt                   <  1:1021109/151‑1 (MQ=255)
                                                                                                                              cccAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGtt                  <  1:1188847/151‑1 (MQ=255)
                                                                                                                               ccAGTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGttt                 >  2:675231/1‑151 (MQ=255)
                                                                                                                                  gTCCGGCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTc               >  2:1181873/1‑150 (MQ=255)
                                                                                                                                      ggCATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCtt          >  2:704872/1‑151 (MQ=255)
                                                                                                                                        cATTTTCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCt        <  1:1132892/151‑1 (MQ=255)
                                                                                                                                             tCGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGc    <  1:1210888/150‑1 (MQ=255)
                                                                                                                                              cGTTCACACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGCCg  <  1:874560/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACCCAGTCCGGCATTTTCGTTCCCACGGCGACAAGTTGCAGCTTCACGCATTAACTCCAGAGTTTTTCCAGTTCATACAGGCGACGGCTCTCTTCCTGCATGACATGGACAATCACATCGCCCAAATCCACGACAATCCAGTCGGCGCTGTTTTCACCTTCTACGCCG  >  NC_000913/668544‑668836

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: