Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,260 A→T 5.3% intergenic (‑37/‑156) grpE ← / → recN heat shock protein GrpE/DNA repair/recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,2600AT 5.3% 236.0 / 2.3 76intergenic (‑37/‑156)grpE/recNheat shock protein GrpE/DNA repair/recombination protein
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base T (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

GGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGA  >  NC_002947/5378122‑5378406
                                                                                                                                                     |                                                                                                                                                  
ggCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCtt                                                                                                                                          >  3:388749/1‑149 (MQ=255)
  ccAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCCTGCGCAATTTTGTGGACTGCGAGCCTGGCTCGcagtc                                                                                                                                         <  1:73379/149‑5 (MQ=255)
  ccAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTCTTGTGGACTGCGAGCCAGGCTCGCCTTg                                                                                                                                         <  1:24651/149‑1 (MQ=255)
     gCTGTTCCTCGAGCTCCAGAAGACGTGTATCAGCTGCATATCATGCACCTGTTTCTCACACTCACTACCCCNCCAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCTGGCTCg                                                                                                                                              <  1:131590/149‑1 (MQ=255)
         ttCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTc                                                                                                                                   <  2:113003/147‑1 (MQ=255)
            cTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTg                                                                                                                                         <  6:358848/138‑1 (MQ=255)
            cTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTg                                                                                                                                         >  5:358848/1‑138 (MQ=255)
            cTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCg                                                                                                                                  >  7:240027/1‑145 (MQ=255)
            cTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCg                                                                                                                                  <  8:240027/145‑1 (MQ=255)
             tCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCta                                                                                                                              >  2:202408/1‑148 (MQ=255)
                   tCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATgggg                                                                                                                       >  4:332329/1‑149 (MQ=255)
                    ccAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATTGGGc                                                                                                                      >  5:365609/1‑149 (MQ=255)
                     cAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCACGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGcc                                                                                                                     >  1:113258/1‑149 (MQ=255)
                               gTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAAc                                                                                                             >  2:191786/1‑147 (MQ=255)
                               gTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCa                                                                                                           <  3:135977/149‑1 (MQ=255)
                                    aGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGtt                                                                                                       >  7:103170/1‑148 (MQ=255)
                                      cTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTgg                                                                                                                                        <  5:420702/113‑1 (MQ=255)
                                      cTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTc                                                                                                      >  1:105737/1‑147 (MQ=255)
                                      cTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATAAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTgg                                                                                                                                        >  6:420702/1‑113 (MQ=255)
                                         cATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAgg                                                                                                  <  7:371989/148‑1 (MQ=255)
                                         cATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAg                                                                                                   <  2:105737/147‑1 (MQ=255)
                                               cTGCACCTGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAgg                                                                                           <  5:344497/149‑1 (MQ=255)
                                               cTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAgg                                                                                           >  4:327412/1‑149 (MQ=255)
                                                       gTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGtttt                                                                                    >  4:345039/1‑148 (MQ=255)
                                                          tcttcCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCaa                                                                                 <  4:388749/148‑1 (MQ=255)
                                                              ccA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACt                                                                             <  1:191786/148‑1 (MQ=255)
                                                                a‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACt                                                                             >  6:204409/1‑146 (MQ=255)
                                                                a‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACt                                                                             <  5:204409/146‑1 (MQ=255)
                                                                   cGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCaa                                                                         <  8:103170/149‑1 (MQ=255)
                                                                                    tttCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGc                                                               >  6:52459/1‑149 (MQ=255)
                                                                                      tctcGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGc                                                               <  3:345039/147‑1 (MQ=255)
                                                                                       ctcGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCaa                                                                                 <  5:85690/128‑1 (MQ=255)
                                                                                       ctcGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCaa                                                                                 >  6:85690/1‑128 (MQ=255)
                                                                                       ctcGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCg                                                              >  5:190767/1‑147 (MQ=255)
                                                                                          gTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACa                                                                       >  2:16576/1‑135 (MQ=255)
                                                                                          gTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACa                                                                       <  1:16576/135‑1 (MQ=255)
                                                                                          gTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAggg                                                         <  1:202408/149‑1 (MQ=255)
                                                                                          gTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAggg                                                         <  6:190767/149‑1 (MQ=255)
                                                                                          gTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAggg                                                         <  3:332329/149‑1 (MQ=255)
                                                                                              aGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAgg                                                                                                         <  2:399040/97‑1 (MQ=255)
                                                                                                ctgctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTg                                                   >  4:376385/1‑149 (MQ=255)
                                                                                                ctgctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTg                                                   >  3:399238/1‑149 (MQ=255)
                                                                                                   ctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGc                                                 <  3:376385/148‑1 (MQ=255)
                                                                                                      atcaGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAgg                                              >  6:228593/1‑148 (MQ=255)
                                                                                                      atcaGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGc                                             >  5:263906/1‑149 (MQ=255)
                                                                                                      atcaGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGGTCAAGGGGCAAGGTAGTTTTCAAGACTGCAAAATAATATGCGAAGGGGCATTGGCAGGc                                             >  7:187683/1‑149 (MQ=255)
                                                                                                            cATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaa                                        <  6:365609/148‑1 (MQ=255)
                                                                                                             aTGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGTCTATATGGGGCCGGAAAAACCAGGTTCAAgg                                                                                                  <  7:359773/89‑1 (MQ=255)
                                                                                                             aTGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAgg                                                                                                  >  8:359773/1‑89 (MQ=255)
                                                                                                             aTGGGGTCTCTCCTGAGCAATTTTGAGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaa                                       >  2:336412/1‑148 (MQ=255)
                                                                                                                gggTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaaaa                                    <  5:52459/148‑1 (MQ=255)
                                                                                                                gggTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaaaa                                    <  3:327412/148‑1 (MQ=255)
                                                                                                                     tctcCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTg                              >  8:116292/1‑149 (MQ=255)
                                                                                                                       tcCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGta                            >  3:263494/1‑149 (MQ=255)
                                                                                                                             cgcAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaa                                      >  6:112041/1‑134 (MQ=255)
                                                                                                                             cgcAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaa                                      <  5:112041/134‑1 (MQ=255)
                                                                                                                              gcAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTataaat                       >  3:331356/1‑148 (MQ=255)
                                                                                                                                aaTTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAAc                    <  2:113258/149‑1 (MQ=255)
                                                                                                                                  ttttGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCCATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTataaa                        >  1:157456/1‑143 (MQ=255)
                                                                                                                                  ttttGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCCATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTataaa                        <  2:157456/143‑1 (MQ=255)
                                                                                                                                           aaggCGAGCCTGGCTCGCCTTGGATGTCGGCTATATGGCGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCa                  >  7:37234/4‑141 (MQ=255)
                                                                                                                                            aggCGAGCCTGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCATGGGGGAAGGGAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCAATGGCATGCCTGAAAAAAGTTCTGTATAAATAAACAGACGTGAtt         >  5:246608/3‑149 (MQ=255)
                                                                                                                                             tGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACGAGCCCTGAttt        >  6:215226/1‑149 (MQ=255)
                                                                                                                                             tGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGAttt        >  6:318237/1‑149 (MQ=255)
                                                                                                                                             tGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGAttt        >  5:23872/1‑149 (MQ=255)
                                                                                                                                               cGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaaaa                                    <  1:361877/119‑1 (MQ=255)
                                                                                                                                               cGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaaaa                                    >  2:361877/1‑119 (MQ=255)
                                                                                                                                                gAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaaaa                                    <  2:396190/118‑1 (MQ=255)
                                                                                                                                                gAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaaaa                                    >  1:396190/1‑118 (MQ=255)
                                                                                                                                                gAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaa                                      <  3:210153/116‑1 (MQ=255)
                                                                                                                                                gAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaa                                      >  4:210153/1‑116 (MQ=255)
                                                                                                                                                gAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCg     <  7:116292/149‑1 (MQ=255)
                                                                                                                                                  gCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCggg   >  8:210291/1‑149 (MQ=255)
                                                                                                                                                  gCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCgg    >  6:173305/1‑148 (MQ=255)
                                                                                                                                                   ccAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGa  <  4:331356/149‑1 (MQ=255)
                                                                                                                                                   ccAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGa  <  4:399238/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCA‑‑CGTTA‑‑‑‑‑‑‑AGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGA  >  NC_002947/5378122‑5378406

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: