Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	136,530	C→T	6.2%	intergenic (‑66/+40)	speE ← / ← yacC	spermidine synthase (putrescine aminopropyltransferase)/PulS_OutS family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	136,530	0	C	T	6.2%	293.9 / 4.8	96	intergenic (‑66/+40)	speE/yacC	spermidine synthase (putrescine aminopropyltransferase)/PulS_OutS family protein
Reads supporting (aligned to +/- strand):  ref base C (57/33);  new base T (3/3);  total (60/36)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

ATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACC  >  NC_000913/136384‑136674                                                                                                                                                   |                                                                                                                                                 aTGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTGGGCAATGGGTTGTAACCACGTTGTTTAACACCCGAAAAAAAAGGGGCCAAAAAAATAGATAAATTTAACCCCggt                                                                                                                                               >  1:1275597/1‑149 (MQ=255)       cAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACACGTAAAAAAAGTGCGCAACATAATAGCTAACATTGACCGCGGATGCACg                                                                                                                                         >  2:1116777/1‑150 (MQ=255)          aaCGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTc                                                                                                                                      <  2:931651/150‑1 (MQ=255)          aaCGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTc                                                                                                                                      <  1:362154/150‑1 (MQ=255)           aCGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCa                                                                                                                                     <  2:237612/150‑1 (MQ=255)              ttATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAatg                                                                                                                                  <  2:1216296/150‑2 (MQ=255)               tATCTACCGCACAGTACTGCCGAAACTGGTCGTGTCGCGTTTCATGCCAATTTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGCCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAatat                                                                                                                                 <  2:260627/150‑1 (MQ=255)               tATCTACCGCAAAGTAATGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAatgt                                                                                                                                 <  1:501453/150‑3 (MQ=255)                aTCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATAtt                                                                                                                                <  2:811303/150‑1 (MQ=255)                aTCTACCGCAAAGTACTGACCAAACTGGTCGTGTAGCGTTTCATTCCACAGTTTTTTTTCGGCCAAGGGTTGATACCTACTTTGTTAAAACCCGGAAAAAAAGGGCGCAACATAATAGCAAACATTTACAGCTGATGCACGGTCAataat                                                                                                                                >  2:545198/1‑148 (MQ=255)                  cTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTTTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAAAATTGACCGCGGATGCACGGTAAATAttta                                                                                                                              >  1:1622400/1‑149 (MQ=255)                  atACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGTTCAATATTTc                                                                                                                              <  1:545198/149‑1 (MQ=255)                   tACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCt                                                                                                                             <  2:1059239/150‑1 (MQ=255)                    gccGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTa                                                                                                                            >  1:379789/2‑150 (MQ=255)                    aCCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAATATTGACCGTGCATCCGCGGTCAATGttagct                                                                                                                            <  2:1206818/150‑5 (MQ=255)                    aCCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTa                                                                                                                            >  2:649803/1‑150 (MQ=255)                    aCCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTAGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTa                                                                                                                            >  2:1288054/1‑150 (MQ=255)                     ccGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTAc                                                                                                                           >  1:862108/1‑150 (MQ=255)                       gCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCCTGCCACTGTTTTTTTTCGGCCATGGGTTGATAACTCCTTTGTTAACACCCGTAAAAAAAGGGCGAAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaa                                                                                                                         >  2:1282238/1‑150 (MQ=255)                       gCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATTGCTAACATTGACCGCGGATGCACGGGCAATTTTTCTACaa                                                                                                                         >  2:1671610/1‑150 (MQ=255)                         aaaGTACTGCCCAAACTGGTCGTGTAGCTTTTCATGCCACTGTTTTTTTTCGTCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCTGATGCACGGTCAATATTTCTACaaaa                                                                                                                       <  1:1269193/150‑1 (MQ=255)                            ttacttaccAAACTGGTCGTGTATCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAgg                                                                                                                    <  2:1622400/143‑1 (MQ=255)                              aCTGACCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTa                                                                                                                  <  1:1329651/150‑1 (MQ=255)                               cTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTAt                                                                                                                 <  2:1407362/150‑1 (MQ=255)                                 gCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACAATTACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTAtaa                                                                                                               >  2:593194/1‑148 (MQ=255)                                  cccAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCATAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAg                                                                                                              >  1:418142/1‑150 (MQ=255)                                   ccAAACTGGTCGTGTAGCGTTTCATGCCACTGGTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAgg                                                                                                             <  2:302881/150‑1 (MQ=255)                                      aaCTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAATAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTTCAAAAAGGTATCAGGGat                                                                                                          <  1:1516165/150‑1 (MQ=255)                                        cTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGatta                                                                                                        <  2:1011431/150‑1 (MQ=255)                                          ggTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGattatt                                                                                                      >  2:1362052/1‑150 (MQ=255)                                          ggTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGattatt                                                                                                      >  2:516605/1‑150 (MQ=255)                                          ggTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACAACCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGattatt                                                                                                      >  2:396867/1‑150 (MQ=255)                                            tCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGTGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCTCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTg                                                                                                    <  1:1528971/150‑1 (MQ=255)                                            gcGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTGGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGATCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTg                                                                                                    <  2:446221/149‑1 (MQ=255)                                              gtgtAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGAACGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGAc                                                                                                  >  1:1492273/1‑150 (MQ=255)                                                    cGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGc                                                                                            >  1:710260/1‑150 (MQ=255)                                                          aTGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCCCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAg                                                                                      >  2:422706/1‑150 (MQ=255)                                                           tGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAgg                                                                                     <  2:465534/150‑1 (MQ=255)                                                             ccACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCt                                                                                   <  1:1446134/150‑1 (MQ=255)                                                             ccACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCt                                                                                   <  2:1653638/150‑1 (MQ=255)                                                               aCTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTa                                                                                 >  2:780079/1‑150 (MQ=255)                                                                cTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTaa                                                                                >  1:855604/1‑150 (MQ=255)                                                                   ttttttttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGa                                                                             >  2:1299160/1‑150 (MQ=255)                                                                   ttttttttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGGCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCCGGATTAAGGa                                                                             >  2:852882/1‑150 (MQ=255)                                                                      tttttCGGCCATGGGTTGATACCTCCTTTGTTACCACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCCCGGTCAATATTTCTACAAAAAGGTATCACGGATTATTTGACGTAGGCAATCAGGCTTAAGGAATc                                                                          >  1:1469193/1‑150 (MQ=255)                                                                      tttttCGGACATGGGTTGATACCTCCTTTGTTAACACCAGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATc                                                                          <  1:99941/150‑1 (MQ=255)                                                                       ttttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATcg                                                                         >  2:629973/1‑150 (MQ=255)                                                                        tttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATcgc                                                                        <  1:649803/150‑1 (MQ=255)                                                                        tttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATcgc                                                                        >  1:743379/1‑150 (MQ=255)                                                                          tCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCgg                                                                      >  1:1608488/1‑150 (MQ=255)                                                                          tCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAACGTCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCgg                                                                      <  1:482139/150‑1 (MQ=255)                                                                               cATGGGTTGATTCCTCCTTTGTTAACACCCGT‑AAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGCCTTAAGGAATCGCGGGCCAgg                                                                >  1:709363/1‑150 (MQ=255)                                                                               cATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGCCCGCGGATGCACGGTCAATATTTCTACATAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGATTCGCGGGCCAg                                                                 >  1:822386/1‑150 (MQ=255)                                                                                 tGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAggg                                                               <  1:1014552/150‑1 (MQ=255)                                                                                  gggTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGTGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGGCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGc                                                              <  1:1085754/150‑1 (MQ=255)                                                                                   ggTTGATACCTCCTTTGTTAACCCCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGG‑CAGGGCtt                                                            >  1:665140/1‑150 (MQ=255)                                                                                    gTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCtt                                                            >  2:1553574/1‑150 (MQ=255)                                                                                      tGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGACAATATTTCTAC‑AAAAGGTATAAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGATTTGc                                                         >  2:1649702/1‑150 (MQ=255)                                                                                      tGATACCTCCTTTGATAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTg                                                          <  2:1192472/150‑1 (MQ=255)                                                                                      tGATACCCCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACAAAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATATGACGTAGGCAAGCAGGCTTAAGGACTCGCGGGCCAGGGCTATg                                                          >  2:776374/1‑150 (MQ=255)                                                                                       gATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGc                                                         <  2:862108/150‑1 (MQ=255)                                                                                          aCCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtt                                                      >  2:978813/1‑150 (MQ=255)                                                                                          aCCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtt                                                      >  2:191232/1‑150 (MQ=255)                                                                                           cctcctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAttt                                                     >  2:1208967/1‑150 (MQ=255)                                                                                              cctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtttttt                                                  <  2:433488/150‑1 (MQ=255)                                                                                              cctTTGTTAACACACGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtttttt                                                  <  1:676331/150‑1 (MQ=255)                                                                                                  tGTTAACACCCGTAAAAAAAGGGCGCAACATGATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCg                                              >  1:1587916/1‑150 (MQ=255)                                                                                                    ttAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCgcc                                            >  1:1549614/1‑148 (MQ=255)                                                                                                     tAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACg                                           >  2:200548/1‑150 (MQ=255)                                                                                                      aaCACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACgg                                          >  2:42179/1‑150 (MQ=255)                                                                                                       acacCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACggg                                         >  1:1340223/1‑150 (MQ=255)                                                                                                          cccGTAAAAAAAGGGCGCATCATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCACGGCTTTGCATTTTTTAGCGACGGGaat                                      >  2:1193713/1‑150 (MQ=255)                                                                                                            cGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgc                                    <  2:1305020/150‑1 (MQ=255)                                                                                                                aaaaaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaat                                >  1:521780/1‑150 (MQ=255)                                                                                                                aaaaaaGGGCGCAACATAATAGCTAACATTGACAGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaat                                >  1:1595798/1‑150 (MQ=255)                                                                                                                  aaaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGTCAAGCAGGCTTAAGGAATCGCTGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaatgc                              >  1:256337/1‑150 (MQ=255)                                                                                                                   aaaGGGCGTAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaatgcc                             >  1:946182/1‑150 (MQ=255)                                                                                                                   aaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaatgcc                             <  1:629973/150‑1 (MQ=255)                                                                                                                     aGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTGTGCATTTTTTAGCGACGGGAATGCCAATGCCGc                           >  1:1233555/1‑150 (MQ=255)                                                                                                                       ggCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGTCGCTg                         >  2:657439/1‑150 (MQ=255)                                                                                                                               ataataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATcgcg                 <  2:1275597/150‑1 (MQ=255)                                                                                                                                aaaaaGCTAACATTCACCGCGGATGCACGGTCAATATTACTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCgg                <  2:1119807/146‑1 (MQ=255)                                                                                                                                 aataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGt               >  2:1690153/1‑150 (MQ=255)                                                                                                                                    aGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGTGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAga            >  1:936677/1‑149 (MQ=255)                                                                                                                                    aGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGc            >  2:886629/1‑150 (MQ=255)                                                                                                                                    aGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGc            >  2:782481/1‑150 (MQ=255)                                                                                                                                    aGCTAACATTGACCGCGGAAGCACCGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGACTCGCGGGCCAGGGCCTGGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGAGCGCGGTAGc            <  2:1250829/150‑1 (MQ=255)                                                                                                                                     gCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCTGGCTTATGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCAGCTGAGAACGCGGTAgat           >  1:126591/1‑148 (MQ=255)                                                                                                                                     agaaataTTGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCt           >  1:1186015/7‑150 (MQ=255)                                                                                                                                      cTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTg          >  2:1158356/1‑150 (MQ=255)                                                                                                                                      cTAACATTGACCGCGGATGCACGGTCAAGATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCATGCTGGCTTAAGGAATCGCGGGCCAGGGATTTGGAGTTTTTTGCGACGGGAATGCCATTGCGGGTGAGATCGGGGTAGCTg          >  2:1543167/1‑150 (MQ=255)                                                                                                                                       tAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGt         >  1:737776/1‑150 (MQ=255)                                                                                                                                          cATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCtt      >  1:368161/1‑150 (MQ=255)                                                                                                                                           aTTGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACTTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTAAATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGAGGGAGCTGTCtta     >  2:71677/1‑149 (MQ=255)                                                                                                                                            ttGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCa    >  2:532229/1‑150 (MQ=255)                                                                                                                                              gACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCAcc  >  1:663909/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                 ATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACC  >  NC_000913/136384‑136674

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