Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,018,260	(A)8→7	coding (315/444 nt)	fliJ →	flagellar protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,018,253	0	A	.	97.4%	342.2 / 2.0	77	coding (308/444 nt)	fliJ	flagellar protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (31/44);  total (32/45)
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

AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348                                                                                                    |                                                                                                agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAACCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  2:638954/1‑100 (MQ=255)  gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  2:629543/1‑100 (MQ=255)     gTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaa                                                                                             >  2:865210/1‑100 (MQ=255)     gTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaa                                                                                             >  1:35637/1‑100 (MQ=255)      tttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaa                                                                                            >  2:1896593/1‑100 (MQ=255)      tttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaa                                                                                            >  2:775153/1‑100 (MQ=255)       ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           >  1:1423185/1‑100 (MQ=255)       ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           >  2:2229144/1‑100 (MQ=255)       ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           >  1:359545/1‑100 (MQ=255)        tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgagatcggaag                                                                                          >  1:2238901/1‑93 (MQ=255)        tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                          <  1:2576513/100‑2 (MQ=255)        tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                          >  1:2601499/1‑99 (MQ=255)        tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                          >  1:2670509/1‑99 (MQ=255)         aTCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaaca                                                                                         >  2:313753/1‑98 (MQ=255)           ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                       <  2:2463542/100‑5 (MQ=255)            cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACg                                                                                     >  1:604549/1‑100 (MQ=255)            cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACg                                                                                     >  1:1776920/1‑100 (MQ=255)               aCGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGttt                                                                                  <  1:722674/100‑1 (MQ=255)                cGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGAGAG‑AAAAAAACAACGTTTg                                                                                 <  2:876509/100‑1 (MQ=255)                 gCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGc                                                                                <  2:292116/100‑1 (MQ=255)                  cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               >  1:180640/1‑100 (MQ=255)                  cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               <  1:844875/100‑1 (MQ=255)                   tGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               <  2:859508/99‑1 (MQ=255)                    ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAgg                                                                             >  1:2157592/1‑100 (MQ=255)                      aaaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGcc                                                                           <  2:1265446/100‑1 (MQ=255)                       aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCt                                                                          >  1:800164/1‑100 (MQ=255)                        aaaGCCATTACTCATCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAAAAAAAAA‑AAAAAAAACAACGTTTGCAGGCCTg                                                                         <  2:1881076/100‑1 (MQ=255)                        aaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTg                                                                         <  1:2341007/100‑1 (MQ=255)                         aaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTgg                                                                        >  2:58690/1‑100 (MQ=255)                          aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  1:2511306/100‑1 (MQ=255)                           gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCa                                                                      >  1:330946/1‑100 (MQ=255)                            ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAg                                                                     >  2:940826/1‑100 (MQ=255)                            ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAg                                                                     >  1:2215237/1‑100 (MQ=255)                             cATTGCTCAGCATCGCCAGCAACGTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAg                                                                     <  2:1587168/99‑1 (MQ=255)                                tACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                                 >  2:743285/1‑100 (MQ=255)                                  cTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTg                                                               >  2:1717063/1‑100 (MQ=255)                                   tCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  1:48673/100‑1 (MQ=255)                                    cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             >  1:1894928/1‑100 (MQ=255)                                    cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             <  2:2594576/100‑1 (MQ=255)                                     aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  1:876950/100‑1 (MQ=255)                                     aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            >  2:1830459/1‑100 (MQ=255)                                        aTCGCCAGCAACTTAATAAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  2:877239/95‑1 (MQ=255)                                           gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      >  1:2514330/1‑100 (MQ=255)                                           gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      <  1:644477/100‑1 (MQ=255)                                            ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  1:1844719/100‑1 (MQ=255)                                            ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  1:1580863/100‑1 (MQ=255)                                             cAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCa                                                    <  2:2601480/100‑1 (MQ=255)                                              aGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCaa                                                   <  1:629548/100‑1 (MQ=255)                                               gCAACTTAACCAGTGGACGCAGAAAGTTGACACTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAAt                                                  <  1:985998/100‑1 (MQ=255)                                                  aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                               >  1:771254/1‑100 (MQ=255)                                                  aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                               <  1:58691/100‑1 (MQ=255)                                                   cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              <  1:1177325/100‑1 (MQ=255)                                                   cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              <  1:2468635/100‑1 (MQ=255)                                                   cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              <  1:743290/100‑1 (MQ=255)                                                     tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            >  2:1393680/1‑100 (MQ=255)                                                      aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcg                                             >  2:985338/1‑98 (MQ=255)                                                        tCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                         >  1:2746144/1‑100 (MQ=255)                                                         cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                        <  2:1331521/100‑1 (MQ=255)                                                          aGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCa                                       >  1:2295460/1‑100 (MQ=255)                                                            tGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                     <  1:871360/100‑1 (MQ=255)                                                              gACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGc                                   <  2:180637/100‑1 (MQ=255)                                                               aCGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCt                                  <  2:986988/100‑1 (MQ=255)                                                                 gCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTg                                >  2:1254831/1‑100 (MQ=255)                                                                  cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  2:2204411/100‑1 (MQ=255)                                                                  cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  1:313756/100‑1 (MQ=255)                                                                   aGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              >  2:1014088/1‑100 (MQ=255)                                                                   aGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              <  2:1958724/100‑1 (MQ=255)                                                                   aGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              <  1:298353/100‑1 (MQ=255)                                                                     aaaGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGa                            <  1:940832/100‑1 (MQ=255)                                                                        gTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                         <  1:2229162/100‑1 (MQ=255)                                                                          tGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcc                       >  2:2024763/1‑100 (MQ=255)                                                                              aTTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCt                   <  2:2140908/100‑1 (MQ=255)                                                                               ttGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTc                  >  1:2097859/1‑100 (MQ=255)                                                                                 gCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg                 >  1:877721/1‑100 (MQ=255)                                                                                 gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa                >  1:1295501/1‑100 (MQ=255)                                                                                 gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa                >  1:1296544/1‑100 (MQ=255)                                                                                  cccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               >  1:1099929/1‑100 (MQ=255)                                                                                   ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  2:1795287/100‑1 (MQ=255)                                                                                   ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:1830474/100‑1 (MQ=255)                                                                                   catGAACAGTTGGCGAG‑AAAAAAACAACGTCTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  2:702316/98‑1 (MQ=255)                                                                                     tGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAg            <  1:985345/100‑1 (MQ=255)                                                                                     tGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAg            <  2:1045635/100‑1 (MQ=255)                                                                                      gAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa           <  2:981747/100‑1 (MQ=255)                                                                                           gTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg      >  1:1813550/1‑100 (MQ=255)                                                                                           gTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg      >  2:177937/1‑100 (MQ=255)                                                                                           gTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCGCGATCAGAAAAAg      >  1:1789566/1‑100 (MQ=255)                                                                                            ttGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAga     <  2:1396241/100‑1 (MQ=255)                                                                                              ggCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   >  1:1720634/1‑100 (MQ=255)                                                                                              ggCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   <  1:1896608/100‑1 (MQ=255)                                                                                              ggCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   >  2:1673249/1‑100 (MQ=255)                                                                                               gCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg  <  1:181865/100‑1 (MQ=255)                                                                                               gCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg  <  2:2274236/100‑1 (MQ=255)                                                                                                    |                                                                                                AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348

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