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	.	98.9%	405.9 / ‑2.2	88	coding (308/444 nt)	fliJ	flagellar protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (43/44);  total (43/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 = 9.91e-01

GTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATG  >  NC_000913/2018158‑2018347                                                                                                |                                                                                               gTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaa                                                                                            <  2:2765614/100‑1 (MQ=255)    tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                         >  2:318136/1‑99 (MQ=255)    tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                         >  1:2414570/1‑99 (MQ=255)    tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                         <  2:2715591/100‑2 (MQ=255)    tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                         <  2:3332977/100‑2 (MQ=255)     aTCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGAAAGaaaaaaac                                                                                         <  2:1961141/99‑2 (MQ=255)       ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                      >  2:636613/1‑96 (MQ=255)            cGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTg                                                                                >  2:1250744/1‑100 (MQ=255)            cGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTg                                                                                <  1:1420009/100‑1 (MQ=255)            cGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTg                                                                                <  1:1446502/100‑1 (MQ=255)            cGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTg                                                                                <  2:324606/100‑1 (MQ=255)              cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                              >  2:317663/1‑100 (MQ=255)                ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAgg                                                                            <  1:2140940/100‑1 (MQ=255)                 gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAgg                                                                            <  2:1536560/99‑1 (MQ=255)                 gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGc                                                                           <  2:3183073/100‑1 (MQ=255)                 gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGc                                                                           >  1:489258/1‑100 (MQ=255)                 gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGACCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGc                                                                           >  2:1856723/1‑100 (MQ=255)                 gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCACTGGACGAAGAAAGTTGACATTGCCCTGAAAAAAAAAAAA‑AAAAAAAACAACGTTTGCAGGc                                                                           <  2:1458568/100‑1 (MQ=255)                   aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCt                                                                         <  1:2876396/100‑1 (MQ=255)                    aaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTg                                                                        <  2:199481/100‑1 (MQ=255)                      aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGc                                                                      <  2:1157741/100‑1 (MQ=255)                       gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCa                                                                     <  2:3364303/100‑1 (MQ=255)                       gCCACTACTCAGCATCGCGAGCAACTTAATCAATGGACTCAGAAAGTTGACATTGCCCTGAACAGTTGTGGAG‑AATTAAACAACGTTTGCAGGCCTGTCa                                                                     <  2:2659250/100‑1 (MQ=255)                        ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAg                                                                    >  1:2227171/1‑100 (MQ=255)                        ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAg                                                                    <  2:1196451/100‑1 (MQ=255)                         cATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGa                                                                   >  1:332835/1‑100 (MQ=255)                         cATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGa                                                                   >  1:672941/1‑100 (MQ=255)                         cATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGa                                                                   >  1:844593/1‑100 (MQ=255)                          aTTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                  <  1:1687235/100‑1 (MQ=255)                          aTTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                  <  2:2603440/100‑1 (MQ=255)                             aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                               >  1:2652238/1‑100 (MQ=255)                                cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                             <  2:1380655/99‑1 (MQ=255)                                cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                            >  2:2789763/1‑100 (MQ=255)                                  gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                          >  1:1168444/1‑100 (MQ=255)                                   cATCGCCAGCAACTTAATCAGTTTACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGa                                                         <  1:2839384/100‑1 (MQ=255)                                   cATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                          <  2:2778335/99‑1 (MQ=255)                                       gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                     <  2:844593/100‑1 (MQ=255)                                        ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                    <  1:3403379/100‑1 (MQ=255)                                        ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                    <  1:256634/100‑1 (MQ=255)                                        ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                    <  1:2381539/100‑1 (MQ=255)                                         cAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCa                                                   <  1:455335/100‑1 (MQ=255)                                         cAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCa                                                   >  1:2384809/1‑100 (MQ=255)                                          aGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCaa                                                  >  2:2539603/1‑100 (MQ=255)                                           gCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAAt                                                 >  1:608814/1‑100 (MQ=255)                                             aaCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATcc                                               >  1:117507/1‑100 (MQ=255)                                              aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                              >  1:295497/1‑100 (MQ=255)                                               cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                             <  1:2036139/100‑1 (MQ=255)                                                     cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                       >  1:3026357/1‑100 (MQ=255)                                                         ggACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCAAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTg                                   >  1:2403179/1‑100 (MQ=255)                                                         ggACGCAGAAAGTTGACATTGCCCTGAACAGTAGGAGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                    <  2:547027/99‑1 (MQ=255)                                                          gACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGc                                  >  1:1304506/1‑100 (MQ=255)                                                           aCGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCt                                 >  2:1704049/1‑100 (MQ=255)                                                           aCGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCt                                 <  2:870905/100‑1 (MQ=255)                                                            cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCtt                                >  2:2359324/1‑100 (MQ=255)                                                            cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCtt                                >  1:261315/1‑100 (MQ=255)                                                             gCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTg                               <  2:295496/100‑1 (MQ=255)                                                              cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                              <  2:1810326/100‑1 (MQ=255)                                                              cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                              >  2:533712/1‑100 (MQ=255)                                                              cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                              >  1:2820047/1‑100 (MQ=255)                                                              cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                              >  2:2898252/1‑100 (MQ=255)                                                               aGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                             <  2:962035/100‑1 (MQ=255)                                                               aGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                             <  2:2227165/100‑1 (MQ=255)                                                                 aaaGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGa                           <  1:2742509/100‑1 (MQ=255)                                                                  aaGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaa                          >  2:2679884/1‑100 (MQ=255)                                                                    gTTGACAGTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                        >  1:1794840/1‑100 (MQ=255)                                                                      tGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcc                      >  1:185699/1‑100 (MQ=255)                                                                       gACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCg                     >  1:1155323/1‑100 (MQ=255)                                                                       gACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCg                     >  1:1199067/1‑100 (MQ=255)                                                                          aTTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCt                  >  2:930092/1‑100 (MQ=255)                                                                            tGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg                <  1:2359331/100‑1 (MQ=255)                                                                            tGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg                <  1:317664/100‑1 (MQ=255)                                                                            tGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCAGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg                <  2:2384802/100‑1 (MQ=255)                                                                             gCCCTGAGCAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa               >  1:2419075/1‑100 (MQ=255)                                                                             gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa               >  1:2184300/1‑100 (MQ=255)                                                                              cccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt              >  2:407746/1‑100 (MQ=255)                                                                               ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt              <  2:2743662/99‑1 (MQ=255)                                                                               ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc             <  2:2652229/100‑1 (MQ=255)                                                                                cTGAACAGTTGGCGAG‑AAAATAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa            >  2:2081850/1‑100 (MQ=255)                                                                                cTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa            >  1:1378880/1‑100 (MQ=255)                                                                                 tGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAg           <  1:318137/100‑1 (MQ=255)                                                                                 tGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAg           >  2:3365765/1‑100 (MQ=255)                                                                                 tGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAg           >  2:2908399/1‑100 (MQ=255)                                                                                  gAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa          <  2:2653446/100‑1 (MQ=255)                                                                                  gAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa          >  2:903551/1‑100 (MQ=255)                                                                                      aGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaaa      <  2:332834/100‑1 (MQ=255)                                                                                       gTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg     <  1:1937115/100‑1 (MQ=255)                                                                                       gTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg     >  2:1266998/1‑100 (MQ=255)                                                                                        ttGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAga    <  2:71976/100‑1 (MQ=255)                                                                                        ttGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAga    <  1:2539611/100‑1 (MQ=255)                                                                                        ttGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCAGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAga    <  2:1818451/100‑1 (MQ=255)                                                                                         tGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat   <  1:930093/100‑1 (MQ=255)                                                                                         tGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat   >  2:2134144/1‑100 (MQ=255)                                                                                          ggCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg  >  2:3392005/1‑100 (MQ=255)                                                                                          ggCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg  <  1:1856727/100‑1 (MQ=255)                                                                                          ggCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg  >  1:3404019/1‑100 (MQ=255)                                                                                                |                                                                                               GTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATG  >  NC_000913/2018158‑2018347

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: ‑