Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,706,0170CT21.1% 215.4 / 46.9 94intergenic (‑227/+81)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (40/31);  new base T (19/0);  total (63/31)
Fisher's exact test for biased strand distribution p-value = 2.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

GCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTT  >  NC_000913/3705920‑3706115
                                                                                                  |                                                                                                     
gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                                                     >  1:1498831/1‑100 (MQ=255)
 ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                                                    <  2:609643/100‑1 (MQ=255)
 ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                                                    <  2:609744/100‑1 (MQ=255)
  cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                                                   >  1:605013/1‑100 (MQ=255)
  cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                                                   <  2:1505975/100‑1 (MQ=255)
   aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                                                  <  1:228258/100‑1 (MQ=255)
    ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                                                 <  2:490392/100‑1 (MQ=255)
    ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                                                 >  2:1094653/1‑100 (MQ=255)
      gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                                                              >  1:300835/1‑100 (MQ=255)
            aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑Tc                                                                                      <  2:1647620/100‑1 (MQ=255)
            aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑Tc                                                                                      >  2:482641/1‑100 (MQ=255)
               gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTa                                                                                   <  2:1371717/100‑1 (MQ=255)
                 ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAgg                                                                                 >  2:1250555/1‑100 (MQ=255)
                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGcc                                                                               >  1:1096353/1‑100 (MQ=255)
                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGcc                                                                               <  2:1249185/100‑1 (MQ=255)
                    ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCt                                                                              >  1:1627216/1‑100 (MQ=255)
                    ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCt                                                                              <  2:1031849/100‑1 (MQ=255)
                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTg                                                                             >  2:198267/1‑100 (MQ=255)
                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTg                                                                             >  2:152690/1‑100 (MQ=255)
                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATCCATTGCA‑GAAT‑TT‑TCGTAGGCCTGa                                                                            >  1:1352443/1‑100 (MQ=255)
                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGa                                                                            >  1:425821/1‑100 (MQ=255)
                       aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGAt                                                                           <  1:1317016/100‑1 (MQ=255)
                          aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAg                                                                        >  1:795600/1‑100 (MQ=255)
                             tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCgg                                                                     >  1:1238273/1‑100 (MQ=255)
                              gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGa                                                                    >  1:554668/1‑100 (MQ=255)
                               cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAg                                                                   >  1:865909/1‑100 (MQ=255)
                               cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAg                                                                   >  2:193449/1‑100 (MQ=255)
                                  gatttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TT‑ATGTAGGCCTGATAAGCGTAGCGCa                                                               >  2:1115534/10‑100 (MQ=18)
                                   ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCa                                                               >  1:423905/1‑100 (MQ=255)
                                   atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCC‑TTGTAGGCCTGATAAGCGTAGCGCAt                                                              >  2:904568/9‑100 (MQ=18)
                                    tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCAt                                                              >  2:1071092/1‑100 (MQ=255)
                                     gTTTATGCCGGATGCGGCGTTAACGCCTGATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATc                                                             <  2:71491/100‑1 (MQ=255)
                                       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAg                                                           <  2:298353/100‑1 (MQ=255)
                                       gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CA‑T‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                           >  2:1508229/5‑100 (MQ=18)
                                       gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAg                                                           >  2:244991/5‑100 (MQ=18)
                                         aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAACATATTGCA‑ATTC‑TC‑TTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                          >  1:1699646/3‑99 (MQ=18)
                                         aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAACATATTGCA‑‑ATT‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                         >  1:1506162/3‑100 (MQ=18)
                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CA‑T‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                        >  1:767661/2‑100 (MQ=18)
                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑ATA‑TG‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGcga                                                       >  2:953868/2‑98 (MQ=21)
                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑AT‑TGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                        >  1:1085379/2‑100 (MQ=21)
                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑AT‑TGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                        >  2:131046/2‑100 (MQ=21)
                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑CC‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                        >  2:288038/2‑100 (MQ=18)
                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CC‑TCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                        >  2:351287/2‑100 (MQ=21)
                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CT‑GCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                        >  2:1338067/2‑100 (MQ=17)
                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCC‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                      >  1:1455779/1‑100 (MQ=18)
                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GGGA‑CT‑ATGTAGGCCTGATAAGCATAGCGCATCAGGCag                                                       >  2:1235653/1‑99 (MQ=18)
                                            ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGGAC‑TA‑T‑GTAGGCCTGATAAGCATAGCGCATCAGGCagc                                                      >  2:711970/9‑98 (MQ=18)
                                            ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                      >  2:1198727/1‑100 (MQ=255)
                                            ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                      >  2:1528506/1‑100 (MQ=255)
                                             cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑ATCT‑TC‑TGGTAGGCCTGATAAGCATCGCGCATCAGGCaagc                                                     >  2:413237/8‑97 (MQ=18)
                                             cggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TT‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCAttta                                                    >  1:114986/8‑99 (MQ=21)
                                              ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑AT‑TGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                    >  2:681480/7‑96 (MQ=25)
                                              ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CC‑TCGTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                                                    >  2:36398/7‑96 (MQ=21)
                                              ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑AC‑ACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                                                    >  1:1510346/7‑100 (MQ=18)
                                              ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑AC‑ACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                                                    >  2:191830/7‑100 (MQ=18)
                                               gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CC‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                   >  2:1353041/6‑100 (MQ=18)
                                               gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CC‑TCGTAGGCCTGATAAGCGTAGCGCATCAGGCaattct                                                   >  2:1032851/6‑95 (MQ=25)
                                                aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                  <  2:1041634/100‑1 (MQ=255)
                                                  gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                <  2:1130083/100‑1 (MQ=255)
                                                   cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                               >  1:249202/1‑100 (MQ=255)
                                                   cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                               >  1:390453/1‑100 (MQ=255)
                                                        ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCACTACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                            <  1:1081175/98‑1 (MQ=255)
                                                         tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                         >  2:691652/1‑100 (MQ=255)
                                                          aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                        >  1:901263/1‑100 (MQ=255)
                                                           aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                       <  2:255049/100‑1 (MQ=255)
                                                           aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                       <  2:223884/100‑1 (MQ=255)
                                                            cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                      >  2:529809/1‑100 (MQ=255)
                                                            cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                      >  2:1653603/1‑100 (MQ=255)
                                                             gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                     <  1:700412/100‑1 (MQ=255)
                                                              ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                    >  2:1008388/1‑100 (MQ=255)
                                                               cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                   >  1:523525/1‑100 (MQ=255)
                                                                  aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                <  1:833447/100‑1 (MQ=255)
                                                                  aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                <  2:955534/100‑1 (MQ=255)
                                                                    ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                              <  1:260684/100‑1 (MQ=255)
                                                                      ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                            <  2:1361718/100‑1 (MQ=255)
                                                                        ccTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                          <  1:1268324/100‑1 (MQ=255)
                                                                         cTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                         >  1:1590143/1‑100 (MQ=255)
                                                                          tAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                        <  1:244991/100‑1 (MQ=255)
                                                                            c‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                      >  2:238174/1‑100 (MQ=255)
                                                                            c‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                      >  1:1631186/1‑100 (MQ=255)
                                                                              aaaaTCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                     <  2:640780/100‑1 (MQ=255)
                                                                                aaTCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTTTTTTTTTTCAGTCTTGTATGGCTTTTAat                   >  1:108534/1‑100 (MQ=25)
                                                                                  tCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                 <  1:712486/100‑1 (MQ=255)
                                                                                   cGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                <  2:523526/100‑1 (MQ=255)
                                                                                   cGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                 <  1:672559/99‑1 (MQ=255)
                                                                                       cAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa            >  1:606341/1‑100 (MQ=255)
                                                                                        aaaTTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt           <  2:1238279/100‑1 (MQ=255)
                                                                                            tCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag       <  2:1547748/100‑1 (MQ=255)
                                                                                             cAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga      >  2:1429814/1‑100 (MQ=255)
                                                                                             cAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga      >  1:451239/1‑100 (MQ=255)
                                                                                             cAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga      >  1:1597095/1‑100 (MQ=255)
                                                                                              aaTACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc     >  1:212292/1‑100 (MQ=255)
                                                                                                 aCATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtt  <  2:116069/100‑1 (MQ=255)
                                                                                                 aCATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtt  >  2:979394/1‑100 (MQ=255)
                                                                                                  |                                                                                                     
GCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTT  >  NC_000913/3705920‑3706115

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