Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 28.6% intergenic (‑195/+113) 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,9850TC28.6% 448.7 / 207.1 257intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (92/91);  new base C (69/4);  total (161/96)
Fisher's exact test for biased strand distribution p-value = 9.59e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-03

Only 100 of 261 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084
                                                                                                   |                                                                                                             
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                               >  1:70557/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                              >  2:2924243/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGACGTTAACGCCTTa                                                                                                              >  1:605633/1‑100 (MQ=255)
   tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                            <  1:2080456/100‑1 (MQ=255)
    gcCGGACGGCGCTGCCCGCAGAGGCACAGCGCCACCGGGTATACAGACTGTTGACACAGTGCGCTTCGTTCCTGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                           <  1:2712024/100‑1 (MQ=37)
     ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                          >  2:3311636/1‑100 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                           >  1:3058352/1‑98 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                         >  1:2863805/1‑100 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                        <  2:677890/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                       <  2:2483789/100‑1 (MQ=255)
          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                     >  1:1392216/1‑100 (MQ=255)
           ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                                    <  2:1912245/100‑2 (MQ=25)
              cctgtaggaagagGCGCAGTGCCACCGGGTATACAGATTGCTGACAAAGTGAGCTCTGATTATGCCCGATGCGTCGTTGACGCCTTATCCGGCCTAc                                                                                                    >  1:1775629/10‑97 (MQ=17)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                              <  2:2583514/100‑1 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                            <  1:3311611/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                          <  2:2370937/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                          <  2:2634915/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                          >  2:2055967/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                          <  1:2677849/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                         <  1:3608290/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                        <  2:1528133/100‑1 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                   <  1:2930836/100‑1 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                   <  1:903891/100‑1 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                                  >  2:1215504/1‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                   >  1:3240987/1‑99 (MQ=255)
                              caGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACTCCTTATCCGGCCTAC‑AAAATCCTGCAAATTCAAt                                                                                <  1:808522/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                                <  2:1107637/100‑1 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                                <  1:767953/99‑1 (MQ=255)
                                gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                              >  2:504571/1‑100 (MQ=255)
                                  gcCACCGAGTATGCAGATTGTTGACAAAGTGCTCTTTGTGTATGCCCGATGCGGCGTTAACGCCTTAGCAGGCCTAC‑AAGCTCGTGCAAATTCAAt                                                                                <  1:1725959/96‑1 (MQ=17)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                           >  1:2629182/1‑100 (MQ=255)
                                     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                         <  2:1392222/100‑1 (MQ=255)
                                     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                         <  1:316794/100‑1 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                                      <  1:3509354/100‑1 (MQ=255)
                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑At‑t                                                               >  1:1992360/1‑100 (MQ=255)
                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑At‑t                                                               >  1:3164643/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑g                                                      <  1:2975994/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑g                                                      >  2:3229543/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑g                                                      <  2:410958/100‑1 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑Gt                                                     <  2:3275461/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGcc                                                >  2:1344721/1‑100 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGcc                                                <  1:1002677/100‑1 (MQ=255)
                                                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCt                                               <  1:2213313/100‑1 (MQ=255)
                                                       tgACAAAGTGAACTTTGTTTATGCCGGATGCGGCATTAACGCCTTATCCGGCCTAC‑AAAATCGTGAAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAgcc                                                 >  1:1433280/1‑95 (MQ=17)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                           >  1:1070408/1‑100 (MQ=255)
                                                          cAAAGTGCGCTTGGTTTATGCCGGATGCGGCGTTCACGTCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                           <  2:74374/100‑1 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                      <  1:2655214/100‑1 (MQ=255)
                                                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                     <  2:2290456/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                   >  2:2311046/1‑100 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  2:2237953/100‑1 (MQ=255)
                                                                     atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                               >  2:1236679/9‑100 (MQ=18)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                            >  2:3065013/1‑100 (MQ=255)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑‑AT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  2:115973/5‑100 (MQ=21)
                                                                          taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑‑A‑‑AT‑T‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:3142398/4‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:1631766/4‑100 (MQ=18)
                                                                          ctggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAg                            <  1:1322762/97‑1 (MQ=21)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:2006942/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1627716/3‑100 (MQ=18)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                         >  2:3435657/1‑100 (MQ=255)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCA‑ATTCAA‑C‑‑‑TT‑‑‑GTAGGCCTGATAAGCGCAGCGCATCAgag                          >  2:1322768/2‑98 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                         >  2:1118406/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑‑CT‑C‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:1641664/2‑99 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑‑AA‑A‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                         >  1:3572448/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑‑GA‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                         >  1:1164148/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAG‑AT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2084962/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑‑GA‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2442943/2‑100 (MQ=17)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑GT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:2602012/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:365000/1‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑‑A‑‑AT‑T‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                        >  2:350184/1‑98 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑‑A‑‑AT‑T‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                        >  1:2134425/1‑98 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCA‑TG‑‑GG‑C‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:913008/1‑99 (MQ=17)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:2106601/1‑99 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:2698085/1‑100 (MQ=21)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑‑GC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  2:2789724/9‑98 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:2334547/9‑98 (MQ=21)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                      >  2:1014758/9‑100 (MQ=21)
                                                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  2:3079328/1‑100 (MQ=255)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1898626/7‑100 (MQ=17)
                                                                                ggataaggcgTTTACGCAGCATCCGGCCTAC‑AAAAGCATGCAAATTCGATATATTGCA‑GA‑‑GT‑T‑‑‑CAC‑‑GTAAGCCTGATAAGCGTAGCGCATAAGGCAgtt                      >  1:547477/7‑99 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                     >  2:1478204/7‑96 (MQ=25)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                     >  1:3622922/7‑96 (MQ=25)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑‑GA‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  2:64034/7‑100 (MQ=17)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑‑CT‑C‑‑‑GAC‑‑GTAGGCCTGATAAGCATAGAGCATCAAGCActttt                    >  1:1583303/6‑100 (MQ=18)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑‑AA‑A‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgtttt                    >  1:1604488/6‑100 (MQ=18)
                                                                                    gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                 <  1:970168/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAATTTAAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAATCATGATAAGCGGAGTGAATCAGGCATTcttgcct                >  1:2553076/1‑100 (MQ=18)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  1:3470182/1‑100 (MQ=255)
                                                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               >  2:2744096/1‑100 (MQ=255)
                                                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               >  1:3435918/1‑100 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  1:3014201/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:851893/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:359905/100‑1 (MQ=255)
                                                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  1:2484278/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  1:2959480/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  1:1243904/1‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       >  1:3139547/1‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       >  1:1510079/1‑100 (MQ=255)
                                                                                              cGCCGTATCCGGCCTAC‑AAAATCGTGCAAATTCAAGACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:1051884/99‑1 (MQ=255)
                                                                                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  2:1741140/1‑100 (MQ=255)
                                                                                                   tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  2:1867324/1‑100 (MQ=255)
                                                                                                   |                                                                                                             
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084

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: ‑