Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 41.4% 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.4% 376.9 / inf 322intergenic (‑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 (96/92);  new base G (122/11);  total (219/103)
Fisher's exact test for biased strand distribution p-value = 1.13e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-02

Only 100 of 326 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                            
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGCGCGCTGTGTTTATGACGGATGCGGCGTTAACGCCTTa                                                                                                            >  2:2825249/1‑100 (MQ=255)
   ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                           >  2:3816725/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                         <  2:2707694/100‑1 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                         >  2:490017/1‑100 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                       <  2:1055623/100‑1 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                       >  1:146899/1‑100 (MQ=255)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                                >  2:305772/1‑100 (MQ=255)
               gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                              <  2:3473602/100‑1 (MQ=255)
                    aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                          <  2:2531056/99‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                        <  1:2092396/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                       >  1:1012826/1‑100 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                       <  1:2334001/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                       <  1:244550/100‑1 (MQ=255)
                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                 <  1:1548548/100‑1 (MQ=255)
                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                 <  1:1277517/100‑1 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGTGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                                >  1:627757/1‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                              <  1:1799076/100‑1 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                            >  2:180282/1‑100 (MQ=255)
                                  cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCCTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACa                                                                           >  2:2209725/1‑100 (MQ=255)
                                           gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                       <  2:3732767/95‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑Gt                                                    >  1:1320482/1‑100 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTa                                                   >  1:3090983/1‑100 (MQ=255)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGcc                                               >  2:1347261/1‑100 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTg                                             >  2:4065163/1‑100 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGa                                            >  1:3346544/1‑100 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCAGAtcggaag                                     >  1:1856687/1‑94 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACACTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGTCTGATAAGCgggg                                   <  1:2825256/100‑3 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑CT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                  >  2:1645393/1‑99 (MQ=255)
                                                                        ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑AT‑CC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  1:3675378/7‑100 (MQ=18)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            <  2:3456550/100‑1 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            <  2:1601408/100‑1 (MQ=255)
                                                                         tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:1102970/6‑100 (MQ=18)
                                                                         atctggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCa                            <  1:654367/95‑1 (MQ=18)
                                                                          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           <  1:1591871/100‑1 (MQ=255)
                                                                          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           <  1:3412513/100‑1 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑‑A‑TTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:252764/5‑100 (MQ=18)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:1031472/5‑100 (MQ=18)
                                                                           taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          >  2:2482918/4‑100 (MQ=17)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          <  2:2627754/100‑1 (MQ=255)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           <  2:1484001/99‑1 (MQ=255)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:1979436/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:2978924/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:1359021/3‑100 (MQ=21)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑AT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:2197207/3‑100 (MQ=17)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                        <  2:3006328/100‑1 (MQ=255)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  1:3628725/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:3237597/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:3904255/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                         >  1:3585531/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAACATATTGCA‑AT‑TC‑T‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:776051/2‑98 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:472568/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑C‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1888515/2‑100 (MQ=18)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       >  1:2650065/1‑100 (MQ=255)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑AT‑TC‑T‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg                       >  2:3868412/1‑98 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑‑A‑TTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg                       >  2:1423190/1‑98 (MQ=18)
                                                                              gCCGGATAAGGCGTTCGCGCCGCATCCGGCCTAC‑AAAAGTTTGCAAATTCAATAAATTGCA‑GA‑AT‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAt                       >  2:4191360/1‑100 (MQ=17)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat                      >  1:1712766/1‑97 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:1309585/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑A‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCag                       >  2:3626007/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑AT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:1057376/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:1199639/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                         >  2:3690750/1‑98 (MQ=17)
                                                                               ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                       >  1:92542/9‑99 (MQ=21)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCCTTGCGAATTCAATATATTGTATGG‑AA‑A‑‑‑TTC‑‑‑‑AGGCCTGATAAGCGTAGCGCATCAGGCtttt                     >  2:1727708/9‑100 (MQ=25)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:914270/9‑98 (MQ=21)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                      >  1:4028903/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                     >  2:2163276/9‑100 (MQ=21)
                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1190945/1‑100 (MQ=255)
                                                                                cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                     >  2:2109531/8‑97 (MQ=17)
                                                                                cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                     >  1:396157/1‑100 (MQ=255)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                   >  1:2464264/7‑98 (MQ=21)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:3278712/7‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:3983671/7‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                     >  1:3842269/7‑97 (MQ=21)
                                                                                 ggataaggcgTTCACGCCTTATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GA‑GA‑T‑‑‑GGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                    >  1:3540081/7‑100 (MQ=14)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                     >  2:3204871/7‑96 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:861039/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2400924/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:3599403/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2149184/7‑100 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                    >  1:2568786/7‑96 (MQ=25)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGACTGATAAGCGTAGCGCATCAGGCAattc                    >  1:771647/7‑99 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                    >  1:2172604/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCA‑CA‑GC‑T‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                    >  1:1073558/7‑96 (MQ=37)
                                                                                  gctaaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                  >  1:2848181/6‑97 (MQ=21)
                                                                                    taaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattctgt                 >  1:4071370/4‑93 (MQ=25)
                                                                                    tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTTGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 >  1:3771128/1‑100 (MQ=255)
                                                                                        gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttttttt             >  2:2873620/1‑94 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:3091136/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:2731698/100‑1 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           >  2:810154/1‑100 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  1:1497072/100‑1 (MQ=255)
                                                                                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCAGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  2:1374497/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  1:448589/100‑1 (MQ=255)
                                                                                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  1:1913083/1‑99 (MQ=255)
                                                                                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      >  1:493972/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     >  2:2655945/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:3713940/100‑1 (MQ=255)
                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:3298878/100‑1 (MQ=255)
                                                                                                   ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  1:738607/100‑1 (MQ=255)
                                                                                                   |                                                                                                            
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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