Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	741,597	T→G	A320A (GCA→GCC)	dtpD ←	dipeptide and tripeptide permease D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	741,597	0	T	G	99.0%	340.4 / NA	97	A320A (GCA→GCC)	dtpD	dipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (45/51);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGCTA  >  NC_000913/741500‑741693                                                                                                  |                                                                                                 caaggccaagAGCAAATTTCCCCCGGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGc                                                                                                >  1:2823747/1‑100 (MQ=255)     gccaagAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAg                                                                                            >  2:527584/1‑100 (MQ=255)         agagCAAATTTCCCCCCGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCCCCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCAt                                                                                        <  2:1480242/100‑1 (MQ=255)          gagCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCAtt                                                                                       >  1:483393/1‑100 (MQ=255)            gCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTaa                                                                                     <  1:625873/100‑1 (MQ=255)            gCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTaa                                                                                     >  2:2307015/1‑100 (MQ=255)             cAAATTTCCCCCCGATGCGCACGGTACAATTACCCGCGACAATCTCTTTTACCACCCACGCACGGAACACACCGCACAGCATAACGGCGAAGGCATTAAt                                                                                    <  2:1085996/100‑1 (MQ=255)             cAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAAt                                                                                    >  1:3968036/1‑100 (MQ=255)              aaaTTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATc                                                                                   <  2:2142754/100‑1 (MQ=255)               aaTTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCg                                                                                  >  2:235878/1‑100 (MQ=255)                aTTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGa                                                                                 >  2:245742/1‑100 (MQ=255)                 tttCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGAc                                                                                <  2:3674366/100‑1 (MQ=255)                  ttCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACt                                                                               <  2:3124843/100‑1 (MQ=255)                     ccccAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGa                                                                            <  2:1433892/100‑1 (MQ=255)                      tccAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGaa                                                                           <  2:1294248/99‑1 (MQ=255)                      cccAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGan                                                                           <  1:3261284/100‑2 (MQ=255)                       ccAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGaaa                                                                          >  1:1318260/1‑99 (MQ=255)                       ccAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAAc                                                                          <  2:826114/100‑1 (MQ=255)                       ccAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAAc                                                                          <  2:956916/100‑1 (MQ=255)                         aGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACAt                                                                        <  1:2799527/100‑1 (MQ=255)                            tGCGCACGGTACGATTACCCGCGACCCTCTCTTTTACCACCCACGCCAGGACCACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGATCATcgc                                                                     <  2:319954/100‑1 (MQ=255)                            tGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATcgc                                                                     <  1:2958577/100‑1 (MQ=255)                             gcgcACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATcgcg                                                                    >  2:2482993/1‑100 (MQ=255)                              cgcACGGTACGATTACCCGCGACACTCTCTTTTACCTCCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCgg                                                                   <  2:1364396/100‑1 (MQ=255)                              cgcACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCgg                                                                   <  2:827023/100‑1 (MQ=255)                              cgcACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCgg                                                                   <  1:197369/100‑1 (MQ=255)                               gcACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGt                                                                  <  2:2537335/100‑1 (MQ=255)                                cACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTc                                                                 >  2:2525084/1‑100 (MQ=255)                                cACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTc                                                                 >  2:2525119/1‑100 (MQ=255)                                cACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTc                                                                 >  1:1582113/1‑100 (MQ=255)                                 aCGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCg                                                                <  2:3048160/100‑1 (MQ=255)                                  cGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCgg                                                               <  1:1984514/100‑1 (MQ=255)                                  cGGTACGATTACCCGCGACACTCTCTTATACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCgg                                                               <  1:3280715/100‑1 (MQ=255)                                    gTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGaa                                                             >  1:63599/1‑100 (MQ=255)                                       cGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACgg                                                          <  1:3916949/100‑1 (MQ=255)                                       cGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACgg                                                          >  1:3065135/1‑100 (MQ=255)                                         aTTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGta                                                        <  1:1258238/100‑1 (MQ=255)                                          ttACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATACCGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGtat                                                       <  1:1144811/100‑1 (MQ=255)                                           tACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGtata                                                      <  1:1283089/100‑1 (MQ=255)                                            aCCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATaa                                                     >  1:4039595/1‑100 (MQ=255)                                              ccGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAAcc                                                   >  2:1558470/1‑100 (MQ=255)                                               cgcgACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAAGGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGCAACGGTATAACCa                                                  >  1:3680658/1‑100 (MQ=255)                                                  gACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAAc                                               <  1:2405812/100‑1 (MQ=255)                                                  gACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAAc                                               <  1:405088/100‑1 (MQ=255)                                                  gACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAAc                                               <  1:1549883/100‑1 (MQ=255)                                                  gACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAAc                                               <  2:2171376/100‑1 (MQ=255)                                                     acTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACata                                            >  1:3137941/1‑100 (MQ=255)                                                      ctctctTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACatat                                           >  1:3464178/1‑100 (MQ=255)                                                      ctctctTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACatat                                           >  1:2389863/1‑100 (MQ=255)                                                      ctctctTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACatat                                           >  2:3593904/1‑100 (MQ=255)                                                      ctctctTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACatat                                           >  1:1733311/1‑100 (MQ=255)                                                      ctctctTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACatat                                           >  2:1007559/1‑100 (MQ=255)                                                       tctctTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATc                                          <  1:1752022/100‑1 (MQ=255)                                                          ctTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGCAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATcgcg                                       >  2:1430677/1‑100 (MQ=255)                                                          ctTTTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATcgcg                                       >  1:77937/1‑100 (MQ=255)                                                          ctTCTACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATcgcg                                       <  2:131519/100‑1 (MQ=255)                                                           ttttACCACCTACGCCAGTAACACACCGCACAGCATAACGGAGAAGGCATTAATCGTCTGGAACATCACTGTTGGAACGGTATAACCAAACATATCGCgg                                      >  2:602654/1‑100 (MQ=255)                                                           ttttACCACCCACGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCgg                                      <  1:710375/100‑1 (MQ=255)                                                                ccacccacGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAAc                                 >  1:3487184/1‑100 (MQ=255)                                                                ccacccacGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAAc                                 <  2:3352439/100‑1 (MQ=255)                                                                    ncacGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAg                             >  1:3251439/2‑100 (MQ=255)                                                                    ccacGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAg                             >  1:3198223/1‑100 (MQ=255)                                                                     cacGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTTGGAACGGTATAACCAAACATATCGCGGTTAACGAAGc                            <  2:3470725/100‑1 (MQ=255)                                                                     cacGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGc                            <  1:1433831/100‑1 (MQ=255)                                                                     cacGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGc                            <  2:205197/100‑1 (MQ=255)                                                                       cGCCAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCgg                          <  2:3482085/100‑1 (MQ=255)                                                                         ccAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTc                        <  2:3703543/100‑1 (MQ=255)                                                                         ccAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTc                        <  1:3603397/100‑1 (MQ=255)                                                                         ccAGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTc                        >  1:3734868/1‑100 (MQ=255)                                                                          cAGGAACACACCGCACAGCATAACGGCGAAGGCATTACTCGACTGGAACATCGCGGTCGGAACGGTATACCCAAACATATCACGGTTACCGAAGAGGTCg                       >  2:2398509/1‑100 (MQ=255)                                                                          cAGGAACACACCGCAATGCATAACTGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACCTATCGCGGTTAACGAAGCGGTCg                       <  2:3352718/100‑1 (MQ=255)                                                                           aGGAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGa                      <  2:1028020/100‑1 (MQ=255)                                                                            ggAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGat                     <  2:2191546/100‑1 (MQ=255)                                                                            ggAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGat                     <  1:3063897/100‑1 (MQ=255)                                                                            ggAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGat                     <  1:403113/100‑1 (MQ=255)                                                                             gAACACACCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGata                    <  2:3427027/100‑1 (MQ=255)                                                                               acacacCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGAtata                  >  2:3889924/1‑100 (MQ=255)                                                                               acacacCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGAtata                  <  1:4043674/100‑1 (MQ=255)                                                                                cacacCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAATCATATCGCGGTTAACGATGCGGTCGATATaa                 >  2:3131503/1‑100 (MQ=255)                                                                                cacacCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATaa                 <  1:858210/100‑1 (MQ=255)                                                                                 acacCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATaaa                >  2:2671733/1‑100 (MQ=255)                                                                                 acacCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATaaa                >  1:2802073/1‑100 (MQ=255)                                                                                 acacCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATaaa                <  2:2045173/100‑1 (MQ=255)                                                                                  cacCGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAg               >  1:2112545/1‑100 (MQ=255)                                                                                    ccGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACACATCGCGGTTAACGAAGCGGTCGATATAAAGGc             <  2:3602197/100‑1 (MQ=255)                                                                                     cGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCt            <  2:2432527/100‑1 (MQ=255)                                                                                     cGCACAGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCt            >  2:2584368/1‑100 (MQ=255)                                                                                         caGCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATc        >  2:2289766/1‑100 (MQ=255)                                                                                          aGCATAACGGCGAAGGCATTAATTGACTGGAACATCGCGGTCAGAACGGTATAACCAAACCTATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCg       <  1:3557260/100‑1 (MQ=255)                                                                                           gCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGa      >  2:124555/1‑100 (MQ=255)                                                                                           gCATAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGa      >  1:834691/1‑100 (MQ=255)                                                                                             aTAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGc    >  2:2785993/1‑100 (MQ=255)                                                                                             aTAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGc    <  1:2986098/100‑1 (MQ=255)                                                                                             aTAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGc    >  1:1378982/1‑100 (MQ=255)                                                                                             aTAACGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGc    >  1:4034888/1‑100 (MQ=255)                                                                                              tAACGGCGAAGGCAATAATCGACTGGGACATCGCGGTCGGAACGGTAAAACCAAGCATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGCt   >  2:3739357/1‑100 (MQ=255)                                                                                               aaCGGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGCTa  >  1:1494920/1‑100 (MQ=255)                                                                                                  |                                                                                                 CAAGGCCAAGAGCAAATTTCCCCCAGATGCGCACGGTACGATTACCCGCGACACTCTCTTTTACCACCCACGCCAGGAACACACCGCACAGCATAACTGCGAAGGCATTAATCGACTGGAACATCGCGGTCGGAACGGTATAACCAAACATATCGCGGTTAACGAAGCGGTCGATATAAAGGCTAATCGAGCTA  >  NC_000913/741500‑741693

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