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,9660AG97.9% 758.9 / 1.6 238intergenic (‑122/+14)gppA/rhlBguanosine pentaphosphatase/putative ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (115/118);  total (117/121)
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.79e-01

Only 100 of 246 total aligned reads displayed.

GCTGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGGCGC  >  NC_000913/3961876‑3962055
                                                                                          |                                                                                         
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       >  2:1093704/1‑93 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       <  1:1562109/93‑1 (MQ=255)
  tGCATTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGAc                                                                                       >  1:3607742/1‑93 (MQ=255)
     aTTCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtt                                                                                    >  2:1600958/1‑93 (MQ=255)
      ttCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACAttt                                                                                   <  2:2494780/93‑1 (MQ=255)
       tCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtttt                                                                                  <  2:4078250/93‑1 (MQ=255)
       tCTTTCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtttt                                                                                  >  2:1787351/1‑93 (MQ=255)
         tttCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACAtttttt                                                                                >  1:4662558/1‑93 (MQ=255)
          ttCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTttta                                                                               <  2:3687999/93‑1 (MQ=255)
          ttCACCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTttta                                                                               <  2:648400/93‑1 (MQ=255)
             aCCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattt                                                                            <  2:2031813/93‑1 (MQ=255)
             aCCTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattt                                                                            >  2:614785/1‑93 (MQ=255)
              ccTGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTATTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTtttattta                                                                           >  2:582110/1‑93 (MQ=255)
               cttaaTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTaa                                                                          <  2:1225287/90‑1 (MQ=255)
               cTGAATCACTTACTACAGTAAGCACATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTaa                                                                          <  2:2583833/93‑1 (MQ=255)
                tGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAc                                                                         >  2:2446384/1‑93 (MQ=255)
                tGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAc                                                                         <  2:3132748/93‑1 (MQ=255)
                tGAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAc                                                                         <  2:443820/93‑1 (MQ=255)
                 gAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAcc                                                                        <  2:297994/93‑1 (MQ=255)
                 gAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAcc                                                                        <  2:2156377/93‑1 (MQ=255)
                 gAATCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAAcc                                                                        <  1:2806964/93‑1 (MQ=255)
                  aaTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCt                                                                       >  1:4235247/1‑93 (MQ=255)
                  aaTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCt                                                                       <  1:2420625/93‑1 (MQ=255)
                  aaTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCt                                                                       >  2:29066/1‑93 (MQ=255)
                  aaTCACTTACCACAGAAATCTCATCGGCATTCACGCTCTTGCCATCTTAAGACAGTTTGAAGGATTTTTAGTGTGACATTTTTTATTTAACCt                                                                       <  2:4254705/93‑1 (MQ=255)
                   aTCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTg                                                                      <  1:3844626/93‑1 (MQ=255)
                    tCACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGa                                                                     <  2:2409117/93‑1 (MQ=255)
                     cACTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGaa                                                                    <  2:3046358/93‑1 (MQ=255)
                      aCTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAac                                                                   <  1:661137/93‑1 (MQ=255)
                      aCTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAac                                                                   <  1:3828997/93‑1 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  <  1:238064/93‑1 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  >  1:917968/1‑93 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  <  2:130434/93‑1 (MQ=255)
                       cTTACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacg                                                                  >  2:3247887/1‑93 (MQ=255)
                        ttACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATGCAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacga                                                                 <  2:69968/93‑1 (MQ=255)
                        ttACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacga                                                                 <  1:1074325/93‑1 (MQ=255)
                        ttACTACAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacga                                                                 >  1:1417134/1‑93 (MQ=255)
                         tactacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgac                                                                >  1:2054388/1‑93 (MQ=255)
                         tactacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgac                                                                <  1:438327/93‑1 (MQ=255)
                          actacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgacg                                                               <  1:3958209/93‑1 (MQ=255)
                           ctacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgacga                                                              <  2:470440/93‑1 (MQ=255)
                           ctacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgacga                                                              <  1:128001/93‑1 (MQ=255)
                            tacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAacgacgac                                                             >  1:4278176/1‑93 (MQ=255)
                            tacAGTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGCgccacgac                                                             >  2:2298002/1‑93 (MQ=255)
                                gTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAtt                                                         >  1:2684797/1‑93 (MQ=255)
                                gTAAGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGAtt                                                         <  2:480673/93‑1 (MQ=255)
                                  aaGCTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAc                                                       <  1:1487503/93‑1 (MQ=255)
                                     cTCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTAcgcg                                                    <  1:1174074/93‑1 (MQ=255)
                                      tCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCgg                                                   >  1:756236/1‑93 (MQ=255)
                                      tCATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCgg                                                   >  2:2967668/1‑93 (MQ=255)
                                       cATCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGGTTTTGAGCGGGACATTTTTTATTTACCCTGAACGACGACGATTACGCGGa                                                  >  2:3502110/1‑93 (MQ=255)
                                       cATCGGGATTCACTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGa                                                  <  1:3176345/93‑1 (MQ=255)
                                         tCGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgc                                                >  1:2076399/1‑93 (MQ=255)
                                          cGGGATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAgcg                                               >  1:1431827/1‑93 (MQ=255)
                                            ggATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGcc                                             <  1:1872070/93‑1 (MQ=255)
                                            ggATTCATTCTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGcc                                             >  2:3602839/1‑93 (MQ=255)
                                              attcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAg                                           >  1:839110/1‑93 (MQ=255)
                                              attcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAg                                           >  1:4036734/1‑93 (MQ=255)
                                              attcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAg                                           >  2:1522326/1‑93 (MQ=255)
                                              attcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAg                                           <  1:3649410/93‑1 (MQ=255)
                                               ttcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGt                                          >  2:2240231/1‑93 (MQ=255)
                                               ttcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGt                                          >  2:1276947/1‑93 (MQ=255)
                                               ttcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGt                                          >  1:4370114/1‑93 (MQ=255)
                                                tcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTa                                         <  2:4199309/93‑1 (MQ=255)
                                                tcattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTa                                         <  1:1003606/93‑1 (MQ=255)
                                                 tattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTac                                        <  1:1590369/92‑1 (MQ=255)
                                                 cattcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACTACGATTACGCGGAGCGCCAGTac                                        >  2:1824207/1‑93 (MQ=255)
                                                  attcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGGCAGTacg                                       >  1:3770589/1‑93 (MQ=255)
                                                  attcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacg                                       >  2:1463856/1‑93 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      <  2:3784463/93‑1 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      >  1:3000221/1‑93 (MQ=255)
                                                   ttcTCTTGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTacga                                      >  2:1339809/1‑93 (MQ=255)
                                                        ttGCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCgg                                 >  2:839283/1‑93 (MQ=255)
                                                          gCCATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGAc                               >  2:4647390/1‑93 (MQ=255)
                                                            cATCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGAcaa                             >  2:268927/1‑91 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           <  1:1897716/93‑1 (MQ=255)
                                                              tCTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCAtt                           <  1:3201946/93‑1 (MQ=255)
                                                               cTTGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTg                          >  2:4348822/1‑93 (MQ=255)
                                                                ttGATACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGc                         >  1:2203781/1‑93 (MQ=255)
                                                                    tACAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACAACGCGGACAATTGCCtgt                     >  2:17007/1‑93 (MQ=255)
                                                                    tACAGTTTGAATGATTTTGAGTATGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCtgt                     >  1:2345116/1‑93 (MQ=255)
                                                                     aCAGTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCtgtg                    >  1:536667/1‑93 (MQ=255)
                                                                        gTTTGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTgcgc                 >  1:3282581/1‑93 (MQ=255)
                                                                         tttGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGcgcg                >  1:1751726/1‑93 (MQ=255)
                                                                         tttGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGcgcg                >  2:1415919/1‑93 (MQ=255)
                                                                          ttGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCgg               >  2:3760714/1‑93 (MQ=255)
                                                                           tGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCggg              <  2:3317016/93‑1 (MQ=255)
                                                                           tGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCggg              >  2:1191440/1‑93 (MQ=255)
                                                                           tGAATGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCggg              >  2:4666045/1‑93 (MQ=255)
                                                                           tGAATGATGTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCggg              <  1:2286200/93‑1 (MQ=255)
                                                                              aTGATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGgcgc           <  1:899495/93‑1 (MQ=255)
                                                                                gATTTTGAGTGTGACATTTTTTATTTAACCTGAGCGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGt         <  1:4564120/93‑1 (MQ=255)
                                                                                gATTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGt         <  1:1679509/93‑1 (MQ=255)
                                                                                 aTTTTGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTg        >  1:4355224/1‑93 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTgg       >  1:3358055/1‑92 (MQ=255)
                                                                                  ttttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCGGTGCGCGGGCGGGTgg       >  2:1640440/1‑92 (MQ=255)
                                                                                   tttGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAg      <  1:485831/93‑1 (MQ=255)
                                                                                     tGAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgc    <  1:292378/93‑1 (MQ=255)
                                                                                      gAGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgcg   >  1:687412/1‑93 (MQ=255)
                                                                                       aGTGTGACATTTTTTATTTAACCTGAACGACGACGATTACGCGGAGCGCCAGTACGACGCGGACCATTGCCTGTGCGCGGGCGCGTGAGgcgc  <  2:9198/93‑1 (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 < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: