Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,747,477 G→A intergenic (+69/+659) b4540 → / ← grpE hypothetical protein/GrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,747,4770GA86.6% 747.9 / 76.5 278intergenic (+69/+659)b4540/grpEhypothetical protein/GrpE
Reads supporting (aligned to +/- strand):  ref base G (20/17);  new base A (132/107);  total (152/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-03

Only 100 of 280 total aligned reads displayed.

TCCTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGAAGCCTG  >  NC_000913/2747387‑2747567
                                                                                          |                                                                                           
atcTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCAcc                                                                                           <  2:1709198/91‑1 (MQ=255)
atcTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCAcc                                                                                           >  2:1334342/3‑93 (MQ=255)
atcTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCAcc                                                                                           >  2:5440200/3‑93 (MQ=255)
atcTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCAcc                                                                                           >  2:5272055/3‑93 (MQ=255)
  cTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCaa                                                                                         <  1:2605855/93‑1 (MQ=255)
   tGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCGCCAAc                                                                                        >  1:7019495/1‑93 (MQ=255)
   tGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAAc                                                                                        >  2:2145918/1‑93 (MQ=255)
     gAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAActc                                                                                      <  2:331641/93‑1 (MQ=255)
       aaaaaTGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACtctc                                                                                   >  2:2912667/1‑93 (MQ=255)
         aaaTGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTccc                                                                                 >  2:661254/1‑93 (MQ=255)
          aaTGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCa                                                                                <  2:891950/93‑1 (MQ=255)
          aaTGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACCGGAAATCGATTCTGCACCAACTCTCCCa                                                                                <  1:3265748/93‑1 (MQ=255)
           aTGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCAc                                                                               <  1:5994461/93‑1 (MQ=255)
            tGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCAct                                                                              >  1:224230/1‑93 (MQ=255)
            tGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCAct                                                                              >  1:3803337/1‑93 (MQ=255)
             gACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCActc                                                                             >  1:2822518/1‑93 (MQ=255)
             gACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCActc                                                                             <  2:1252419/93‑1 (MQ=255)
             gACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCActc                                                                             >  1:7407073/1‑93 (MQ=255)
             gACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCActc                                                                             >  2:6818906/1‑93 (MQ=255)
              aCCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCActca                                                                            <  2:5334985/93‑1 (MQ=255)
               ccGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCActcac                                                                           >  1:6083687/1‑93 (MQ=255)
                  ttGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCActcacctc                                                                        >  2:1981539/1‑93 (MQ=255)
                   tgatgaCATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCActcacctca                                                                       <  2:5947111/93‑1 (MQ=255)
                    gatgaCATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCActcacctcac                                                                      <  1:3746690/93‑1 (MQ=255)
                       gaCATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACggg                                                                   <  1:112929/93‑1 (MQ=255)
                           tcatGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCa                                                               >  1:3689729/1‑93 (MQ=255)
                            catGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCAt                                                              >  1:6674519/1‑93 (MQ=255)
                             atGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATc                                                             <  2:641929/93‑1 (MQ=255)
                             atGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATc                                                             >  1:4583322/1‑93 (MQ=255)
                             atGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATc                                                             >  1:7016272/1‑93 (MQ=255)
                              tGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCg                                                            >  1:1244908/1‑93 (MQ=255)
                               ggTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATTACTGGGAATCGATTCTGCACCAACTCTCCCCCCTCACCTCACGGGCGCATCg                                                            <  1:6157336/93‑1 (MQ=255)
                               ggTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGt                                                           <  2:3575699/93‑1 (MQ=255)
                                    cgcgCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGctct                                                      <  2:518747/93‑1 (MQ=255)
                                     gcgcAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTa                                                     <  2:6370070/93‑1 (MQ=255)
                                      cgcAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTAc                                                    <  2:6281034/93‑1 (MQ=255)
                                      cgcAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTAc                                                    <  2:960401/93‑1 (MQ=255)
                                       gcAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTAcc                                                   >  2:1957072/1‑93 (MQ=255)
                                       gcAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTAcc                                                   <  2:491831/93‑1 (MQ=255)
                                       gcAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTAcc                                                   >  2:1487008/1‑93 (MQ=255)
                                       gcAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTAcc                                                   >  1:4473327/1‑93 (MQ=255)
                                        cAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCg                                                  >  2:1562006/1‑93 (MQ=255)
                                        cAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCg                                                  >  2:4369228/1‑93 (MQ=255)
                                          gTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTg                                                >  1:1342401/1‑93 (MQ=255)
                                          gTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCACACTCACCTCACGGGCGCATCGTGCTCTACCGTg                                                >  2:1830096/1‑93 (MQ=255)
                                            gAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGAt                                              <  2:3921807/93‑1 (MQ=255)
                                            gAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGAt                                              >  1:5665773/1‑93 (MQ=255)
                                              aaTTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTc                                            <  2:3888670/93‑1 (MQ=255)
                                                ttATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGc                                          <  1:1895129/93‑1 (MQ=255)
                                                 tATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCt                                         <  1:997549/93‑1 (MQ=255)
                                                 tATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCt                                         <  2:4648953/93‑1 (MQ=255)
                                                  aTCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTg                                        >  2:3813895/1‑93 (MQ=255)
                                                  aTCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTg                                        >  1:3177368/1‑93 (MQ=255)
                                                  aTCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTg                                        <  2:2666275/93‑1 (MQ=255)
                                                  aTCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTg                                        >  2:259767/1‑93 (MQ=255)
                                                   tCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTgg                                       >  1:934125/1‑93 (MQ=255)
                                                    cGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGga                                      <  1:4980023/93‑1 (MQ=255)
                                                    cGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGga                                      <  1:6889733/93‑1 (MQ=255)
                                                    cGGTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGga                                      <  2:4116329/93‑1 (MQ=255)
                                                      gTATTGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGgacg                                    >  2:7329138/1‑93 (MQ=255)
                                                         ttGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGacgacg                                 >  2:7465140/1‑93 (MQ=255)
                                                          tGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGc                                <  1:6483914/93‑1 (MQ=255)
                                                          tGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGc                                >  1:5311554/1‑93 (MQ=255)
                                                          tGATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGCACGACGc                                >  2:2318840/1‑93 (MQ=255)
                                                           gATATCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGcc                               <  1:4373439/93‑1 (MQ=255)
                                                            atatCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCa                              <  2:5705616/93‑1 (MQ=255)
                                                              atCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATc                            <  2:142611/93‑1 (MQ=255)
                                                               tCAACGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCa                           <  1:2069190/93‑1 (MQ=255)
                                                               tCAACGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCa                           >  2:6444736/1‑93 (MQ=255)
                                                                 aaCGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGt                         >  1:6572791/1‑93 (MQ=255)
                                                                  aCGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTa                        <  1:2385347/93‑1 (MQ=255)
                                                                   cGATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTAt                       <  1:7214208/93‑1 (MQ=255)
                                                                    gATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTAtg                      >  1:4081980/1‑93 (MQ=255)
                                                                    gATGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTAtg                      >  1:4252636/1‑93 (MQ=255)
                                                                     aTGACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTAtgc                     <  2:6102382/93‑1 (MQ=255)
                                                                       gACTGGAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTAtgctg                   <  1:1758144/93‑1 (MQ=255)
                                                                           ggAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAg               >  2:7211021/1‑93 (MQ=255)
                                                                           ggAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAg               <  1:3377554/93‑1 (MQ=255)
                                                                           ggAAATCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAg               <  2:6756685/93‑1 (MQ=255)
                                                                              aaTCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTAc            <  2:813578/93‑1 (MQ=255)
                                                                              aaTCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTAc            <  1:3640789/93‑1 (MQ=255)
                                                                               aTCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACg           >  1:3110036/1‑93 (MQ=255)
                                                                               aTCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACg           >  1:153493/1‑93 (MQ=255)
                                                                               aTCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACg           <  1:1103919/93‑1 (MQ=255)
                                                                               aTCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACg           <  1:1847300/93‑1 (MQ=255)
                                                                               aTCGATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACg           <  1:198330/93‑1 (MQ=255)
                                                                                  gATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGa        >  2:6057355/1‑93 (MQ=255)
                                                                                  gATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGa        <  1:5470321/93‑1 (MQ=255)
                                                                                  gATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGa        >  1:5747309/1‑93 (MQ=255)
                                                                                  gATTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGa        <  1:643240/93‑1 (MQ=255)
                                                                                   aTTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGaa       <  1:3882589/93‑1 (MQ=255)
                                                                                   aTTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGaa       <  2:2786440/93‑1 (MQ=255)
                                                                                   aTTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGaa       >  2:6648981/1‑93 (MQ=255)
                                                                                   aTTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGaa       >  1:7509786/1‑93 (MQ=255)
                                                                                   aTTCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGaa       >  1:4667750/1‑93 (MQ=255)
                                                                                    ttCTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGAAg      >  1:5113370/1‑93 (MQ=255)
                                                                                     tCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGAAGc     <  2:5620775/93‑1 (MQ=255)
                                                                                      cTGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGAAGcc    <  2:6521278/93‑1 (MQ=255)
                                                                                       tGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGAAGCCt   <  1:1427373/93‑1 (MQ=255)
                                                                                       tGCACCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGAAGCCt   >  1:4436011/1‑93 (MQ=255)
                                                                                          |                                                                                           
TCCTGGAAAAAATGACCGTTGATGACATCATGGTGCCGCGCAGTGAAATTATCGGTATTGATATCAACGATGACTGGAAATCGATTCTGCGCCAACTCTCCCACTCACCTCACGGGCGCATCGTGCTCTACCGTGATTCGCTGGACGACGCCATCAGTATGCTGCGAGTACGTGAAGCCTG  >  NC_000913/2747387‑2747567

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