Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,882,902 Δ2 bp intergenic (‑96/+520) atpI ← / ← gidB F0F1 ATP synthase subunit I/glucose‑inhibited division protein; b3740_1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,882,9000T.99.3% 1403.6 / 1.9 276intergenic (‑94/+523)atpI/gidBF0F1 ATP synthase subunit I/glucose‑inhibited division protein; b3740_1
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (133/141);  total (134/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01
*REL6063,882,9010A.99.3% 1401.8 / 1.8 276intergenic (‑95/+522)atpI/gidBF0F1 ATP synthase subunit I/glucose‑inhibited division protein; b3740_1
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (133/141);  total (134/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 306 total aligned reads displayed.

TCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCC  >  REL606/3882825‑3882969
                                                                           ||                                                                    
tCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAtt                                                                       <  1:3320804/76‑1 (MQ=255)
tCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAtt                                                                       <  1:3027548/76‑1 (MQ=255)
 cccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATta                                                                      <  2:2019609/76‑1 (MQ=255)
 cccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATta                                                                      <  2:2322428/76‑1 (MQ=255)
 cccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATta                                                                      <  1:1913355/76‑1 (MQ=255)
 cccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATta                                                                      >  1:6972597/1‑76 (MQ=255)
 cccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATta                                                                      >  1:6420057/1‑76 (MQ=255)
 cccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATta                                                                      <  1:5275383/76‑1 (MQ=255)
  ccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATtac                                                                     >  1:2217032/1‑75 (MQ=255)
  ccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATtac                                                                     <  1:6081868/76‑2 (MQ=255)
         ggTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTgcgc                                                            >  2:5553291/1‑76 (MQ=255)
           tATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGccc                                                          >  2:446654/1‑76 (MQ=255)
           tATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGccc                                                          >  2:3355599/1‑76 (MQ=255)
            aTGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGcccc                                                         >  1:4498478/1‑76 (MQ=255)
             tGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGccccc                                                        <  2:3789219/76‑1 (MQ=255)
              gCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCg                                                       <  1:4184735/76‑1 (MQ=255)
               ccGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAATCGGTCAAAT‑‑TACGGTGCGCCCCCGt                                                      >  1:6155757/1‑76 (MQ=255)
               ccGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGt                                                      <  2:1192347/76‑1 (MQ=255)
                cgcgTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTg                                                     <  1:7515938/76‑1 (MQ=255)
                  cgTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGAt                                                   >  2:3529968/1‑76 (MQ=255)
                    tgtCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGAttt                                                 >  1:5327866/1‑76 (MQ=255)
                     gtCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTc                                                <  1:5675894/76‑1 (MQ=255)
                     gtCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTc                                                >  1:6884100/1‑76 (MQ=255)
                     gtCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTc                                                >  2:48636/1‑76 (MQ=255)
                     gtCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTc                                                >  2:1995941/1‑76 (MQ=255)
                     gtCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTc                                                <  1:4675435/76‑1 (MQ=255)
                     gtCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTc                                                <  1:3550780/76‑1 (MQ=255)
                       cGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCaa                                              >  1:12139/1‑76 (MQ=255)
                       cGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCaa                                              <  2:578925/76‑1 (MQ=255)
                        gTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCaaa                                             <  2:544416/76‑1 (MQ=255)
                         tataAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAAc                                            <  1:1995919/76‑1 (MQ=255)
                         tataAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAAc                                            >  2:4321089/1‑76 (MQ=255)
                         tataAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAAc                                            >  1:3364453/1‑76 (MQ=255)
                         tataAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAAc                                            >  2:2993820/1‑76 (MQ=255)
                          ataAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACa                                           >  1:371979/1‑76 (MQ=255)
                          ataAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACa                                           >  1:5798458/1‑76 (MQ=255)
                          ataAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACa                                           >  1:1833897/1‑76 (MQ=255)
                          ataAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACa                                           >  1:4850352/1‑76 (MQ=255)
                           taATACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACaa                                          >  1:2059541/1‑76 (MQ=255)
                           taAAACTTTCTTTAAGGCTTAGATTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACaa                                          <  2:7514054/76‑1 (MQ=255)
                           taAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACaa                                          >  2:5194104/1‑76 (MQ=255)
                            aaaaCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAAt                                         >  1:2878076/1‑76 (MQ=255)
                            aaaaCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAAt                                         <  2:4362418/76‑1 (MQ=255)
                            aaaaCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAAt                                         <  1:7051213/76‑1 (MQ=255)
                            aaaaCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAAt                                         >  2:1916294/1‑76 (MQ=255)
                             aaaCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATa                                        <  1:4633602/76‑1 (MQ=255)
                              aaCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATaa                                       <  2:5021323/76‑1 (MQ=255)
                                ctttctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGt                                     <  1:1009704/76‑1 (MQ=255)
                                ctttctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGt                                     >  1:424933/1‑76 (MQ=255)
                                 tttctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTa                                    <  1:664487/76‑1 (MQ=255)
                                 tttctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTa                                    >  1:3970396/1‑76 (MQ=255)
                                 tttctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTa                                    >  1:7298719/1‑76 (MQ=255)
                                 tttctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTa                                    >  1:2291491/1‑76 (MQ=255)
                                  ttctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAg                                   <  2:1267478/76‑1 (MQ=255)
                                  ttctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAg                                   >  2:6549063/1‑76 (MQ=255)
                                   tctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGc                                  >  1:2438144/1‑76 (MQ=255)
                                    ctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGcc                                 <  1:5506040/76‑1 (MQ=255)
                                     tttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCa                                <  1:1456264/76‑1 (MQ=255)
                                     tttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCa                                <  1:5146711/76‑1 (MQ=255)
                                      ttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaa                               <  2:5913731/76‑1 (MQ=255)
                                      ttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaa                               >  1:761363/1‑76 (MQ=255)
                                       tAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaa                              <  2:6852009/76‑1 (MQ=255)
                                       tAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaa                              >  2:6274606/1‑76 (MQ=255)
                                        aaGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaaa                             >  2:7479596/1‑76 (MQ=255)
                                        aaGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaaa                             <  2:7697213/76‑1 (MQ=255)
                                        aaGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaaa                             <  1:2120348/76‑1 (MQ=255)
                                        aaGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaaa                             >  1:3160297/1‑76 (MQ=255)
                                        aaGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaaa                             <  1:4439821/76‑1 (MQ=255)
                                         aGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTAT‑‑GGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAg                            >  1:226771/1‑76 (MQ=255)
                                         aGGCTTAGAATCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAg                            <  2:6602864/76‑1 (MQ=255)
                                          ggCTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAgg                           <  2:5105965/76‑1 (MQ=255)
                                            cTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTg                         >  1:1113967/1‑76 (MQ=255)
                                            cTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTg                         <  1:1357819/76‑1 (MQ=255)
                                            cTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTg                         >  2:2430018/1‑76 (MQ=255)
                                            cTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTg                         >  2:1459281/1‑76 (MQ=255)
                                            cTTAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTg                         <  1:4762376/76‑1 (MQ=255)
                                             ttAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGa                        >  1:500367/1‑76 (MQ=255)
                                             ttAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGa                        <  2:1731667/76‑1 (MQ=255)
                                              tAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGaa                       >  2:2642266/1‑76 (MQ=255)
                                              tAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGaa                       >  2:2823235/1‑76 (MQ=255)
                                              tAGAGTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGaa                       >  1:3452032/1‑76 (MQ=255)
                                                gagTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATa                     <  1:5083707/76‑1 (MQ=255)
                                                gagTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATa                     <  2:4109968/76‑1 (MQ=255)
                                                gagTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATa                     <  2:1027071/76‑1 (MQ=255)
                                                 agTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATaa                    <  1:5742066/76‑1 (MQ=255)
                                                 agTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATaa                    <  1:1660164/76‑1 (MQ=255)
                                                 agTCAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATaa                    <  1:2531897/76‑1 (MQ=255)
                                                    cAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATg                 <  1:238626/76‑1 (MQ=255)
                                                    cAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATg                 >  1:1570589/1‑76 (MQ=255)
                                                    cAAGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATg                 >  1:741568/1‑76 (MQ=255)
                                                     aaGCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGt                >  2:7316833/1‑76 (MQ=255)
                                                       gCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGttt              <  2:6056647/76‑1 (MQ=255)
                                                       gCATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGttt              >  2:2155749/1‑76 (MQ=255)
                                                        cATCAAAAAGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTa             <  2:7027541/76‑1 (MQ=255)
                                                            aaaaaGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAata         >  2:4609842/1‑76 (MQ=255)
                                                             aaaaGCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAatat        <  1:6547577/76‑1 (MQ=255)
                                                                 gCGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAttttt    <  2:4026387/76‑1 (MQ=255)
                                                                  cGGTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTc   <  2:3147953/76‑1 (MQ=255)
                                                                   ggTCAAAT‑‑TACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTcc  <  1:5990975/76‑1 (MQ=255)
                                                                     tCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTcc  <  2:1628622/76‑1 (MQ=255)
                                                                           ||                                                                    
TCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCC  >  REL606/3882825‑3882969

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