Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,632,949	C→T	A20A (GCC→GCT)	ynfT →	Qin prophage; protein YnfT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,632,949	0	C	T	100.0%	201.6 / NA	84	A20A (GCC→GCT)	ynfT	Qin prophage; protein YnfT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/40);  total (44/40)

CGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCAATA  >  NC_000913/1632818‑1633051                                                                                                                                    |                                                                                                       cGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATaaa                                                                                                >  1:996618/1‑140 (MQ=25)  gAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATaaaa                                                                                               >  2:559020/1‑140 (MQ=25)   aaCTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAAt                                                                                              >  2:378786/1‑140 (MQ=25)   aaCTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAAt                                                                                              >  2:246923/1‑140 (MQ=25)    aCTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAAtt                                                                                             <  1:237368/140‑1 (MQ=25)        tAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATTAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAggg                                                                                         <  1:1106055/140‑1 (MQ=255)        tAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAggg                                                                                         >  2:185329/1‑140 (MQ=25)         aaGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAGAATTAGGGa                                                                                        <  2:563863/140‑1 (MQ=255)         aaGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGa                                                                                        <  2:247454/140‑1 (MQ=21)         aaGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGa                                                                                        <  2:1159673/140‑1 (MQ=25)              cacCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGc                                                                                   >  2:419201/1‑140 (MQ=25)               acCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCa                                                                                  <  2:738162/140‑1 (MQ=21)               acCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCa                                                                                  <  2:200519/140‑1 (MQ=25)                cccAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCat                                                                                 <  2:384653/140‑1 (MQ=37)                cccAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCat                                                                                 >  2:126724/1‑140 (MQ=25)                          aaaTTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAAc                                                                       <  1:1007240/140‑1 (MQ=17)                          aaaTTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAAc                                                                       <  1:696704/140‑1 (MQ=14)                               gTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGa                                                                  >  1:428177/1‑140 (MQ=18)                                     tataTCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAAc                                                            <  1:378786/140‑1 (MQ=14)                                     tataTCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCATTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAAc                                                            <  1:401535/140‑1 (MQ=9)                                          cATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATg                                                       <  1:47156/140‑1 (MQ=17)                                            tatataACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGta                                                     <  1:419201/140‑1 (MQ=14)                                               ataACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATaa                                                  >  1:673024/1‑140 (MQ=18)                                                taACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAg                                                 <  1:559020/140‑1 (MQ=14)                                                taACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAg                                                 <  1:593000/140‑1 (MQ=14)                                                taACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAg                                                 <  2:996618/140‑1 (MQ=14)                                                taACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAg                                                 <  1:128604/140‑1 (MQ=17)                                                 aaCACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCg                                                          >  2:1398393/1‑130 (MQ=21)                                                 aaCACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCg                                                          <  1:1398393/130‑1 (MQ=14)                                                 aaCACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGt                                                >  1:189788/1‑140 (MQ=17)                                                 aaCACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGt                                                >  1:1350108/1‑140 (MQ=18)                                                 aaCACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGt                                                >  2:1106206/1‑140 (MQ=18)                                                 aaCACATTACTAATTCGAGGCTATATGAACAACATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGt                                                >  1:186465/1‑140 (MQ=17)                                                   cacaTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCa                                              >  2:1306043/1‑140 (MQ=18)                                                     caTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAAc                                            >  2:648544/1‑140 (MQ=18)                                                      aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGta                                                     <  1:1238365/130‑1 (MQ=14)                                                      aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGta                                                     >  2:1238365/1‑130 (MQ=21)                                                      aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACa                                           >  1:508288/1‑140 (MQ=18)                                                      aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACa                                           >  2:148047/1‑140 (MQ=18)                                                      aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACa                                           >  1:786049/1‑140 (MQ=18)                                                      aTTACTAATTCGAGGCTATATGAACAACATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACa                                           >  2:590417/1‑140 (MQ=16)                                                        tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACata                                         >  1:1368088/1‑140 (MQ=18)                                                        tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACata                                         >  2:268381/1‑140 (MQ=18)                                                        tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACata                                         >  2:1321439/1‑140 (MQ=18)                                                          cTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATc                                       <  2:1352489/140‑1 (MQ=17)                                                          cTAATTCGAGGCTATATGAACAACATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATc                                       <  1:1480311/140‑1 (MQ=17)                                                             aTTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAg                                                 <  2:1116183/127‑1 (MQ=14)                                                             aTTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAg                                                 >  1:1116183/1‑127 (MQ=21)                                                             aTTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTg                                    >  1:1366159/1‑140 (MQ=18)                                                              ttCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGa                                   <  1:1464079/140‑1 (MQ=17)                                                                 gAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATc                                <  1:805389/140‑1 (MQ=17)                                                                 gAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATc                                <  2:673024/140‑1 (MQ=17)                                                                  aGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCa                               >  1:303788/1‑140 (MQ=18)                                                                  aGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCa                               >  2:49865/1‑140 (MQ=18)                                                                  aGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCa                               >  2:77371/1‑140 (MQ=18)                                                                   ggCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAg                              <  2:1368088/140‑1 (MQ=17)                                                                   ggCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAg                              <  1:648544/140‑1 (MQ=17)                                                                   ggCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAg                              <  1:182427/140‑1 (MQ=17)                                                                     cTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGAc                            <  2:896031/140‑1 (MQ=17)                                                                      tataTGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACa                           <  1:1106206/140‑1 (MQ=17)                                                                      tataTGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACa                           <  2:303788/140‑1 (MQ=17)                                                                        tatGAACAGCTTACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATa                         <  1:543329/140‑1 (MQ=17)                                                                             aCAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAg                              <  2:1138070/130‑1 (MQ=17)                                                                             aCAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAg                              >  1:1138070/1‑130 (MQ=18)                                                                                gCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCg                 >  1:1098030/1‑140 (MQ=18)                                                                                gCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCg                 >  2:1488605/1‑140 (MQ=18)                                                                                 cATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAAt                                 <  2:1359580/123‑1 (MQ=17)                                                                                 cATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAAt                                 >  1:1359580/1‑123 (MQ=14)                                                                                 cATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGc                >  2:559405/1‑140 (MQ=17)                                                                                 cATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGc                >  2:31849/1‑140 (MQ=18)                                                                                    aCTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAAt             >  2:291420/1‑140 (MQ=18)                                                                                    aCTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAAt             >  2:924031/1‑140 (MQ=18)                                                                                    aCTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAAt             >  2:1238817/1‑140 (MQ=18)                                                                                      tGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACa                           >  2:1106843/1‑124 (MQ=18)                                                                                      tGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACa                           <  1:1106843/124‑1 (MQ=17)                                                                                        aTAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGaaa         <  1:717835/140‑1 (MQ=17)                                                                                          aaTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAAt       <  2:508288/140‑1 (MQ=17)                                                                                           aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATc      >  1:84543/1‑140 (MQ=18)                                                                                           aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATc      >  2:320474/1‑140 (MQ=18)                                                                                             cacaTCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCaa    >  2:641716/1‑140 (MQ=17)                                                                                             cacaTCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCaa    >  1:373659/1‑140 (MQ=18)                                                                                               caTCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATGTCCAATAACCCAGACAAGACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCaata  >  2:1356558/1‑140 (MQ=2)                                                                                               caTCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCaata  <  1:294881/140‑1 (MQ=14)                                                                                               caTCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCaata  <  1:1306043/140‑1 (MQ=14)                                                                                                                                    |                                                                                                       CGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCAATA  >  NC_000913/1632818‑1633051

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