Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,520,352	C→T	P31S (CCT→TCT)	mcbR →	DNA‑binding transcriptional dual regulator McbR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,520,352	0	C	T	100.0%	269.7 / NA	75	P31S (CCT→TCT)	mcbR	DNA‑binding transcriptional dual regulator McbR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/43);  total (32/43)

TCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCG  >  NC_000913/1520215‑1520486                                                                                                                                          |                                                                                                                                       tctaccATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCt                                                                                                                                      >  2:533734/5‑140 (MQ=255) tCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCt                                                                                                                                      >  1:87515/1‑140 (MQ=255)              caggGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGa                                                                                                                         <  1:443252/140‑1 (MQ=255)               aggGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGAt                                                                                                                        <  1:842924/140‑1 (MQ=255)                   ttGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACt                                                                                                                    <  2:87515/140‑1 (MQ=255)                       aCAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTaaaa                                                                                                                >  2:760485/1‑140 (MQ=255)                           aaaaCACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCt                                                                                                            >  2:546132/1‑140 (MQ=255)                           aaaaCACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCt                                                                                                            >  2:253525/1‑140 (MQ=255)                           aaaaCACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCt                                                                                                            >  2:252312/1‑140 (MQ=255)                           aaaaCACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCt                                                                                                            >  2:199669/1‑140 (MQ=255)                           aaaaCACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCt                                                                                                            >  1:197406/1‑140 (MQ=255)                           aaaaCACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCt                                                                                                            >  1:823677/1‑140 (MQ=255)                           aaaaCACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCt                                                                                                            >  1:539135/1‑140 (MQ=255)                              acacGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGc                                                                                                         <  2:827481/140‑1 (MQ=255)                               cacGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCg                                                                                                        <  2:724188/140‑1 (MQ=255)                                acGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCgg                                                                                                       >  1:66966/1‑140 (MQ=255)                                 cgcgAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGa                                                                                                      <  1:1111697/140‑1 (MQ=255)                                 cgcgAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGa                                                                                                      <  1:760485/140‑1 (MQ=255)                                   cgAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGc                                                                                                    <  1:151838/140‑1 (MQ=255)                                     aaGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCaa                                                                                                  >  1:93573/1‑140 (MQ=255)                                      aGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAAt                                                                                                 <  1:1057933/140‑1 (MQ=255)                                      aGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAAt                                                                                                 <  2:735593/140‑1 (MQ=255)                                      aGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAAt                                                                                                 <  1:252312/140‑1 (MQ=255)                                          aTTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAgg                                                                                             >  2:1224553/1‑140 (MQ=255)                                          aTTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAgg                                                                                             >  1:857168/1‑140 (MQ=255)                                          aTTAAAATGCCTGGAACGGG‑AAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGt                                                                                            <  2:1144199/140‑1 (MQ=255)                                           ttAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGt                                                                                            <  2:539135/140‑1 (MQ=255)                                           ttAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGt                                                                                            <  2:767816/140‑1 (MQ=255)                                             aaaaTGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTAt                                                                                          <  1:1024433/140‑1 (MQ=255)                                             aaaaTGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTAt                                                                                          <  1:958175/140‑1 (MQ=255)                                             aaaaTGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTAt                                                                                          <  2:66966/140‑1 (MQ=255)                                             aaaaTGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTAt                                                                                          <  2:831816/140‑1 (MQ=255)                                              aaaTGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATg                                                                                         <  1:1224553/140‑1 (MQ=255)                                              aaaTGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATg                                                                                         <  1:199669/140‑1 (MQ=255)                                              aaaTGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATg                                                                                         <  1:644138/140‑1 (MQ=255)                                              aaaTGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATg                                                                                         <  1:546132/140‑1 (MQ=255)                                                aTGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAg                                                                                       <  1:1018927/140‑1 (MQ=255)                                                  gCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTa                                                                                     >  2:1181543/1‑140 (MQ=255)                                                     tGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAAt                                                                                                 >  1:858545/1‑125 (MQ=255)                                                     tGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAAt                                                                                                 <  2:858545/125‑1 (MQ=255)                                                     tGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTa                                                                                  >  1:137294/1‑140 (MQ=255)                                                             gAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATg                                                                                         <  1:161446/125‑1 (MQ=255)                                                             gAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATg                                                                                         >  2:161446/1‑125 (MQ=255)                                                                 aaTGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTg                                                                      <  2:576767/140‑1 (MQ=255)                                                                    gAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAg                                                                   <  2:1368365/140‑1 (MQ=255)                                                                       aCATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCat                                                                >  2:70812/1‑140 (MQ=255)                                                                           gTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCattatt                                                            >  2:814223/1‑140 (MQ=255)                                                                                tttGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTc                                                       <  1:814223/140‑1 (MQ=255)                                                                                               tGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATg                                        >  1:811554/1‑140 (MQ=255)                                                                                                 agaACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATgcg                                      <  2:823677/140‑1 (MQ=255)                                                                                                  gaACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATgcgc                                     >  1:1350644/1‑140 (MQ=255)                                                                                                   aaCGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCt                                    <  1:1241014/140‑1 (MQ=255)                                                                                                    aCGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCtt                                   <  2:464048/140‑1 (MQ=255)                                                                                                       aCCTGAAACATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCt                                >  2:322712/1‑140 (MQ=255)                                                                                                            atgtttCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGc                           >  1:68474/6‑140 (MQ=255)                                                                                                            aaaCATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGc                           <  1:253525/140‑1 (MQ=255)                                                                                                             aaCATCAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCa                          <  2:857168/140‑1 (MQ=255)                                                                                                                  cAGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGc                     >  2:691158/1‑140 (MQ=255)                                                                                                                   aGCTTAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGca                    <  1:1181543/140‑1 (MQ=255)                                                                                                                      ttAGTATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACaa                 <  1:295102/140‑1 (MQ=255)                                                                                                                          tATTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCAt             >  1:453797/1‑140 (MQ=255)                                                                                                                           aTTGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCAtt            <  2:197406/140‑1 (MQ=255)                                                                                                                             tGGCGCATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTa          <  2:706914/140‑1 (MQ=255)                                                                                                                               gcgcATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACa        <  2:1350644/140‑1 (MQ=255)                                                                                                                               gcgcATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACa        >  2:566402/1‑140 (MQ=255)                                                                                                                                  cATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACa        >  2:427533/1‑137 (MQ=255)                                                                                                                                  cATTAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACa        <  1:427533/137‑1 (MQ=255)                                                                                                                                    ttAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTa          >  2:636948/1‑133 (MQ=255)                                                                                                                                    ttAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTa          <  1:636948/133‑1 (MQ=255)                                                                                                                                    ttAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTcc   >  1:859494/1‑140 (MQ=255)                                                                                                                                    ttAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTcc   >  1:259577/1‑140 (MQ=255)                                                                                                                                    ttAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTcc   >  2:398895/1‑140 (MQ=255)                                                                                                                                     tAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCg  <  2:1019608/140‑1 (MQ=255)                                                                                                                                     tAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCg  <  1:1157101/140‑1 (MQ=255)                                                                                                                                     tAAAATCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCg  <  1:1339655/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       TCTGCCATGCAGGCAGGGTTGGACAGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGACTCTGCAGGTTGAGAACGACCTGAAACATCAGCTTAGTATTGGCGCATTAAAACCTGGCGCACGCCTGATTACTAAAAATCTGGCGGAGCAATTAGGTATGAGTATTACACCTGTGCGTGAAGCATTATTACGTCTGGTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCG  >  NC_000913/1520215‑1520486

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