Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,961,966 A→G intergenic (‑122/+14) gppA ← / ← rhlB guanosine pentaphosphatase/putative ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,961,9660AG99.9% 2151.1 / ‑5.4 653intergenic (‑122/+14)gppA/rhlBguanosine pentaphosphatase/putative ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (306/345);  total (307/346)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

Only 100 of 670 total aligned reads displayed.

GCTGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGGCGC  >  NC_000913/3961876‑3962055
                                                                                          |                                                                                         
 cTGCATTCTTTCACCTGAATGACTTACTACAGTAAGATCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGCATGa                                                                                        >  2:2205807/1‑93 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       >  2:8955119/1‑93 (MQ=255)
    cATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAt                                                                                     <  1:8503997/93‑1 (MQ=255)
     aTTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtt                                                                                    >  2:9239995/1‑93 (MQ=255)
       tCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtttt                                                                                  <  1:914536/93‑1 (MQ=255)
        cTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAttttt                                                                                 >  1:4518664/1‑93 (MQ=255)
           tCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttat                                                                              >  2:4605042/1‑93 (MQ=255)
            cACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttatt                                                                             >  2:9817287/1‑93 (MQ=255)
             cccTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattt                                                                            <  2:8284780/92‑1 (MQ=255)
              ccTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattta                                                                           >  2:8504283/1‑93 (MQ=255)
              ccTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattta                                                                           >  1:3207087/1‑93 (MQ=255)
              ccTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattta                                                                           <  2:1586170/93‑1 (MQ=255)
               cTGAATCACTTAGTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTaa                                                                          <  1:3415082/93‑1 (MQ=255)
               cTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTaa                                                                          <  2:1939338/93‑1 (MQ=255)
               cTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTaa                                                                          <  2:9092534/93‑1 (MQ=255)
               cTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTaa                                                                          <  1:4037404/93‑1 (MQ=255)
                tGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAc                                                                         <  2:7792544/93‑1 (MQ=255)
                 gAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAcc                                                                        <  2:1194351/93‑1 (MQ=255)
                   aTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTg                                                                      <  1:5458713/93‑1 (MQ=255)
                    tCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGa                                                                     <  2:2432629/93‑1 (MQ=255)
                    tCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGa                                                                     >  1:8910712/1‑93 (MQ=255)
                     cACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGaa                                                                    <  1:2215424/93‑1 (MQ=255)
                     cACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGaa                                                                    <  1:8132155/93‑1 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  <  2:228036/93‑1 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  <  1:4295006/93‑1 (MQ=255)
                        ttACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacga                                                                 >  1:9235914/1‑93 (MQ=255)
                         tactacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgac                                                                <  1:8002031/93‑1 (MQ=255)
                         tactacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgac                                                                <  1:9641055/93‑1 (MQ=255)
                         tactacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgac                                                                <  2:4175319/93‑1 (MQ=255)
                         tactacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgac                                                                >  2:7234395/1‑93 (MQ=255)
                           ctacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgacga                                                              <  1:5683549/93‑1 (MQ=255)
                              cAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAAcgacgacga                                                           <  2:3509269/93‑1 (MQ=255)
                               aGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAt                                                          >  2:6528656/1‑93 (MQ=255)
                               aGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAt                                                          >  2:5587573/1‑93 (MQ=255)
                                gTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAtt                                                         >  1:5259192/1‑93 (MQ=255)
                                 gaaGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTa                                                        <  2:2039016/92‑1 (MQ=255)
                                  aaGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAc                                                       <  1:9883731/93‑1 (MQ=255)
                                  aaGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAc                                                       <  2:8087522/93‑1 (MQ=255)
                                   aGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTatg                                                      >  2:28188/1‑91 (MQ=255)
                                    gCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAcgc                                                     >  2:6093993/1‑93 (MQ=255)
                                       cATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGa                                                  <  2:706684/93‑1 (MQ=255)
                                       cATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGa                                                  <  1:4952656/93‑1 (MQ=255)
                                        aTCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAg                                                 >  1:4231505/1‑93 (MQ=255)
                                         tCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgc                                                >  1:3562928/1‑93 (MQ=255)
                                         tCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgc                                                <  2:5741851/93‑1 (MQ=255)
                                         tCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgc                                                >  1:7527994/1‑93 (MQ=255)
                                          cGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgcg                                               >  1:1040998/1‑93 (MQ=255)
                                           gggATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgcgc                                              >  1:1705066/1‑93 (MQ=255)
                                            ggATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGcc                                             <  2:5200214/93‑1 (MQ=255)
                                             gATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCa                                            <  1:9603280/93‑1 (MQ=255)
                                             gATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCa                                            <  1:6837563/93‑1 (MQ=255)
                                             gATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCa                                            <  2:7923773/93‑1 (MQ=255)
                                              attcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAg                                           >  1:7703172/1‑93 (MQ=255)
                                              attcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAg                                           >  2:4763826/1‑93 (MQ=255)
                                               ttcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAATGACGACGATTACGCGGAGCGCCAga                                          >  2:7317450/1‑92 (MQ=255)
                                               ttcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGt                                          >  2:7646217/1‑93 (MQ=255)
                                               ttcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGt                                          >  1:2519759/1‑93 (MQ=255)
                                                tcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTa                                         <  2:6928665/93‑1 (MQ=255)
                                                tcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTa                                         <  1:4629172/93‑1 (MQ=255)
                                                tcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTa                                         >  1:3672918/1‑93 (MQ=255)
                                                 cattcTCTTGCCATCTTGATACAGTTTGAATGGTTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTac                                        >  2:4019516/1‑93 (MQ=255)
                                                  attcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacg                                       <  1:2994015/93‑1 (MQ=255)
                                                  attcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacg                                       >  1:1465631/1‑93 (MQ=255)
                                                     ctctTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacgacg                                    <  1:1857587/93‑1 (MQ=255)
                                                      tctTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGAcgc                                   >  1:9432280/1‑93 (MQ=255)
                                                       ctTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGAcgcg                                  >  2:5450077/1‑93 (MQ=255)
                                                       ctTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGAcgcg                                  <  2:3895654/93‑1 (MQ=255)
                                                        ttGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCgg                                 >  1:864975/1‑93 (MQ=255)
                                                         tGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGa                                <  2:8790126/93‑1 (MQ=255)
                                                           ccATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGAcc                              <  2:8612968/93‑1 (MQ=255)
                                                             aTCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAt                            >  2:1050635/1‑93 (MQ=255)
                                                             aTCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAt                            <  1:2753215/93‑1 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           >  2:6370130/1‑93 (MQ=255)
                                                               cTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTg                          >  1:3765869/1‑93 (MQ=255)
                                                                ttGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGc                         >  2:8414458/1‑93 (MQ=255)
                                                                ttGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGc                         >  1:6284062/1‑93 (MQ=255)
                                                                ttGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGc                         >  1:7889980/1‑93 (MQ=255)
                                                                ttGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGc                         >  2:4294772/1‑93 (MQ=255)
                                                                  gATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCt                       >  2:6710595/1‑93 (MQ=255)
                                                                    tACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCtgt                     <  1:6741341/93‑1 (MQ=255)
                                                                    tACAGATTGAATGGTTTTGAGTGTGACGTTTTTTATTTAACCTCAACGACGACAATGTCGCGGCGCGCCAGTACGACGCGGACCATTAGCtgt                     >  2:7492966/1‑93 (MQ=255)
                                                                      cAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgc                   <  1:6950089/93‑1 (MQ=255)
                                                                      cAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgc                   >  1:4870734/1‑93 (MQ=255)
                                                                       aGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCGTTGCCTGTgcg                  <  2:5861766/93‑1 (MQ=255)
                                                                         tttGAATGATTTTGAGTGTCACATTTTTTATTTAGCCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGcgcg                <  1:646572/93‑1 (MQ=255)
                                                                           tGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCggg              >  1:7327701/1‑93 (MQ=255)
                                                                           tGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCggg              <  1:5938967/93‑1 (MQ=255)
                                                                           tGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCggg              <  2:3275266/93‑1 (MQ=255)
                                                                            gAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGgc             >  1:8250262/1‑93 (MQ=255)
                                                                             aaTGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGgcg            <  2:9524301/93‑1 (MQ=255)
                                                                               tGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGcgcg          >  2:2987743/1‑93 (MQ=255)
                                                                                 aTTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTg        <  2:3718145/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       <  1:1256571/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       <  2:9666135/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       <  1:2422083/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       >  2:2744793/1‑93 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       >  1:7140111/1‑93 (MQ=255)
                                                                                   tttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAg      <  2:256476/93‑1 (MQ=255)
                                                                                     tGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgc    >  1:6065945/1‑93 (MQ=255)
                                                                                       aGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgcgc  >  2:1763738/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCTGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGGCGC  >  NC_000913/3961876‑3962055

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