Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,018,260	(A)8→7	100%	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	.	94.9%	329.5 / 8.1	79	coding (308/444 nt)	fliJ	flagellar protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base . (36/39);  total (37/42)
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                                                                                                    |                                                                                                agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  1:2362247/1‑100 (MQ=255)  gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  2:2054945/1‑100 (MQ=255)  gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                <  1:1748933/100‑1 (MQ=255)     gTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATAAGTGGACGCAGAAAGGTGACATTGACCTGAACAGTTGGCGAGaaaaa                                                                                             >  2:1758754/1‑100 (MQ=255)     gTTTATCCAGACGCTGGAAAAAGCCAGTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaa                                                                                             >  1:1580462/1‑100 (MQ=255)       ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           >  2:2240610/1‑100 (MQ=255)           ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                       <  2:242830/100‑5 (MQ=255)             aGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAaaaaaaaaaagt                                                                                     >  1:479023/1‑98 (MQ=255)                  cTGGGAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               >  2:1502960/1‑100 (MQ=255)                  cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               <  2:2186374/100‑1 (MQ=255)                  cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               >  1:1899909/1‑100 (MQ=255)                    ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAgg                                                                             >  1:801636/1‑100 (MQ=255)                     gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               <  2:451413/97‑1 (MQ=255)                     gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGc                                                                            >  1:1713051/1‑100 (MQ=255)                       aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCt                                                                          >  1:1455893/1‑100 (MQ=255)                       aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCt                                                                          >  2:1319028/1‑100 (MQ=255)                       aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCt                                                                          >  1:2427011/1‑100 (MQ=255)                       aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCt                                                                          <  1:2240623/100‑1 (MQ=255)                       aaaaGCCATTACTCAGCATCGCCAGCAACTGAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACCGTTGGCGAGAAAAAAAACAACGTTTGCAGGcc                                                                           <  1:2249535/100‑1 (MQ=255)                        aaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTg                                                                         <  1:2054956/100‑1 (MQ=255)                        aaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTg                                                                         <  1:2259891/100‑1 (MQ=255)                          aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  1:1582037/100‑1 (MQ=255)                          aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  2:2422135/100‑1 (MQ=255)                          aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  1:1247778/100‑1 (MQ=255)                           gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCa                                                                      <  1:653275/100‑1 (MQ=255)                               ttACTCAGATTCGCTAGCACCTTAATCAGTCGACGGCGAGCGTTGAAATTACGCTGAACAGTGGGTGAG‑AAAAAAACAACGTTTGCAGCCCGGGCa                                                                      <  2:1126707/96‑1 (MQ=255)                                   tCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  1:2282200/100‑1 (MQ=255)                                    cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCGGAACAGTGGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  2:1105827/99‑1 (MQ=255)                                     aGCATCGCCAGCAACTTAATCAGTGTACGCAGAAAGTTGACAGTGCCCTGAACAGTTGGCTAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  2:875294/100‑1 (MQ=255)                                     aGCATCGACAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTAAAAAAAAAAAAA‑AAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  1:2346417/100‑1 (MQ=255)                                      gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                           >  2:315842/1‑100 (MQ=255)                                      gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                           >  1:559903/1‑100 (MQ=255)                                         tCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAAc                                                        >  1:876224/1‑100 (MQ=255)                                         tCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAAc                                                        <  1:470008/100‑1 (MQ=255)                                          cGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACg                                                       >  1:1050615/1‑100 (MQ=255)                                           gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      >  1:689051/1‑100 (MQ=255)                                           gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      <  2:152248/100‑1 (MQ=255)                                            ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:1153205/100‑1 (MQ=255)                                            ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:19940/100‑1 (MQ=255)                                            ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:2064228/100‑1 (MQ=255)                                             cAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCGCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCa                                                    >  2:1606574/1‑100 (MQ=255)                                              aGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCaa                                                   <  1:766910/100‑1 (MQ=255)                                                cAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATc                                                 <  1:1232719/100‑1 (MQ=255)                                                  aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAACCCa                                               >  1:964329/1‑100 (MQ=255)                                                   cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              >  1:1693650/1‑100 (MQ=255)                                                    ttAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACCACGTTTGCAGTCCTGGCAGACACTGCAGGAACGGCAATTCAcg                                             >  1:378425/1‑100 (MQ=255)                                                     tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            >  2:1064221/1‑100 (MQ=255)                                                      aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           >  1:2378826/1‑100 (MQ=255)                                                        tCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                         >  2:425570/1‑100 (MQ=255)                                                         cAGTGGACGCAGAAAGTGGACATCGCCCGGAACAGTTGACGAG‑AAAAAAACAACGTTTGCAGGCCTAGCAGACACTGCAGGAACGGCAATCCAcggc                                           <  2:1009797/97‑1 (MQ=255)                                                             ggACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTg                                    >  1:2031776/1‑100 (MQ=255)                                                             ggACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTg                                    >  2:1717228/1‑100 (MQ=255)                                                             ggACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTCCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                     >  2:1554949/1‑99 (MQ=255)                                                              gACGCAGAAAGTTGACATTGCCCTAAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                     <  2:834534/98‑1 (MQ=255)                                                                ctcAGAAAGTTGACATCACCCTGAACAGTTGCCGAG‑AAAAAACCAACGTTTGCAGGCCTGGCAGGCGCTGCAGGAACGGCAATCCACGGCCGCACTGCtt                                 >  2:307842/3‑100 (MQ=255)                                                                cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACAGCt                                  >  2:260141/1‑99 (MQ=255)                                                                 gCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTg                                <  1:1221087/100‑1 (MQ=255)                                                                 gCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTg                                >  1:2513824/1‑100 (MQ=255)                                                                   aGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              <  1:1731888/100‑1 (MQ=255)                                                                   aGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              >  2:2512195/1‑100 (MQ=255)                                                                      aaGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaa                           <  1:662444/100‑1 (MQ=255)                                                                       aGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGCCACTGCAGGGACGGCAATCCACGGCGGCACTGCTTGCAGaaa                          >  2:1823892/1‑100 (MQ=255)                                                                       aGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaa                          >  1:458043/1‑100 (MQ=255)                                                                         ttGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAc                        >  1:1353166/1‑100 (MQ=255)                                                                            aCATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGc                     <  1:1341638/100‑1 (MQ=255)                                                                             cATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                    >  1:169620/1‑100 (MQ=255)                                                                              aTTGCCCTGAAGAGATGGCGA‑CAAAAGAACAACGTTTGCAGGCCTGGCACACAATGCAGGAACGTCATTCTATGTCGGCACTGCGGGCAGAAAACCg                      <  2:1050993/97‑1 (MQ=255)                                                                               ttGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTc                  <  1:497128/100‑1 (MQ=255)                                                                                 gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa                >  1:693014/1‑100 (MQ=255)                                                                                 gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa                >  1:2078388/1‑100 (MQ=255)                                                                                  cccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               <  2:2439862/100‑1 (MQ=255)                                                                                  cccTGAACAGTTGGCGAG‑‑AAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:790936/100‑1 (MQ=255)                                                                                  cccTGAACAGTTGGCCAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               <  1:2564331/100‑1 (MQ=255)                                                                                   ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              >  2:898079/1‑100 (MQ=255)                                                                                   ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  2:1867014/100‑1 (MQ=255)                                                                                   ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:1402933/100‑1 (MQ=255)                                                                                    cTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             <  1:2068107/100‑1 (MQ=255)                                                                                     tGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAg            <  1:1829400/100‑1 (MQ=255)                                                                                      gAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa           <  2:842638/100‑1 (MQ=255)                                                                                      gAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa           >  1:1345722/1‑100 (MQ=255)                                                                                         cAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaa        <  1:70902/100‑1 (MQ=255)                                                                                           gTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg      <  2:2426998/100‑1 (MQ=255)                                                                                           gTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACAGCAATCCACGGCGGCACTGCTTGCAGCAAACCGCCACGATCAGAAAAAg      >  2:2388235/1‑100 (MQ=255)                                                                                             tGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat    >  2:2421570/1‑100 (MQ=255)                                                                                              ggCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   <  2:559900/100‑1 (MQ=255)                                                                                               gCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgattg  <  1:1769791/100‑3 (MQ=255)                                                                                                    |                                                                                                AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348

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