Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 41.3% intergenic (‑194/+114) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG41.3% 281.1 / inf 242intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (70/72);  new base G (89/11);  total (159/83)
Fisher's exact test for biased strand distribution p-value = 3.94e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-02

Only 100 of 245 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                         
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         <  1:1709789/100‑1 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGCCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                       >  2:519938/1‑100 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCTGCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCTGATAAGGTGTTAACGCCGCATCCgg                                                                                                    <  2:1614414/99‑1 (MQ=25)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                 <  2:1683930/100‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                 <  2:743460/100‑3 (MQ=37)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑a                                                                                              <  1:2517435/100‑1 (MQ=25)
            ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑Ac                                                                                             >  1:1625871/1‑100 (MQ=255)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑a                                                                                           <  2:1529982/100‑1 (MQ=255)
                  gCACGCAGAGGCACAGCGCCATCGGGTAGACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                                         <  2:2831839/100‑11 (MQ=18)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                                     <  2:2057814/100‑1 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                                     <  1:2295274/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                   >  2:613041/1‑100 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                   <  1:1370036/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                  <  2:3582598/100‑1 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCa                                                                                 <  1:206555/100‑1 (MQ=255)
                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                            <  2:3524788/100‑1 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                            <  1:3035695/99‑1 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCa                                                                           >  1:3140673/1‑100 (MQ=255)
                              gcaGCGCCATCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCaa                                                                          >  1:3643810/2‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                         <  2:310388/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                         <  2:671520/100‑1 (MQ=255)
                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCTGCCT‑‑AC‑AAAATCGTGCAAATTCAATa                                                                        <  1:2219017/100‑1 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATAc                                                                       >  1:2019660/1‑100 (MQ=255)
                                  cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACa                                                                      <  1:213036/100‑1 (MQ=255)
                                   gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAt                                                                     >  1:134892/1‑100 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                    >  1:950880/1‑100 (MQ=255)
                                     cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTg                                                                   <  1:3361132/100‑1 (MQ=255)
                                             ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t‑t                                                        >  1:2943577/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCg                                                     <  2:647631/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCg                                                     <  1:3068839/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCg                                                     >  1:1918879/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCg                                                     <  1:1469346/100‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTa                                                   >  2:3064913/1‑100 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTa                                                   >  2:3407641/1‑100 (MQ=255)
                                                    gcataaaCAAAGCGCACTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAgg                                                 >  2:1160602/7‑100 (MQ=18)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGcc                                               <  1:3311824/100‑1 (MQ=255)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGcc                                               <  1:2983164/100‑1 (MQ=255)
                                                      ataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGcc                                               >  1:2819664/5‑100 (MQ=18)
                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCt                                              <  2:1056185/100‑1 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTg                                             >  1:355915/1‑100 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGc                                       <  1:3213095/100‑1 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGAA‑GAAT‑T‑TTCGTGGGACTGATatg                                        >  2:3579623/1‑97 (MQ=17)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCg                                      >  1:115850/1‑100 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCg                                      <  2:1865346/100‑1 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCgg                                     <  1:2445447/100‑1 (MQ=255)
                                                                 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGa                                    <  1:805363/100‑1 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGCCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGACCTGATAAGCgg                                     <  2:2551137/98‑1 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAgcg                                 >  1:446404/1‑100 (MQ=255)
                                                                    agttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAAAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑A‑CACGTAGGCCTGATAAGCGAAg                                   >  2:754061/11‑100 (MQ=2)
                                                                     tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAgcgc                                >  2:328768/1‑100 (MQ=255)
                                                                      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCa                               >  1:1578042/1‑100 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                            >  1:1692890/1‑100 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑CTTGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  1:2192925/5‑100 (MQ=18)
                                                                          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑ATCT‑T‑CTGGTAGGCCTGATAAGCATCGCGCATCAg                           >  2:2763245/5‑100 (MQ=18)
                                                                           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑CTTGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  1:3458623/4‑100 (MQ=18)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑TTGGTAGGCCTGATAAGCGCAGCGCATCAgg                          >  1:2786673/4‑100 (MQ=21)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑AC‑AAATTCTTCCAAATTCAATATATTGCA‑GAAA‑T‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:1822416/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAACATATTGCA‑‑ATTCT‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:861865/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑ATCT‑T‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGc                         >  2:488448/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCAGGGAC‑T‑AT‑GTAGGCCTGATAAGCATAGCGCATCAGGc                         >  1:1032579/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCAGGGAC‑T‑AT‑GTAGGCCTGATAAGCATAGCGCATCAGGc                         >  1:549715/3‑100 (MQ=18)
                                                                             ggtcggataaggcgTTCACGCCGCATCCGACCC‑‑AC‑AAAGTCTTGCAAAATCAATAAATTGCA‑GGAA‑C‑GATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2149564/11‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  1:1847174/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2437847/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AGAT‑T‑TTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:145228/2‑100 (MQ=17)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:1986198/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑C‑CACGTAGGCCTGATAAGCGTAGCGCATCAgaga                        >  1:287876/1‑96 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑T‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                      >  2:2222745/1‑100 (MQ=18)
                                                                               ccggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑ATCT‑T‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGCaag                      >  1:1770171/9‑98 (MQ=18)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:2357918/9‑98 (MQ=17)
                                                                                cggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                    >  1:1821163/8‑99 (MQ=21)
                                                                                cggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑CACGTATGCCTGATAAGCGTAGAGCATCAGGCa                        >  1:1318405/8‑97 (MQ=17)
                                                                                cggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                     >  1:1945411/8‑100 (MQ=17)
                                                                                cGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                     <  2:2670477/100‑1 (MQ=255)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAAT‑CCCTTGTAGGCCTGATAAGCGTAGCGCATCAGCCAtttgc                   >  2:694816/7‑98 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1127222/7‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑ATCT‑T‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGCaagct                    >  1:997406/7‑96 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1726801/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2485018/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGAAgttt                    >  1:3546941/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AGAT‑T‑TTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                    >  1:1778313/7‑96 (MQ=17)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                   >  1:2710035/6‑100 (MQ=17)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AGAT‑T‑TTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacgctc                   >  1:3251252/6‑95 (MQ=17)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑C‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttg                   >  1:123204/6‑99 (MQ=18)
                                                                                  gataaggcgTTAACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                   >  1:1550911/6‑100 (MQ=18)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:1065965/100‑1 (MQ=255)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:2571908/100‑1 (MQ=255)
                                                                                     gcggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcn                <  2:135915/100‑2 (MQ=255)
                                                                                     gcggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  2:1505022/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:3140650/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         <  1:3407666/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  1:620544/1‑100 (MQ=255)
                                                                                             aaCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑C‑TTCGTAGGCCTGATAAGCGGATCGCATCAGGCATTtttgcctttgtcct        >  2:2211785/1‑96 (MQ=255)
                                                                                              aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  1:672270/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:1957047/100‑1 (MQ=255)
                                                                                               cGCCTTATCCGGCCT‑‑AC‑AAGATCGTGAAAATTCAATACATTGCA‑GCAT‑T‑TTCGTAGGCCGGATCAGCGGAACGCATTAGGCATTtttgcctttggta         >  2:3022143/1‑94 (MQ=17)
                                                                                               cGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  2:3341122/1‑99 (MQ=255)
                                                                                                gCCTTATCCGGCGT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:1266239/100‑1 (MQ=255)
                                                                                                   ttATCCGGCTT‑‑AC‑AAAATCGTGCAAATTCAATAAATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCAATCAGTCt  <  2:1240889/100‑1 (MQ=255)
                                                                                                   ttATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  1:1135717/100‑1 (MQ=255)
                                                                                                   |                                                                                                         
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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