Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,932,640 +C intergenic (‑107/+43) rluA ← / ← minE bifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,932,6381.C98.3% 201.4 / ‑4.8 61intergenic (‑105/+45)rluA/minEbifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (31/28);  minor base A (1/0);  total (33/28)
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.09e-01

GCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATT  >  NC_002947/1932496‑1932780
                                                                                                                                                 |                                                                                                                                               
gcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGGCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgcg                                                                                                                                           <  3:383588/149‑5 (MQ=255)
 cgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgcga                                                                                                                                          <  7:388811/149‑6 (MQ=255)
     ttCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCCTTGCAGCCCTTCGCGGGCAAGCCCGCACCTACACGGGAAGCATAACCATTGTAGGAGCGGGCATGACCGCGAAggg                                                                                                                                       >  1:95309/1‑149 (MQ=25)
        gAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTg                                                                                                                                   >  8:250809/1‑149 (MQ=255)
           gagTATCTGGACATTCGACAGCGGCTTTGCAGGTATTACTAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAATTATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCaa                                                                                                                                <  6:383315/149‑1 (MQ=37)
              gATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGc                                                                                                                             <  5:236890/149‑1 (MQ=255)
                    gACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGccc                                                                                                                        >  5:331/1‑148 (MQ=255)
                         tCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCCAACCCCCTCCTAAAAAGGGTATTCGTGCCCGGTGTGAGCGGGCGTGGCCGCGCAGGGCAGGCAAGGCGCCCCaaaat                                                                                                                  >  4:13007/1‑147 (MQ=255)
                         tCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGGTATGCGTCCCCGGTTGGAGCGGGGGTGCCCGCGAAGGGCGGGAAGGCAGCCCCACAtt                                                                                                                  >  4:415645/1‑149 (MQ=38)
                              cgCGGCATTGCAGGTCTCACGAAAATAGCAGTAACCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGa                                                                                                             <  3:368106/148‑1 (MQ=255)
                                cGGCATTGCCGGTCTCACGAAAACAGCAGTAAACGCGGGGCCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCGCCTACAAAGGATAGCCGCCCCCTGTAGGAGCGGGCTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                           <  3:95977/149‑1 (MQ=37)
                                  gCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAc                                                                                                         <  3:386235/149‑1 (MQ=255)
                                        cAGGCCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGAt                                                                                                   >  8:248871/1‑149 (MQ=255)
                                           gTCTCACTAAAATAGAAGTAAACGCTGGGGCCGCTTTGCGGCCGTTCGCGGGCAAGCCCGCTCCTACAAAGTATATGCGTCCCCTGTAGGAGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCa                                                                                                 <  4:171300/148‑1 (MQ=37)
                                           gTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAg                                                                                                >  5:394322/1‑149 (MQ=255)
                                              tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGa                                                                                             >  7:337468/1‑149 (MQ=255)
                                               cACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGAt                                                                                            >  8:46860/1‑149 (MQ=255)
                                                 cGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCAACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATcc                                                                                          <  5:216869/149‑1 (MQ=255)
                                                               accgcgGGGGCCGCTTGGGGGCCCTTCGCGGGCACCCCCGCTCAAACAAAGTATAGGCGCCCCCTGTTGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGtt                                                                            <  3:415645/147‑1 (MQ=37)
                                                               aaCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCGGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGGCCCCAATTGTCGACTACGCGGATCAGCGATGCGGCAGGGTGATGtt                                                                            >  1:254958/1‑149 (MQ=37)
                                                                  gcgggggCGGCTGTGGGCCCGTGCGCGGGAACCCCCCCTCATACAACGTATATGCGGCCCCTGTCGGAGGGGGCTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAg                                                                         <  8:174407/146‑1 (MQ=255)
                                                                      gggCCGGTTGGCGGCCCTGCGCGGGCACCCCCGCTCCTACAACGTATATGCGCCCCCTGTCGGCGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTcc                                                                     <  3:201458/149‑1 (MQ=37)
                                                                       ggCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGGTATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTc                                                                      >  2:347189/1‑147 (MQ=255)
                                                                       ggCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCa                                                                    <  7:46860/149‑1 (MQ=255)
                                                                        gCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCACCAACAAAGGATATGCGTCCCCAGTAAGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAATGTAGACAACGCCGATCAGCCATCCCGCAGGGTGATGTTGAGCTGCa                                                                    >  2:344437/1‑148 (MQ=25)
                                                                         ccGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAg                                                                   >  5:419850/1‑148 (MQ=255)
                                                                           gCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAAt                                                                >  3:135354/1‑149 (MQ=255)
                                                                               tGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa                                                             <  8:337468/148‑1 (MQ=255)
                                                                               ggCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTTTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa                                                             <  1:344437/147‑1 (MQ=255)
                                                                                gCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag                                                            <  4:135354/148‑1 (MQ=255)
                                                                                  ggCCGTTCGCGGGCAAGCCCGCTCATACAAAGTATATGCGTCCCGTTTAGTAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagcag                                                         <  6:331/149‑1 (MQ=38)
                                                                                       ttCGCGGGCAAGCCCGCTCCTAAAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGGAAGGCAGCCCCAAATTGTCGACTACGCGGAGAAGGGATGCGGGAGGGTGATGGTCAGCTGCAGAATAGGGCAGCTGcg                                                    >  4:304774/1‑148 (MQ=37)
                                                                                              gCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCa                                                                    >  7:78983/1‑126 (MQ=255)
                                                                                              gCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCa                                                                    <  8:78983/126‑1 (MQ=255)
                                                                                              gCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGttt                                             <  1:347189/149‑1 (MQ=255)
                                                                                               cAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtttt                                            >  3:86074/1‑149 (MQ=255)
                                                                                                   cccGCTCCTACAAAGGATATGCGTCCCCTGGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGTTCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                               >  1:272394/1‑127 (MQ=38)
                                                                                                   cccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGTTCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                               <  2:272394/126‑1 (MQ=255)
                                                                                                    ccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGc                                       >  1:276170/1‑149 (MQ=255)
                                                                                                          cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGat                                  <  7:248871/148‑1 (MQ=255)
                                                                                                          cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGat                                  <  6:394322/148‑1 (MQ=255)
                                                                                                            aCAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                               >  4:109437/1‑149 (MQ=255)
                                                                                                              aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc                                                                           <  4:102216/103‑1 (MQ=255)
                                                                                                              aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc                                                                           >  3:102216/1‑103 (MQ=255)
                                                                                                              aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGca                             >  2:118968/1‑149 (MQ=255)
                                                                                                                    tatGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                           <  6:367994/129‑1 (MQ=255)
                                                                                                                    tatGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                           >  5:367994/1‑129 (MQ=255)
                                                                                                                       gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatcat                       <  7:273521/146‑1 (MQ=255)
                                                                                                                       gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatcat                       >  8:273521/1‑146 (MQ=255)
                                                                                                                       gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGt                    >  1:352441/1‑149 (MQ=255)
                                                                                                                       gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGt                    >  3:245559/1‑149 (MQ=255)
                                                                                                                       gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGt                    >  8:12591/1‑149 (MQ=255)
                                                                                                                          tcctttGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg                  <  8:77993/143‑1 (MQ=255)
                                                                                                                            cccTGTCGGCGCGGGCGTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCg               <  3:13007/149‑1 (MQ=39)
                                                                                                                            cccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTACCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGcc                >  2:154809/1‑148 (MQ=255)
                                                                                                                              cTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGcac                            <  5:327430/134‑1 (MQ=255)
                                                                                                                              cTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGcac                            >  6:327430/1‑134 (MQ=255)
                                                                                                                                     gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGca                             <  3:420733/127‑1 (MQ=255)
                                                                                                                                     gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGca                             >  4:420733/1‑127 (MQ=255)
                                                                                                                                       gCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTcaca      >  7:356442/1‑148 (MQ=255)
                                                                                                                                        cGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca        >  6:77025/1‑145 (MQ=255)
                                                                                                                                        cGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca        <  5:77025/145‑1 (MQ=255)
                                                                                                                                          ggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAtt  >  8:220638/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATT  >  NC_002947/1932496‑1932780

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