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.C96.5% 180.7 / ‑0.3 57intergenic (‑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 . (2/0);  new base C (29/26);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932491‑1932781
                                                                                                                                                      |                                                                                                                                                
tGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGc                                                                                                                                                 <  1:259682/149‑1 (MQ=255)
  gcggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcc                                                                                                                                                <  2:379360/148‑1 (MQ=255)
  gcggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcc                                                                                                                                                <  2:374378/148‑1 (MQ=255)
    ggcgTCTTCGAAGAGGAGCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCCGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgc                                                                                                                                             <  5:260726/149‑4 (MQ=255)
     gcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgc                                                                                                                                             >  1:102607/1‑145 (MQ=255)
          ttCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAggg                                                                                                                                       >  8:200577/1‑149 (MQ=255)
               agagGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAGGGGCTGcc                                                                                                                                  >  8:110491/1‑148 (MQ=255)
                   gATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCGCCTACAAAGGATATGCGTCCCCGGTAGGAGCGGGGTTGCCCGCGAAGGGCTGCAAGGc                                                                                                                              >  3:32893/1‑149 (MQ=38)
                    aTCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCAGCGAAGGGCTGCAAgg                                                                                                                               >  7:303101/1‑147 (MQ=255)
                       tGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGcc                                                                                                                          >  2:110464/1‑149 (MQ=255)
                            aTTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCGACAAAGGATATGCGTCCCCTGTAGGAGCGGGCGTGGCCGCGAAGGGCTGCCAGGCAGCCCCaaa                                                                                                                     >  8:391904/1‑149 (MQ=38)
                             ttCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAGGGGCTGCCTTGCAGcccttcgc                                                                                                                    <  7:110491/149‑6 (MQ=37)
                                  cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGc                                                                                                                                   >  4:305458/1‑129 (MQ=255)
                                  cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGc                                                                                                                                   <  3:305458/129‑1 (MQ=255)
                                     cgcaaTTGCAGGTCTCACCAAAATAGCAGTAAACGCGGGGGCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCTCCTACAAAGCAAATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                            <  7:391904/145‑1 (MQ=38)
                                      ggCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTa                                                                                                           >  5:190482/1‑149 (MQ=255)
                                                  ctcACGAAAATCGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAACGCATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCg                                                                                               <  8:379457/149‑1 (MQ=38)
                                                   tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGa                                                                                              >  1:260723/1‑149 (MQ=255)
                                                       gAAAATAGCAGTAACCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTCGGCGCGGGCTTGCCCGCGAAGGGCTGCaa                                                                                                                                 <  2:265585/110‑1 (MQ=25)
                                                       gAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGGCCCTTCGCGGGCAACCCCGCTCCTACAAAGGATTTGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCaa                                                                                                                                 >  1:265585/1‑111 (MQ=37)
                                                        aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAATGATATGCGTCCCCTGTAGGAGCGGGCATGCCCGCCAAGGGCTGCAACGCAGCCCCATATTGTCGACTACGCGGGTCAGCGaaacgg                                                                                         >  1:101634/1‑144 (MQ=37)
                                                            tAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTGCCCGGTAGGAGCGGGCTTGCCCGCGAAGGGCTGGAAGGCAGCCCCAAATTGTCGCCTACGCGGATCCGCGATCCggagct                                                                                     >  3:58365/1‑145 (MQ=37)
                                                             aGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCCATCCGGCAggg                                                                                    >  7:225281/1‑149 (MQ=255)
                                                                  tAACCGCGGGGGCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                               <  8:225281/149‑1 (MQ=38)
                                                                         tGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCt                                                                        <  5:78515/149‑1 (MQ=255)
                                                                         gggggCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCATACAAATGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCt                                                                        <  2:102607/148‑1 (MQ=255)
                                                                                gCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGTATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAAt                                                                 <  3:183821/149‑1 (MQ=255)
                                                                                gCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAAt                                                                 <  6:190482/149‑1 (MQ=255)
                                                                                 ttttGCGGCCCTTCGCGGGCAATCCCGGTCCTACAAAGGATATGCGTCCCCTGTAGGAGCTGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                                <  3:217132/148‑1 (MQ=38)
                                                                                        ggCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCgg                                                                                                      >  8:322086/1‑105 (MQ=255)
                                                                                        ggCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCgg                                                                                                      <  7:322086/105‑1 (MQ=255)
                                                                                          cccTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCAGTATGAGCGGGCGTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAATGTCGACTACGCGGATAAGCGCTCCGGCAGGGTGATGTTGAGCTCCAGAATCGAGCAGCt                                                        >  6:2176/1‑149 (MQ=38)
                                                                                                gcgGGCTAGCCCGCTCCTACAAAGTATATGCGTGCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCt                                                   <  8:323950/148‑1 (MQ=39)
                                                                                                gcgGGCAAGCCCGCTCCTACAAAGGATATGCGTACCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTTTCGACTACGCAGAACAGCGATCCGGCAGGGTGATGTTAAGCTCCAGAATAGAACAACTGCCCTg                                                  >  8:183054/1‑149 (MQ=38)
                                                                                                  gggCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGccc                                                    <  4:233220/145‑1 (MQ=255)
                                                                                                  gggCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGccc                                                    >  3:233220/1‑145 (MQ=255)
                                                                                                  gggCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGcc                                                     >  3:141030/1‑144 (MQ=255)
                                                                                                  gggCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGcc                                                     <  4:141030/144‑1 (MQ=255)
                                                                                                      aaGCCCGCTCCTACAATGTATATGCGTCCCCTTTAGGCGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                            <  5:2176/149‑1 (MQ=255)
                                                                                                      aaGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                            >  6:100195/1‑149 (MQ=255)
                                                                                                        gCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCa                                          <  1:110464/149‑1 (MQ=255)
                                                                                                         cccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                                         >  1:202086/1‑149 (MQ=255)
                                                                                                          ccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGc                                        >  6:432610/1‑149 (MQ=255)
                                                                                                           cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCt                                       >  7:129559/1‑149 (MQ=255)
                                                                                                             cTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCTGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCg                                     <  8:303101/149‑1 (MQ=255)
                                                                                                             cTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCg                                     <  2:202086/149‑1 (MQ=255)
                                                                                                              tCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagc                                                            <  1:164177/125‑1 (MQ=255)
                                                                                                              tCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa                                    >  5:194131/1‑149 (MQ=255)
                                                                                                                 tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                            <  7:124728/138‑1 (MQ=255)
                                                                                                                 tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                            >  8:124728/1‑138 (MQ=255)
                                                                                                                   cAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                                 <  7:200577/147‑1 (MQ=255)
                                                                                                                    aatggATATGCTTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGca                              <  2:101634/146‑1 (MQ=255)
                                                                                                                        gATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGAAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAAACGAGCAGCTGCCCAGGTATTCCAGCTCGATATGCAcat                           >  2:64255/1‑148 (MQ=255)
                                                                                                                                tCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCCGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                  >  8:329427/1‑149 (MQ=255)
                                                                                                                                  cccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGcc                 >  2:146946/1‑148 (MQ=255)
                                                                                                                                    cTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt              >  8:60795/1‑149 (MQ=255)
                                                                                                                                     tGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATg             >  3:76485/1‑149 (MQ=255)
                                                                                                                                      gtcggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGt            <  3:117274/146‑1 (MQ=255)
                                                                                                                                       tAGGAGCGGGCTTGCCCGCGAAGGGCTGCAACGCAGCCCCAAATTGTCGAATACGCGGATCACCGATCCCGCAGGGTGATCTTCAGCTCCAGAATCGACCAGCTGACCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGt            >  2:222012/1‑148 (MQ=38)
                                                                                                                                         ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca         <  5:100195/149‑1 (MQ=255)
                                                                                                                                         ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca         >  8:69827/1‑149 (MQ=255)
                                                                                                                                                gCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAttt  <  8:129559/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932491‑1932781

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: