Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 3,914,856 T→G 8.6% intergenic (‑70/‑176) PP_3454 ← / → PP_3455 putative transcriptional regulatory protein/Multidrug efflux RND membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154513,914,8560TG 8.6% 149.8 / 4.8 44intergenic (‑70/‑176)PP_3454/PP_3455putative transcriptional regulatory protein/Multidrug efflux RND membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base G (3/2);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCT  >  AE015451/3914635‑3915057
                                                                                                                                                                                                                             |                                                                                                                                                                                                         
aGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCgg                                                                                                                                                                                             <  1:711211/235‑1 (MQ=255)
                        cGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGcc                                                                                                                                                                          >  1:8276/1‑231 (MQ=255)
                        cGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCa                                                                                                                                                                         >  2:539513/1‑232 (MQ=255)
                        cGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCa                                                                                                                                                                         <  1:539514/232‑1 (MQ=255)
                        cGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAg                                                                                                                                                                        >  2:262439/1‑233 (MQ=255)
                        cGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGcac                                                                                                                                                                     >  1:716871/1‑236 (MQ=255)
                            ggTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAgg                                                                                                                                                                  <  1:408623/235‑1 (MQ=255)
                             gTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGGAAACCCGCTGCTGCAGGGCCCCg                                                                                                                                                                                                      >  1:150597/1‑198 (MQ=255)
                                 tCGTCCAGCGCATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGccc                                                                                                                                                               <  2:716870/233‑1 (MQ=255)
                                           caTGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAg                                                                                                                                                  >  1:420383/1‑236 (MQ=255)
                                               aTCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGa                                                                                                                                                 >  2:497257/1‑233 (MQ=38)
                                                    cccGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAAc                                                                                                                                          >  1:523147/1‑235 (MQ=255)
                                                              cACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGccc                                                                                                                                                               <  1:205582/204‑1 (MQ=255)
                                                                    gACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGccc                                                                                                                                                               <  1:714945/198‑1 (MQ=255)
                                                                          gAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAgg                                                                                                                                                                  <  2:8276/189‑1 (MQ=255)
                                                                            aGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                                                                   >  2:669029/1‑170 (MQ=255)
                                                                            aGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                                                                   <  1:669030/170‑1 (MQ=255)
                                                                              tCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACtt                                                                                                               >  1:82379/1‑236 (MQ=255)
                                                                              gcgttgcgttgCAGGTAGCTGGCGATAAGCTCGTACAGGGCGCTGTCGTCTTCGACCAGGAAAAGGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTccc                                                                                                                                                                                                        <  2:35902/146‑1 (MQ=255)
                                                                                 ttgcgttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCgg                                                                                                                                                                                             <  2:566136/155‑1 (MQ=255)
                                                                                 ttgcgttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                                                                            >  1:725569/1‑236 (MQ=255)
                                                                                    cgttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                                                                                         <  1:214878/236‑1 (MQ=255)
                                                                                    cgttgcAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                                                                                         <  1:456959/236‑1 (MQ=255)
                                                                                            gTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                                                                                                  <  1:497258/235‑1 (MQ=38)
                                                                                            gTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                                                                                                  <  1:429734/235‑1 (MQ=255)
                                                                                                   gCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                                                                                             <  1:577789/169‑1 (MQ=255)
                                                                                                        cAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                                                              <  2:479417/147‑1 (MQ=255)
                                                                                                              ggACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCaa                                                                                                                                           <  1:1166/176‑1 (MQ=255)
                                                                                                                            cgtcTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCt                                                                 >  1:326389/1‑236 (MQ=255)
                                                                                                                                tgcGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGGAAGCCTGTGGCGAGCGCTGCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGGTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACttttt                                                                                                            <  2:18768/187‑1 (MQ=255)
                                                                                                                                 tCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATc                                                                                                                            >  2:738439/1‑172 (MQ=255)
                                                                                                                                 tCGACCAGGAAAAGGGTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCGTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGcaca                                                                                                                                                                    <  2:34753/132‑1 (MQ=255)
                                                                                                                                     caaGGAGAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGGAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGGCCCCGTAGGAGCGGGGTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                                                                                             <  2:48146/133‑1 (MQ=38)
                                                                                                                                         gAAAATGTTCGGCAGAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAg                                                                      <  2:626025/218‑1 (MQ=255)
                                                                                                                                              tGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                                                                              <  1:313097/109‑1 (MQ=255)
                                                                                                                                                                  gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGTCATCATCGGAAGATCCGCa                           >  1:587125/1‑236 (MQ=255)
                                                                                                                                                                  gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCACGCAGCGTCTGCCTTAGCATATACACGGTCATCATCGGAAGATCCGCa                           >  1:100528/1‑236 (MQ=255)
                                                                                                                                                                       tGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGTCATCATCGGAAGATCCGCATGcc                       >  1:722528/1‑235 (MQ=255)
                                                                                                                                                                                           gCGAGCGCTTCGCGGGTGAACCCGCTCCTACGGGGACCTGCAGCAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGcc   >  1:710678/1‑235 (MQ=38)
                                                                                                                                                                                           gCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                                                                                         >  2:625998/1‑133 (MQ=255)
                                                                                                                                                                                           gCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                                                                                         <  1:625999/133‑1 (MQ=255)
                                                                                                                                                                                            cGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCt  >  1:476801/1‑235 (MQ=255)
                                                                                                                                                                                                           tAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                                                                                         >  1:552206/1‑117 (MQ=255)
                                                                                                                                                                                                           tAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                                                                                         <  2:552205/117‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                         
AGGATGACCAGGTCAGGCCTTTGCCGACGGTACTCGTCCAGCACATGATCGCCCCGGGCGATCACCTGGACATGGAAGTCGTTGCGTTGCAGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCAGGCAGCGTCTGCCTTAGCATATACACGGTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCT  >  AE015451/3914635‑3915057

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