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.2% 1565.9 / ‑1.2 493intergenic (‑122/+14)gppA/rhlBguanosine pentaphosphatase/putative ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (250/239);  total (252/241)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 505 total aligned reads displayed.

GCTGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGGCGC  >  NC_000913/3961876‑3962055
                                                                                          |                                                                                         
gCTGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGtgtg                                                                                         >  2:2570193/1‑93 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       <  1:2382822/93‑1 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       <  2:6318474/93‑1 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       >  2:62512/1‑93 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       <  2:5884888/93‑1 (MQ=255)
   gCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACa                                                                                      <  1:4170521/93‑1 (MQ=255)
   gCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACa                                                                                      >  1:3376338/1‑93 (MQ=255)
    cATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAt                                                                                     >  2:2238334/1‑93 (MQ=255)
      ttCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAttt                                                                                   >  1:7317584/1‑93 (MQ=255)
       tCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtttt                                                                                  >  1:3894343/1‑93 (MQ=255)
        cTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAttttt                                                                                 >  2:763957/1‑93 (MQ=255)
          ttCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTttta                                                                               <  1:6066181/93‑1 (MQ=255)
            cACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttatt                                                                             <  2:1454868/93‑1 (MQ=255)
            cACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttatt                                                                             <  1:3186929/93‑1 (MQ=255)
             aCCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattt                                                                            <  2:1815976/93‑1 (MQ=255)
              ccTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattta                                                                           >  2:5757496/1‑93 (MQ=255)
                tGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAc                                                                         >  1:4006866/1‑93 (MQ=255)
                 gAATCATTTTCTACAGTTTTCTCATTGGGTGTCATTCTCTTTTTATCTTGATACTGTTTGTAGGATTTCGTGTGTGACATTTTTTATTTTAcc                                                                        >  2:6005261/1‑93 (MQ=255)
                 gAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAcc                                                                        <  1:4700801/93‑1 (MQ=255)
                 gAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAcc                                                                        >  2:3747397/1‑93 (MQ=255)
                 gAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAcc                                                                        >  2:1377823/1‑93 (MQ=255)
                  aaTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCt                                                                       <  2:1102819/93‑1 (MQ=255)
                   aTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTg                                                                      <  1:5489948/93‑1 (MQ=255)
                     cACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGaa                                                                    <  2:5373682/93‑1 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  >  2:7839705/1‑93 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  >  2:549163/1‑93 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  >  1:6431855/1‑93 (MQ=255)
                        ttACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacga                                                                 >  1:1450975/1‑93 (MQ=255)
                         tactacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgac                                                                <  2:551118/93‑1 (MQ=255)
                           ctacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgacga                                                              <  1:1206826/93‑1 (MQ=255)
                           ctacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgacga                                                              <  1:588408/93‑1 (MQ=255)
                              cAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAAcgacgacga                                                           <  1:177763/93‑1 (MQ=255)
                                gTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAtt                                                         >  2:1030979/1‑93 (MQ=255)
                                gTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAtt                                                         <  1:6302318/93‑1 (MQ=255)
                                gTAAGCTCAGCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTGAACCTGAACGACGTCGAtt                                                         >  2:7037150/1‑93 (MQ=255)
                                 tAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAGCGACGACGATTa                                                        <  1:5197557/93‑1 (MQ=255)
                                   aGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAcg                                                      <  2:2735978/93‑1 (MQ=255)
                                    gCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAcgc                                                     <  1:2020068/93‑1 (MQ=255)
                                    gCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAcgc                                                     >  1:4403070/1‑93 (MQ=255)
                                      tCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGCCGATTACGCgg                                                   >  1:1834094/1‑93 (MQ=255)
                                      tCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCgg                                                   <  2:2154486/93‑1 (MQ=255)
                                      tCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCgg                                                   >  2:4579663/1‑93 (MQ=255)
                                       cATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGa                                                  <  1:5408940/93‑1 (MQ=255)
                                       cATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGa                                                  <  1:567672/93‑1 (MQ=255)
                                        aTCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAg                                                 >  1:7192981/1‑93 (MQ=255)
                                         tCGGTATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgc                                                <  2:7521099/93‑1 (MQ=255)
                                         tCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGTACGACGACGATTACGCGGAgc                                                >  2:4940946/1‑93 (MQ=255)
                                         tCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgc                                                <  1:524525/93‑1 (MQ=255)
                                         tCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgc                                                >  2:1627129/1‑93 (MQ=255)
                                          cGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACAATTACGCGAAgca                                               >  1:2090423/1‑92 (MQ=255)
                                           gggATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgcgc                                              <  2:4262629/93‑1 (MQ=255)
                                             gATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCa                                            >  1:1594480/1‑93 (MQ=255)
                                              attcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAg                                           <  2:6189765/93‑1 (MQ=255)
                                               ttcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGt                                          >  1:6563888/1‑93 (MQ=255)
                                               ttcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGt                                          >  1:4567066/1‑93 (MQ=255)
                                               ttcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGt                                          >  1:248739/1‑93 (MQ=255)
                                                 cattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTac                                        >  2:4003853/1‑93 (MQ=255)
                                                 cattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTac                                        <  1:6698282/93‑1 (MQ=255)
                                                  attcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacg                                       >  2:4132855/1‑93 (MQ=255)
                                                  attcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacg                                       >  2:6546904/1‑93 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGCTTTTGAGTGTGACATTTTTTATTTAACCTGAACGATGACGATTACGCGGAGCGCCAGTacga                                      >  1:2603014/1‑93 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      >  2:6746848/1‑93 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      >  1:1023464/1‑93 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      >  2:836038/1‑93 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      <  1:5056155/93‑1 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      >  1:2287805/1‑93 (MQ=255)
                                                    tctcTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacgac                                     >  1:5584872/1‑93 (MQ=255)
                                                    tctcTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacgac                                     <  1:4275203/93‑1 (MQ=255)
                                                      tctTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGAcgc                                   >  2:3370691/1‑93 (MQ=255)
                                                           ccATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGAcc                              >  2:5615438/1‑93 (MQ=255)
                                                           ccATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGAcc                              >  2:712596/1‑93 (MQ=255)
                                                             aTCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAt                            >  1:1321025/1‑93 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           >  1:2852620/1‑93 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           >  2:4712449/1‑93 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           >  1:2974074/1‑93 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           >  1:4828497/1‑93 (MQ=255)
                                                                ttGATACAGTTTTAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGc                         >  2:3307129/1‑93 (MQ=255)
                                                                ttGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGc                         >  2:1995044/1‑93 (MQ=255)
                                                                   aTACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCtg                      <  2:6863980/93‑1 (MQ=255)
                                                                   aTACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCtg                      >  2:3562174/1‑93 (MQ=255)
                                                                    tACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCtgt                     <  2:4427300/93‑1 (MQ=255)
                                                                     aCAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCtgtg                    <  1:3540803/93‑1 (MQ=255)
                                                                      cAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgc                   >  2:3844730/1‑93 (MQ=255)
                                                                      cAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgc                   >  1:1936463/1‑93 (MQ=255)
                                                                        gTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgcgc                 >  1:6895725/1‑93 (MQ=255)
                                                                        atttGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgcgc                 <  2:7350935/92‑1 (MQ=255)
                                                                         tttGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGcgcg                >  1:7744472/1‑93 (MQ=255)
                                                                         tttGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGcgcg                >  1:823178/1‑93 (MQ=255)
                                                                         tttGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGcgcg                >  2:1297577/1‑93 (MQ=255)
                                                                           tGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCggg              <  1:6733872/93‑1 (MQ=255)
                                                                            gAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGgc             >  1:7700671/1‑93 (MQ=255)
                                                                            gAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGgc             >  2:2875052/1‑93 (MQ=255)
                                                                             aaTGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGgcg            >  2:5203596/1‑93 (MQ=255)
                                                                              aTGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGgcgc           <  2:2342028/93‑1 (MQ=255)
                                                                               tGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGcgcg          >  1:582515/1‑93 (MQ=255)
                                                                                gATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGt         <  2:15513/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       <  1:2709289/93‑1 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGa       >  2:3049206/1‑93 (MQ=255)
                                                                                     tGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgc    <  1:758600/93‑1 (MQ=255)
                                                                                      gAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCGTGGGCGCGGGCGCGTGAGgcg   >  1:6979136/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCTGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGGCGC  >  NC_000913/3961876‑3962055

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