Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,407,859 A→G 100% intergenic (‑64/‑191) ydaM ← / → ydaN diguanylate cyclase, csgD regulator/putative Zn(II) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,407,8590AG98.3% 381.4 / ‑3.3 115intergenic (‑64/‑191)ydaM/ydaNdiguanylate cyclase, csgD regulator/putative Zn(II) transporter
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (62/51);  total (63/52)
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.01e-01

Only 100 of 115 total aligned reads displayed.

GTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGA  >  NC_000913/1407762‑1407952
                                                                                                 |                                                                                             
gTAAGTCCAGGGTATTTAAGTTGTGCGTAATCATTGAGATCCCGTTTTTTTTATTAATTTTTGTTGCTATCCGTTCGCAAAAAAGCCACGTCATATGGTc                                                                                             <  2:1146039/100‑1 (MQ=255)
 tAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCa                                                                                            <  1:2305652/100‑1 (MQ=255)
   aGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGa                                                                                          >  1:820334/1‑100 (MQ=255)
    gTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGAt                                                                                         >  1:2863425/1‑100 (MQ=255)
    gTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGAt                                                                                         <  1:1726779/100‑1 (MQ=255)
     tCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGa                                                                                          >  2:2270853/1‑98 (MQ=255)
       cAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAAt                                                                                      <  1:484428/100‑1 (MQ=255)
       cAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAAt                                                                                      >  2:2354219/1‑100 (MQ=255)
       cAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAAt                                                                                      >  2:1124527/1‑100 (MQ=255)
        ngggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAAt                                                                                      >  2:1407738/2‑99 (MQ=255)
         gggTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATAAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTc                                                                                    <  1:996797/100‑1 (MQ=255)
          ggTATTGAAGTTGTGCGTAATCATTGAGTTCTCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCt                                                                                   <  1:894779/100‑1 (MQ=255)
           gTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCt                                                                                   <  2:131157/99‑1 (MQ=255)
           gTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTg                                                                                  <  1:2722711/100‑1 (MQ=255)
           gTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTg                                                                                  <  1:1246038/100‑1 (MQ=255)
           gTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTg                                                                                  <  2:1403188/100‑1 (MQ=255)
             aTTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGAt                                                                                <  1:458052/100‑1 (MQ=255)
             aTTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGAt                                                                                >  1:953108/1‑100 (MQ=255)
                gAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAAt                                                                             >  2:2669896/1‑100 (MQ=255)
                gAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAAt                                                                             >  2:225438/1‑100 (MQ=255)
                  aGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGa                                                                           <  2:2853612/100‑1 (MQ=255)
                     tgtgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAg                                                                        <  2:629951/100‑1 (MQ=255)
                     tgtgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAg                                                                        >  2:1864369/1‑100 (MQ=255)
                       tgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGAc                                                                      >  2:2476693/1‑100 (MQ=255)
                           tAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTa                                                                  >  1:2159917/1‑100 (MQ=255)
                           tAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTa                                                                  >  1:178671/1‑100 (MQ=255)
                           tAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTa                                                                  >  2:605995/1‑100 (MQ=255)
                           naaTCATTGAGATCCCGTTATTTTTTTTTTTTTTTTTTGATATCCGATCGCAAAAAATCTATGTTATATGGTCAGATAATTCTTATATTGATTGAc                                                                      >  2:1310882/2‑96 (MQ=255)
                              tCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTAttt                                                               <  2:431535/100‑1 (MQ=255)
                              tCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTAttt                                                               <  2:1688142/100‑1 (MQ=255)
                              tCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTAttt                                                               >  2:43913/1‑100 (MQ=255)
                               cATTGAGATCCCGTTATTTTTATCAATTTATGTTGCTATCCGATAGCAAAAAAGCCAAGTCATATAGTCAGATAATTCTGATAAAGATAGACGCTATTTa                                                              >  2:2267096/1‑100 (MQ=255)
                                aTTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTaa                                                             >  2:1684433/1‑100 (MQ=255)
                                     gaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACtt                                                        >  1:537671/1‑100 (MQ=255)
                                     gaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACtt                                                        >  1:2611741/1‑100 (MQ=255)
                                      aTCCTGTTATTTTTATCAATTTTTGTTTCTCTCCGATCGCAAAAAAGCCACGTCATATGGTCATATAATTCTGATATTTATAGTCGCTATTTAACACTTc                                                       <  2:1924931/100‑1 (MQ=255)
                                      aTCCCTTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTc                                                       <  2:379849/100‑1 (MQ=255)
                                       cttccgatcttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTca                                                      <  2:814321/91‑1 (MQ=255)
                                        cccGTTATTTTTATCAATTTTTTTTTCTTTTCTATCGCAAAAAAGCAACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTcac                                                     >  2:1674275/1‑100 (MQ=255)
                                         ccGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTcaca                                                    >  2:2769089/1‑100 (MQ=255)
                                           gTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCacac                                                   <  1:2841523/99‑1 (MQ=255)
                                           gTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACg                                                  >  1:2061614/1‑100 (MQ=255)
                                             tATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGt                                                >  1:1320909/1‑100 (MQ=255)
                                             tATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGt                                                >  1:2156910/1‑100 (MQ=255)
                                              aTTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGtt                                               >  1:794797/1‑100 (MQ=255)
                                               tttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGttt                                              >  1:1613944/1‑100 (MQ=255)
                                                ttttATCAATTTTTTTTTTTATTCGTTCGTAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTg                                             <  1:1028537/100‑1 (MQ=255)
                                                ttttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTg                                             >  2:1924364/1‑100 (MQ=255)
                                                 tttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCATGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGt                                            <  2:2276903/100‑1 (MQ=255)
                                                 tttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTcaagatcggaag                                            >  1:814325/1‑90 (MQ=255)
                                                 tttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGt                                            >  1:453907/1‑100 (MQ=255)
                                                  ttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGATGCTATTTAACACTTCGCACGGTTTgga                                           <  2:1588004/100‑3 (MQ=255)
                                                  ttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGta                                           <  1:1600745/100‑1 (MQ=255)
                                                  ttATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGta                                           <  1:1272744/100‑1 (MQ=255)
                                                   tATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTAACACGGTTTGtat                                          >  2:876298/1‑100 (MQ=255)
                                                    aTCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGtata                                         >  2:1607193/1‑100 (MQ=255)
                                                     tCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATAc                                        <  1:1060636/100‑1 (MQ=255)
                                                      cAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATAc                                        >  1:354386/1‑99 (MQ=255)
                                                      cAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACg                                       >  2:2518259/1‑100 (MQ=255)
                                                       aaTTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACgg                                      <  2:453904/100‑1 (MQ=255)
                                                       aaTTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACgg                                      >  1:2662106/1‑100 (MQ=255)
                                                       aaTTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACgg                                      >  1:971155/1‑100 (MQ=255)
                                                       aaTTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACgg                                      <  1:2427154/100‑1 (MQ=255)
                                                           tttGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGaaa                                   <  2:1244129/99‑1 (MQ=255)
                                                            ttgttgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAg                                 <  2:1480830/100‑1 (MQ=255)
                                                            ttgttgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAg                                 <  2:1988739/100‑1 (MQ=255)
                                                              gttgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCa                               >  2:2008205/1‑100 (MQ=255)
                                                               ttgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAt                              >  1:921019/1‑100 (MQ=255)
                                                                tgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATACTGATAGACGCTATTTACAACTTCACACGGTTTGTATACGGAAAAGCAt                              >  2:2787106/1‑99 (MQ=255)
                                                                tgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtt                             >  2:3018744/1‑100 (MQ=255)
                                                                tgCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtt                             >  1:1112622/1‑100 (MQ=255)
                                                                ctcTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtt                             <  2:2008810/98‑1 (MQ=255)
                                                                 gCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                            <  2:1730611/100‑1 (MQ=255)
                                                                 gCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                            <  2:1112612/100‑1 (MQ=255)
                                                                 gCTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                            >  1:1036738/1‑100 (MQ=255)
                                                                  cTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtttt                           >  1:2314740/1‑100 (MQ=255)
                                                                  cTATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttag                           >  1:2008833/1‑98 (MQ=255)
                                                                   tATCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTg                          <  2:2176647/100‑1 (MQ=255)
                                                                    agctgATAGCAAAAAAGCTAAGTCAAATTGTCAGATAATTCTGATAATGATATACGCTCTTTAGCACTTCTCACGGTTTGTATACGGAAAAGCAttta                           <  2:261923/94‑2 (MQ=255)
                                                                     tCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCt                        >  1:2143129/1‑100 (MQ=255)
                                                                     tCCGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCt                        >  1:845552/1‑100 (MQ=255)
                                                                      ccGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCtt                       <  1:2496143/100‑1 (MQ=255)
                                                                      ccGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCtt                       <  1:1877996/100‑1 (MQ=255)
                                                                      ccGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCtn                       <  1:1877867/100‑2 (MQ=255)
                                                                       cGATCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCttt                      >  2:1910772/1‑100 (MQ=255)
                                                                          tCGCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACAGAAAAGCAAAATGCttttc                    <  2:86602/99‑2 (MQ=255)
                                                                           cGCAAAAAAGCCCCGTCATATGGTCAGATAATTCTGATAATGATATACGCTATTTAACACTTCACACGTTTTGTATACGGAAAAGCAAAATGCttttccg                  <  2:1503611/100‑4 (MQ=255)
                                                                            gCAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTa                 >  2:209193/1‑100 (MQ=255)
                                                                             cAAAAAAGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTAt                >  2:487544/1‑100 (MQ=255)
                                                                              aaaaaaGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAAAATGCttttccg                  <  1:227301/97‑4 (MQ=255)
                                                                                  aaGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAAt           <  2:1720557/100‑1 (MQ=255)
                                                                                   aGCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAAtt          <  2:1968217/100‑1 (MQ=255)
                                                                                    gCCACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAAttt         <  1:2493864/100‑1 (MQ=255)
                                                                                      cACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAg       <  1:2911254/100‑1 (MQ=255)
                                                                                      cACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAg       <  1:1779179/100‑1 (MQ=255)
                                                                                      cACGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTATGCTTTTTGTATTCAATTTAg       >  1:2558027/1‑100 (MQ=255)
                                                                                       aCGTCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGa      <  1:2076596/100‑1 (MQ=255)
                                                                                          tCATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACa    >  1:611322/1‑99 (MQ=255)
                                                                                           cATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTCTACGGACAAGCATTTTGCTTTTTGTATTCAATTTAGACAg   >  2:1433684/1‑99 (MQ=255)
                                                                                           cATATGGTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGa  >  2:385582/1‑100 (MQ=255)
                                                                                                 |                                                                                             
GTAAGTCCAGGGTATTGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGA  >  NC_000913/1407762‑1407952

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