Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 16.0% intergenic (‑181/+127) 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,9710CA16.0% 670.7 / 102.4 217intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (94/88);  new base A (32/3);  total (126/91)
Fisher's exact test for biased strand distribution p-value = 4.17e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

Only 100 of 226 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                     
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGTGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                    >  2:2436922/1‑100 (MQ=255)
   aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                   <  2:3017077/100‑1 (MQ=255)
    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGTACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                  <  1:327274/100‑1 (MQ=255)
     tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                 <  1:246775/100‑1 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                <  1:872903/100‑1 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGGCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                               >  1:1897366/1‑100 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                             >  1:920051/1‑100 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                           >  2:1474505/1‑100 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                           <  2:3364374/100‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                           <  2:3271094/100‑1 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                        >  1:943263/1‑100 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                        <  2:633553/100‑1 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                        >  1:1215480/1‑100 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                       >  2:38623/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                     >  1:151164/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                     >  2:1244809/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                     <  2:842695/100‑1 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgn                                                                                   <  1:2039563/100‑2 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                   >  2:1687999/1‑100 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                 <  2:3153067/100‑1 (MQ=255)
                      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                <  2:2714226/100‑3 (MQ=37)
                       tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                               >  2:2918913/1‑100 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                            <  2:2203751/100‑1 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                              >  1:2172834/1‑99 (MQ=255)
                          cgcTGCACGCAGAGGCACATCGCCACCGGGTATACAGATTGTTAACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                            >  2:734554/1‑99 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACGAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                           >  1:1928874/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                           >  1:967196/1‑100 (MQ=255)
                           gttGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                          >  2:2338309/3‑100 (MQ=255)
                            cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                         >  2:2641531/1‑100 (MQ=255)
                            cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                         >  2:378204/1‑100 (MQ=255)
                             tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                        >  1:338329/1‑100 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                      <  2:2563166/100‑1 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                      <  2:3343962/100‑1 (MQ=255)
                                aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                     <  1:1138382/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                    >  1:457585/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                    <  2:1902549/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                    >  2:252702/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                    <  1:847513/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                   <  2:1358289/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                   >  1:2269018/1‑100 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                            <  2:1145506/100‑1 (MQ=255)
                                                 cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                    >  1:2231321/1‑100 (MQ=255)
                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                   <  1:1988304/100‑1 (MQ=255)
                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                   >  1:2659629/1‑100 (MQ=255)
                                                       gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAt                                              <  2:2238163/100‑1 (MQ=255)
                                                          taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAt‑ttt                                          >  2:2394041/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCg                                        >  2:201333/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCg                                        >  2:1981182/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTAGTTTATGCCGGATACGGCGTTAACGCCTTATCCGTCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAt‑ttt                                          >  1:79632/1‑98 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGt                                       <  2:2502859/100‑1 (MQ=255)
                                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTa                                      <  2:1603235/100‑1 (MQ=255)
                                                               aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAg                                     <  1:1630935/100‑1 (MQ=255)
                                                               aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAg                                     >  1:2877278/1‑100 (MQ=255)
                                                                ttgttgCCAAAGTGAGCTTTCTTTATGCCTGATGAGGCGTTAGCGGCTTATCCGGCGTAC‑GAAGTCGTGCAAATTCATTACATTGCAGAAT‑TGTCGtcag                                    >  1:548177/1‑97 (MQ=255)
                                                                ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAgg                                    <  1:780896/100‑1 (MQ=255)
                                                                 cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGc                                   >  1:10325/6‑100 (MQ=18)
                                                                 cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGc                                   >  2:1061401/6‑100 (MQ=18)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGcc                                  <  2:128359/100‑1 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCt                                 >  1:732815/1‑100 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGa                               >  1:3166985/1‑100 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGAt                              <  2:754842/100‑1 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTGATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGAt                              >  1:2416055/1‑100 (MQ=255)
                                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCg                         <  2:435/100‑1 (MQ=255)
                                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCg                         <  2:2467419/100‑1 (MQ=255)
                                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCg                         >  1:3461884/1‑100 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCgg                        >  1:96130/1‑100 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGa                       <  1:2338290/100‑1 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAg                      <  2:3317055/100‑1 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAg                      >  2:2892393/1‑100 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGGCTTATCCGGCCGAC‑AAAATCGTGCAAATTCAATACATTACAGCAT‑TTTCGTCGGCATGATAAGCGGAgc                     >  2:1468100/1‑100 (MQ=18)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAgc                     >  1:410848/1‑100 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAgc                     <  2:943270/100‑1 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAgcg                    <  1:63398/100‑1 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAgcg                    >  2:876262/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAgcg                    <  2:2875981/100‑1 (MQ=255)
                                                                                 tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATCGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAgcgc                   <  1:2283542/100‑1 (MQ=255)
                                                                                   tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGAGT‑TCACGTAGGCCTGATAAGCGTAGCGCAt                 >  2:159091/8‑100 (MQ=17)
                                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCAGGAT‑ACGTGTAGGCCTGATAAGCGTAGCGCAt                 >  1:295602/8‑100 (MQ=18)
                                                                                    ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTT‑ATTGGTAGGCCTGATAAGCGCAGCGCATc                >  2:2006981/7‑100 (MQ=21)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                >  1:2009318/1‑100 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa               <  2:2599090/100‑1 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa               <  1:2812776/100‑1 (MQ=255)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg              <  1:1474492/100‑1 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:2085240/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGAGT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:1300739/4‑100 (MQ=17)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg             >  2:532454/1‑100 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:3111786/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1027473/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGAGT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1201162/3‑100 (MQ=17)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCAGGAA‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGc            >  1:2598759/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc            >  2:1264129/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:1266443/3‑100 (MQ=18)
                                                                                        aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc            >  2:3124770/1‑100 (MQ=255)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGGAG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:771773/2‑100 (MQ=17)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGGAG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:2733161/1‑100 (MQ=21)
                                                                                           ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         >  2:2185395/1‑100 (MQ=255)
                                                                                           ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         <  2:2777647/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:1492103/100‑1 (MQ=255)
                                                                                                tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  1:2964701/100‑1 (MQ=255)
                                                                                                tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  2:2255934/100‑1 (MQ=255)
                                                                                                  |                                                                                                     
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070

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