Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,658,893	G→A	*176* (TAG→TAA)	gadE →	gad regulon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,658,893	0	G	A	100.0%	300.2 / NA	92	*176* (TAG→TAA)	gadE	gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (46/46);  total (46/46)

GTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAGTAGAAC  >  NC_000913/3658776‑3659027                                                                                                                       |                                                                                                                                       gTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACCGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGaa                                                                                                                    >  2:797437/1‑139 (MQ=255)   cGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAAtt                                                                                                                 >  2:1329273/1‑139 (MQ=255)     aCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGa                                                                                                               >  2:359445/1‑139 (MQ=255)      cGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGc                                                                                                                           <  1:694463/126‑1 (MQ=255)      cGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGc                                                                                                                           >  2:694463/1‑126 (MQ=255)        cTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                            <  2:1416071/139‑1 (MQ=255)           aaTATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggctg                                                                                                         >  1:483120/1‑139 (MQ=255)                  cGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTTCAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGaaa                                                                                                                   >  1:748254/1‑122 (MQ=255)                  cGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGaaa                                                                                                                   <  2:748254/122‑1 (MQ=255)                   gCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                                 <  1:1329273/139‑1 (MQ=255)                   gCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                                 <  2:759387/139‑1 (MQ=255)                   gCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                                 >  2:1032079/1‑139 (MQ=255)                    caacaaACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACCGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                                 <  1:797437/139‑1 (MQ=255)                          acaTTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATg                                                                                          <  1:437569/139‑1 (MQ=255)                            aTTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCa                                                                                        >  1:676187/1‑139 (MQ=255)                                aaGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATg                                                                                    <  2:572545/139‑1 (MQ=255)                                     tCAGCAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATg                                                                                                    <  2:1159846/118‑1 (MQ=255)                                     tCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGACTGGCATg                                                                                                    >  1:1159846/1‑118 (MQ=255)                                      cAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAtt                                                                              >  1:195990/1‑139 (MQ=255)                                        gAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTc                                                                            >  2:198108/1‑139 (MQ=255)                                                      tGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGAc                                                                                                  <  2:484439/103‑1 (MQ=255)                                                      tGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGAc                                                                                                  >  1:484439/1‑103 (MQ=255)                                                      tGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCg                                                                  <  1:978667/135‑1 (MQ=255)                                                      tGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCg                                                                  >  2:978667/1‑135 (MQ=255)                                                      tGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGa                                                              >  1:931366/1‑139 (MQ=255)                                                       gTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                             <  2:195990/139‑1 (MQ=255)                                                       gTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                             <  2:483120/139‑1 (MQ=255)                                                       gTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                             <  2:931366/139‑1 (MQ=255)                                                           aaGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCAt                                                                          <  2:867361/122‑1 (MQ=255)                                                           aaGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCAt                                                                          >  1:867361/1‑122 (MQ=255)                                                             gCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTaa                                                       >  1:506054/1‑139 (MQ=255)                                                             gCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTaa                                                       >  2:549151/1‑139 (MQ=255)                                                             gCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTaa                                                       >  2:1174492/1‑139 (MQ=255)                                                                gAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGa                                                                                >  1:525570/1‑111 (MQ=255)                                                                gAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGa                                                                                <  2:525570/111‑1 (MQ=255)                                                                 aaaCTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                                   >  1:1675718/1‑139 (MQ=255)                                                                 aaaCTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGACGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                                   >  1:782741/1‑139 (MQ=255)                                                                  aaCTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGa                                                                                                >  2:729055/1‑93 (MQ=255)                                                                  aaCTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGa                                                                                                <  1:729055/93‑1 (MQ=255)                                                                  aaCTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                               <  1:3547/126‑1 (MQ=255)                                                                  aaCTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                               >  2:3547/1‑126 (MQ=255)                                                                  aaCTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGgag                                                  <  2:1675718/139‑1 (MQ=255)                                                                  aaCTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAACga                                                   >  1:1102244/1‑138 (MQ=255)                                                                  aaCTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAACga                                                   <  2:1102244/138‑1 (MQ=255)                                                                       aGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                               >  1:1159570/1‑105 (MQ=255)                                                                       aGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                               <  2:1159570/105‑1 (MQ=255)                                                                       aGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACt                                                         >  2:885440/1‑127 (MQ=255)                                                                       aGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACt                                                         <  1:885440/127‑1 (MQ=255)                                                                       aGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAgg                                                    <  1:1387620/132‑1 (MQ=255)                                                                       aGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAgg                                                    >  2:1387620/1‑132 (MQ=255)                                                                         a‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggct                                                                                                          >  1:1270106/1‑76 (MQ=255)                                                                         a‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggct                                                                                                          <  2:1270106/76‑1 (MQ=255)                                                                         a‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                                   <  1:1023194/131‑1 (MQ=255)                                                                         a‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                                   >  2:1023194/1‑131 (MQ=255)                                                                             tATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                                 <  1:1453630/82‑1 (MQ=255)                                                                             tATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                                 >  2:1453630/1‑82 (MQ=255)                                                                              aTGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATTCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACg                                       >  2:757994/1‑139 (MQ=255)                                                                              aTGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                             <  1:1042869/117‑1 (MQ=255)                                                                              aTGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                             >  2:1042869/1‑117 (MQ=255)                                                                               tGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACt                                                         <  2:177964/120‑1 (MQ=255)                                                                               tGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACt                                                         >  1:177964/1‑120 (MQ=255)                                                                               tGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGa                                      >  1:333050/1‑139 (MQ=255)                                                                                    gACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                                 <  1:198108/139‑1 (MQ=255)                                                                                    gACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                                 <  2:506054/139‑1 (MQ=255)                                                                                      cATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTAtt                               >  1:359891/1‑139 (MQ=255)                                                                                       aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                                 <  1:1325438/136‑1 (MQ=255)                                                                                       aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                                 >  2:1325438/1‑136 (MQ=255)                                                                                       aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTA‑‑CCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTccc                            <  2:676187/139‑1 (MQ=255)                                                                                        tCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATTCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTcc                             <  1:757994/139‑1 (MQ=255)                                                                                             aCCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTgg                        >  1:1372814/1‑139 (MQ=255)                                                                                              cccTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGt                       <  2:231889/139‑1 (MQ=255)                                                                                              cccTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGt                       <  2:359891/139‑1 (MQ=255)                                                                                                cTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTAt                     <  1:1174492/139‑1 (MQ=255)                                                                                                cTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTAt                     <  2:333050/139‑1 (MQ=255)                                                                                                cTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTAt                     >  2:350491/1‑139 (MQ=255)                                                                                                 tGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAt                                                                                                              >  2:389801/1‑49 (MQ=255)                                                                                                 tGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAt                                                                                                              <  1:389801/49‑1 (MQ=255)                                                                                                   ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGa                                                                                <  1:756280/77‑1 (MQ=255)                                                                                                   ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGa                                                                                >  2:756280/1‑77 (MQ=255)                                                                                                     tATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                            >  1:491916/1‑47 (MQ=255)                                                                                                     tATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                            <  2:491916/47‑1 (MQ=255)                                                                                                       tCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATc              <  2:1071134/139‑1 (MQ=255)                                                                                                          caTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATccac           <  1:359445/139‑1 (MQ=255)                                                                                                             cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGa                                                                                            >  1:202727/1‑55 (MQ=255)                                                                                                             cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGa                                                                                            <  2:202727/55‑1 (MQ=255)                                                                                                             cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                               <  1:1787734/84‑1 (MQ=255)                                                                                                             cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTg                                                               >  2:1787734/1‑84 (MQ=255)                                                                                                             cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                                 >  2:377070/1‑114 (MQ=255)                                                                                                             cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                                 <  1:377070/114‑1 (MQ=255)                                                                                                               tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAGTa      >  1:1381906/1‑139 (MQ=255)                                                                                                                 tttttAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAGTAGa    >  2:243908/1‑139 (MQ=255)                                                                                                                   tttAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAGTAGAAc  <  2:1372814/139‑1 (MQ=255)                                                                                                                       |                                                                                                                                       GTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGA‑CGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAGTAGAAC  >  NC_000913/3658776‑3659027

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