Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,366,935	C→T	*326* (TAG→TAA)	pspF ←	psp operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,366,935	0	C	T	98.7%	252.4 / NA	79	*326* (TAG→TAA)	pspF	psp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (40/38);  minor base G (0/1);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCT  >  NC_000913/1366807‑1367065                                                                                                                                 |                                                                                                                                   aGCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGt                                                                                                                          <  1:1587763/139‑1 (MQ=255)   cGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGc                                                                                                                        >  2:437504/1‑139 (MQ=255)           aGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTcaa                                                                                                                >  1:653330/1‑139 (MQ=255)                  cGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcg                                                                                                         >  1:1791391/1‑139 (MQ=255)                  cGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcg                                                                                                         >  1:1581303/1‑139 (MQ=255)                    cAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcgcg                                                                                                       <  2:653330/139‑1 (MQ=255)                            tCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACtgat                                                                                               >  2:290463/1‑139 (MQ=255)                            tCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACt                                                                                                  >  2:1484423/1‑136 (MQ=255)                            tCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACt                                                                                                  <  1:1484423/136‑1 (MQ=255)                             cTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACtgatg                                                                                              <  1:437504/139‑1 (MQ=255)                               gATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGc                                                                                                                        <  1:926633/111‑1 (MQ=255)                               gATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGc                                                                                                                        >  2:926633/1‑111 (MQ=255)                                 tgcgctgcgcTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATaa                                                                                          >  1:914124/1‑139 (MQ=255)                                    gctgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGtt                                                                                       <  2:1396493/139‑1 (MQ=255)                                    gctgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGtt                                                                                       <  1:290463/139‑1 (MQ=255)                                      tgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTcaaca                                                                                                              <  1:1418774/114‑1 (MQ=255)                                      tgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTcaaca                                                                                                              >  2:1418774/1‑114 (MQ=255)                                      tgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcgcg                                                                                                       <  2:581724/121‑1 (MQ=255)                                      tgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcgcg                                                                                                       >  1:581724/1‑121 (MQ=255)                                      tgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTaa                                                                                     <  2:914124/139‑1 (MQ=255)                                          ctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAAcccc                                                                                 <  2:575674/139‑1 (MQ=255)                                             aTCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGaa                                                                                                    >  1:1724654/1‑117 (MQ=255)                                             aTCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGaa                                                                                                    <  2:1724654/117‑1 (MQ=255)                                                  gCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTc                                                                         >  1:1389596/1‑139 (MQ=255)                                                     tGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGc                                                                      >  1:1190469/1‑139 (MQ=255)                                                       cAAACGACGTATTGATTATGTATGCCGAATAAGGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcg                                                                    >  1:1066365/1‑139 (MQ=255)                                                       cAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACt                                                                                                  <  1:225162/109‑1 (MQ=255)                                                       cAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACt                                                                                                  >  2:225162/1‑109 (MQ=255)                                                       cAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTgatgat                                                                                            <  2:1475783/115‑1 (MQ=255)                                                       cAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTgatgat                                                                                            >  1:1475783/1‑115 (MQ=255)                                                       cAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcg                                                                    >  1:127303/1‑139 (MQ=255)                                                         aaCGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcgcg                                                                  <  1:1788168/139‑1 (MQ=255)                                                              cgTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCtt                                                             <  2:127303/139‑1 (MQ=255)                                                                 aTTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTGTTTCAACAACGCGCGGAACTgatgat                                                                                            >  1:1098900/1‑105 (MQ=255)                                                                 aTTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTCCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTgatgat                                                                                            <  2:1098900/105‑1 (MQ=255)                                                                     aTTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTta                                                      <  2:1190469/139‑1 (MQ=255)                                                                      ttATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTgg                                                         <  1:1248115/135‑1 (MQ=255)                                                                      ttATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTgg                                                         >  2:1248115/1‑135 (MQ=255)                                                                         tgtatgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATAtt                                                  >  1:884138/1‑139 (MQ=255)                                                                             tgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTgg                                                         >  2:17341/1‑128 (MQ=255)                                                                             tgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTgg                                                         <  1:17341/128‑1 (MQ=255)                                                                             tgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTatat                                                   >  1:1772338/1‑134 (MQ=255)                                                                             tgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTatat                                                   <  2:1772338/134‑1 (MQ=255)                                                                             tgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTccc                                              >  1:1700247/1‑139 (MQ=255)                                                                              gCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCt                                             >  2:3767/1‑139 (MQ=255)                                                                               ccGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcg                                                                                                         <  1:1760586/78‑1 (MQ=255)                                                                               ccGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcg                                                                                                         >  2:1760586/1‑78 (MQ=255)                                                                               ccGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcgcg                                                                  <  2:1107861/117‑1 (MQ=255)                                                                               ccGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCtt                                                             >  2:1580162/1‑122 (MQ=255)                                                                               ccGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCtt                                                             <  1:1580162/122‑1 (MQ=255)                                                                               ccGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTCAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcgcg                                                                  >  1:1107861/1‑117 (MQ=255)                                                                                   aTAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCt                                                                                                                             <  1:1570587/54‑1 (MQ=255)                                                                                   aTAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCt                                                                                                                             >  2:1570587/1‑54 (MQ=255)                                                                                   aTAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACtg                                                                                                 <  2:803820/82‑1 (MQ=255)                                                                                   aTAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACtg                                                                                                 >  1:803820/1‑82 (MQ=255)                                                                                   aTAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcgcg                                                                  <  1:1262677/113‑1 (MQ=255)                                                                                   aTAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcgcg                                                                  >  2:1262677/1‑113 (MQ=255)                                                                                       aGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAAc                                                                                                   <  1:1713421/76‑1 (MQ=255)                                                                                       aGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAAc                                                                                                   >  2:1713421/1‑76 (MQ=255)                                                                                         cATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGtt                                                                                       >  2:111553/1‑86 (MQ=255)                                                                                         cATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGtt                                                                                       <  1:111553/86‑1 (MQ=255)                                                                                            tCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcg                                                                    >  2:1331282/1‑102 (MQ=255)                                                                                            tCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcg                                                                    <  1:1331282/102‑1 (MQ=255)                                                                                              aCGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcg                                                                                                         >  1:1211588/1‑63 (MQ=255)                                                                                              aCGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcg                                                                                                         <  2:1211588/63‑1 (MQ=255)                                                                                              aCGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGtt                             <  2:1791391/139‑1 (MQ=255)                                                                                               cgccgcATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTg                            >  1:649198/1‑139 (MQ=255)                                                                                                gccgcATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTgc                           >  1:1110872/1‑139 (MQ=255)                                                                                                    cATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTgatgat                                                                                            >  1:766266/1‑70 (MQ=255)                                                                                                    cATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTgatgat                                                                                            <  2:766266/70‑1 (MQ=255)                                                                                                    cATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTgcagca                       >  2:1234562/1‑139 (MQ=255)                                                                                                         ggCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAActct                  <  2:649198/139‑1 (MQ=255)                                                                                                               ttGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTcc            <  2:1700247/139‑1 (MQ=255)                                                                                                                   aTTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCctgct        <  1:3767/139‑1 (MQ=255)                                                                                                                   aTTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCctgct        <  2:1110872/139‑1 (MQ=255)                                                                                                                       cTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAAcccc                                                                                 >  2:298499/1‑62 (MQ=255)                                                                                                                       cTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAAcccc                                                                                 <  1:298499/62‑1 (MQ=255)                                                                                                                        tCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATc   >  1:966080/1‑139 (MQ=255)                                                                                                                         cAACTTCGGTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCTCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCt  <  2:1066365/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                   AGCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCT  >  NC_000913/1366807‑1367065

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