Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,591,982 C→A 20.3% intergenic (‑89/‑114) PP_1395 ← / → PP_1396 AraC family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,591,9820CA20.3% 124.6 / 25.7 63intergenic (‑89/‑114)PP_1395/PP_1396AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base A (8/5);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

GCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTT  >  NC_002947/1591838‑1592129
                                                                                                                                                |                                                                                                                                                   
gCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCgcgc                                                                                                                                                 >  3:213569/1‑149 (MQ=255)
  cAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTTGATTGCCTGATCATGCCGAGCTg                                                                                                                                               >  8:105537/1‑149 (MQ=255)
  cAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGAAGGCTTTCCTCTGCTGACCGCCGCAGCATAGGTTTGGGCAAAGCCGGGCGCTCGTCGCAAGCTGCATGTGGATTGCCTGATCCTGCCgcgcgg                                                                                                                                               >  8:109890/1‑147 (MQ=255)
   aGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATATGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGc                                                                                                                                              <  7:68732/149‑1 (MQ=255)
        cAGGCGCTGCAGTTCCAGGGGCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGAt                                                                                                                                         <  8:208579/149‑1 (MQ=255)
           gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAg                                                                                                                                      <  5:199797/149‑1 (MQ=255)
           gcgcTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAg                                                                                                                                      <  8:17878/149‑1 (MQ=255)
                       cAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACAt                                                                                                                           <  4:281407/148‑1 (MQ=255)
                       cAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACAt                                                                                                                           >  3:281407/1‑148 (MQ=255)
                       cAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATg                                                                                                                          >  4:161973/1‑149 (MQ=255)
                       cAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATg                                                                                                                          >  3:233650/1‑149 (MQ=255)
                       cAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACAt                                                                                                                           <  4:54231/148‑1 (MQ=255)
                                ttcggtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCATAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATgg                                                                                                                 <  4:213569/149‑1 (MQ=255)
                                  cggtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGc                                                                                                                <  4:233650/148‑1 (MQ=255)
                                       gACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATc                                                                                                                                <  8:659/127‑1 (MQ=255)
                                       gACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATc                                                                                                                                >  7:659/1‑127 (MQ=255)
                                               ggtggtCAAGATGGCTTTCATCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAAcgcgccgcgc                                                                                                  >  8:45322/1‑149 (MQ=255)
                                                        gatgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAgcg                                                                                         >  8:103648/1‑149 (MQ=255)
                                                            gccttCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCt                                                                                      <  3:267343/145‑1 (MQ=255)
                                                                tCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAAc                                                                                   <  8:157009/147‑1 (MQ=255)
                                                                    ctgctgACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGGAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACACGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTgg                                                                              >  7:161028/1‑148 (MQ=255)
                                                                            cGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGc                                                                                                                >  8:48901/1‑106 (MQ=255)
                                                                            cGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGc                                                                                                                <  7:48901/106‑1 (MQ=255)
                                                                            cGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGc                                                                     >  7:242375/1‑149 (MQ=255)
                                                                                gcagcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCACGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCtttt                                                                 >  1:253796/1‑149 (MQ=255)
                                                                                 cagcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtc                                                                <  7:45322/149‑1 (MQ=255)
                                                                                         gTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGcc                                                        >  7:244375/1‑149 (MQ=255)
                                                                                         gTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGcc                                                        >  8:18946/1‑149 (MQ=255)
                                                                                               gCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGTTCCTGGCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGGCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTgg                                                                              <  4:98416/121‑1 (MQ=255)
                                                                                               gCAAACACGGGCGCTAGTCGAAAGATGCATGTGCATTCCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTgg                                                                              >  3:98416/1‑121 (MQ=255)
                                                                                                cAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCAGTCCGCCTGCTTc                                                  >  1:62756/1‑148 (MQ=255)
                                                                                                 aaaCACGGTCGCTAGTCGATAGCTGCATGCGGATTGCCTGATCCTGCCTCGATGCAGGATCAGGCAATCGACATGc                                                                                                                         <  5:282006/76‑1 (MQ=255)
                                                                                                  aaGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAg                                               <  7:107032/149‑1 (MQ=255)
                                                                                                  aaGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAg                                               <  7:159174/149‑1 (MQ=255)
                                                                                                          cgcTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGAc                                                                                                                             <  6:257133/63‑1 (MQ=255)
                                                                                                          cgcTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGAc                                                                                                                             >  5:257133/1‑63 (MQ=255)
                                                                                                           gcTAGTCGGAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACACGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATtgtg                                      <  8:161028/149‑1 (MQ=255)
                                                                                                            cgaggcGAAAGCTGCATGTGGATTGCGTGATCCTGCCGCGCTGCAGGATGAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCACGGAGATTGTGc                                     <  7:109890/144‑1 (MQ=255)
                                                                                                             tAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGACTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGga                                             >  6:242027/1‑140 (MQ=255)
                                                                                                             tAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGACTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGga                                             <  5:242027/140‑1 (MQ=255)
                                                                                                                gccattacTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAAt                                  >  5:326092/8‑148 (MQ=37)
                                                                                                                  cattacTGCATGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACATGCAGGAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGGGCAATgcc                               >  6:41495/6‑147 (MQ=25)
                                                                                                                   aaaGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGAcc                              >  5:295002/1‑149 (MQ=255)
                                                                                                                      gCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCt                            >  1:191699/1‑148 (MQ=255)
                                                                                                                           aTGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTa                          >  7:174141/1‑145 (MQ=255)
                                                                                                                           aTGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTa                          <  8:174141/145‑1 (MQ=255)
                                                                                                                           aTGTCGATTGCCTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaac                       >  5:181008/1‑148 (MQ=255)
                                                                                                                              tgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaacaac                    <  7:18946/148‑1 (MQ=255)
                                                                                                                                gATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCg                  <  3:161973/148‑1 (MQ=255)
                                                                                                                                 aTTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGAccc                             >  1:67924/1‑136 (MQ=255)
                                                                                                                                 aTTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGAccc                             <  2:67924/136‑1 (MQ=255)
                                                                                                                                     ccTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTaa                                                                                                            >  6:240552/1‑53 (MQ=255)
                                                                                                                                     ccTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTaa                                                                                                            <  5:240552/53‑1 (MQ=255)
                                                                                                                                      cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccc           <  8:242375/149‑1 (MQ=255)
                                                                                                                                      cTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacc            >  6:309594/1‑148 (MQ=255)
                                                                                                                                      cTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGAcc                              >  3:46944/1‑130 (MQ=255)
                                                                                                                                      cTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGAcc                              <  4:46944/130‑1 (MQ=255)
                                                                                                                                      cTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaac                       <  4:292878/137‑1 (MQ=255)
                                                                                                                                      cTGATCCTGCAGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaac                       >  3:292878/1‑137 (MQ=255)
                                                                                                                                      cTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccc           >  1:43906/1‑149 (MQ=255)
                                                                                                                                      cTGATCCTGCAGCGCGGCAGGATCAGGCAATCCACATGCAGCTTTCGACTA‑‑GCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcacccac         >  2:170842/1‑149 (MQ=25)
                                                                                                                                          tCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTaaa                        <  1:144637/132‑1 (MQ=255)
                                                                                                                                          tCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTaaa                        >  2:144637/1‑132 (MQ=255)
                                                                                                                                               ccGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCtt  <  2:45579/149‑1 (MQ=255)
                                                                                                                                                cgcgCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCtt  <  7:103648/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCCAGCTCCAGGCGCTGCAGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCCACCTTGCTT  >  NC_002947/1591838‑1592129

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