Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,738,534	C→T	5.1%	R346R (CGC→CGT)	malS →	alpha‑amylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,738,534	0	C	T	5.1%	273.8 / 3.5	79	R346R (CGC→CGT)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand):  ref base C (34/41);  new base T (2/2);  total (36/43)
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.88e-01

GATAGCGCACATCAGCGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGA  >  NC_000913/3738388‑3738683                                                                                                                                                   |                                                                                                                                                      gATAGCGCACATCAGCGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGTTgg                                                                                                                                                    >  1:1137766/1‑150 (MQ=255) gATAGCGCACATCAGCGAGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGAAGAAGTGAAAAAATCTCTGGGTGAACGTTgg                                                                                                                                                    >  1:1137769/1‑150 (MQ=255)   tAGCGCACATCAGCGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAg                                                                                                                                                  <  2:787364/150‑1 (MQ=255)      cgcACATCAGCGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGa                                                                                                                                               <  1:211016/150‑1 (MQ=255)       gcACATCAGCGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGAc                                                                                                                                              >  2:854196/1‑150 (MQ=255)             cAGCGCGTTATTCGTTTTCTCTTTGATGTGGTGATGAACCACACCGGCTATGCCACGCTGGCGTATATGCAGGAGTATCTGTTTGGCGCGTTATATCTTTCTGGGGACGATGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGaaa                                                                                                                                        <  1:607479/150‑1 (MQ=255)             cAGCGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGaaa                                                                                                                                        >  1:188468/1‑150 (MQ=255)             cAGCGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGaaa                                                                                                                                        <  1:178900/150‑1 (MQ=255)              aGCGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAAc                                                                                                                                       >  1:42897/1‑150 (MQ=255)               gcgcGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAAcc                                                                                                                                      <  2:157044/150‑1 (MQ=255)                        tcgtattcTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCgg                                                                                                                             <  2:1346829/150‑1 (MQ=255)                        tcgtattcTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCgg                                                                                                                             <  2:1346934/150‑1 (MQ=255)                         cgtattcTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCggg                                                                                                                            <  1:373853/150‑1 (MQ=255)                         cgtattcTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCggg                                                                                                                            <  1:1077845/150‑1 (MQ=255)                           tattcTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCa                                                                                                                          >  1:748394/1‑150 (MQ=255)                             ttcTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCaaa                                                                                                                        >  1:145898/1‑150 (MQ=255)                                tctTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCt                                                                                                                     >  2:954441/1‑150 (MQ=255)                                 ctTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTg                                                                                                                    <  1:795699/150‑1 (MQ=255)                                  tttGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTgg                                                                                                                   <  2:74485/150‑1 (MQ=255)                                   ttGATGTCGTGCTGAACCACACCGGCTATGCCACGCGGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTGTCTGGGGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCGGGc                                                                                                                  <  2:465912/150‑1 (MQ=255)                                   ttGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGc                                                                                                                  <  2:135440/150‑1 (MQ=255)                                   ttGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGc                                                                                                                  >  1:671323/1‑150 (MQ=255)                                         tCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCt                                                                                                            >  2:1015718/1‑150 (MQ=255)                                              atgaACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAAc                                                                                                       <  1:1082451/150‑1 (MQ=255)                                                   ccACACCGGCTATGCCACGCTGGCGGATATGCTGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTa                                                                                                  <  2:1180360/150‑1 (MQ=255)                                                   ccACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTa                                                                                                  <  1:1377754/150‑1 (MQ=255)                                                    cacaCCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTAc                                                                                                 >  1:374858/1‑150 (MQ=255)                                                     acacCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACa                                                                                                <  1:964709/150‑1 (MQ=255)                                                         cGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTaa                                                                                            >  2:712036/1‑150 (MQ=255)                                                            cTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAAttt                                                                                         >  2:64095/1‑150 (MQ=255)                                                               tGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAg                                                                                      <  1:701105/150‑1 (MQ=255)                                                               tGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAg                                                                                      <  1:1398758/150‑1 (MQ=255)                                                                 ccACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCg                                                                                    >  2:1001876/1‑150 (MQ=255)                                                                  cACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGa                                                                                   <  2:897684/150‑1 (MQ=255)                                                                    cGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACa                                                                                 >  1:1343901/1‑150 (MQ=255)                                                                         gCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACa                                                                            >  2:184693/1‑150 (MQ=255)                                                                         gCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACa                                                                            >  2:184628/1‑150 (MQ=255)                                                                         gCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACa                                                                            >  2:975803/1‑150 (MQ=255)                                                                         gCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACa                                                                            <  1:364113/150‑1 (MQ=255)                                                                          cGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAg                                                                           >  2:711292/1‑150 (MQ=255)                                                                           ggATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAgg                                                                          >  1:521041/1‑150 (MQ=255)                                                                                  cAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGAt                                                                   >  1:244884/1‑150 (MQ=255)                                                                                   aGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATa                                                                  <  2:1369372/150‑1 (MQ=255)                                                                                          tCAGTTTTGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAAtggtg                                                           <  1:1352997/150‑1 (MQ=255)                                                                                          tCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAAtggtg                                                           <  2:493506/150‑1 (MQ=255)                                                                                           cAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAAtggtgg                                                          <  1:1390909/150‑1 (MQ=255)                                                                                           cAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAAtggtgg                                                          <  2:939200/150‑1 (MQ=255)                                                                                           cAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAAtggtgg                                                          <  2:939208/150‑1 (MQ=255)                                                                                             gTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGTAGCGACTGGAAACCTGCCGCCGGGCAAACCTGTCATAGATTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTgggg                                                        <  1:1016491/150‑1 (MQ=255)                                                                                              tttGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGa                                                       >  1:162738/1‑150 (MQ=255)                                                                                              tttGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGa                                                       >  1:162841/1‑150 (MQ=255)                                                                                              tttGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGa                                                       >  2:1029316/1‑150 (MQ=255)                                                                                                    gcgTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAAc                                                 >  1:224840/1‑150 (MQ=255)                                                                                                     cgTTATATCTGTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACt                                                <  2:1209521/150‑1 (MQ=255)                                                                                                      gTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTg                                               <  2:188468/150‑1 (MQ=255)                                                                                                      gTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTg                                               <  1:309960/150‑1 (MQ=255)                                                                                                        tataTCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGa                                             >  2:1047929/1‑150 (MQ=255)                                                                                                         atatCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGAt                                            <  2:1197785/150‑1 (MQ=255)                                                                                                          tatCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATc                                           >  1:1171758/1‑150 (MQ=255)                                                                                                          tatCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATc                                           >  1:558975/1‑150 (MQ=255)                                                                                                          tatCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCGGGGATACAGGGGGGGGACAAAACTGGATc                                           >  2:966257/1‑150 (MQ=255)                                                                                                           atCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCa                                          <  2:850105/150‑1 (MQ=255)                                                                                                                  tGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGa                                   >  1:1402390/1‑150 (MQ=255)                                                                                                                         gAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGc                            >  1:438721/1‑150 (MQ=255)                                                                                                                         gAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGc                            >  2:496305/1‑150 (MQ=255)                                                                                                                          aaGTGAAAAAATCGCTGGGTGAACGTTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCg                           <  2:1137766/150‑1 (MQ=255)                                                                                                                          aaGTGAAAAAATCGCTGGGTGAACGTTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCg                           <  2:1137769/150‑1 (MQ=255)                                                                                                                               aaaaaaTCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTAc                      >  1:18934/1‑150 (MQ=255)                                                                                                                                 aaaaTCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGa                    <  2:747059/150‑1 (MQ=255)                                                                                                                                 aaaaTCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGa                    <  2:1124156/150‑1 (MQ=255)                                                                                                                                 aaaaTCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGa                    >  1:1220284/1‑150 (MQ=255)                                                                                                                                     tCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAAt                >  2:32059/1‑150 (MQ=255)                                                                                                                                     tCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAAt                >  2:32014/1‑150 (MQ=255)                                                                                                                                        cTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCt             <  2:799740/150‑1 (MQ=255)                                                                                                                                          gtgTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTgg           <  2:672180/148‑1 (MQ=255)                                                                                                                                          gggAGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACCATTACATTAATTTCAGCGACAAAACAGGCTGCGATAAAACGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTgg           <  2:462387/150‑1 (MQ=255)                                                                                                                                           ggTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGa          <  1:1009635/150‑1 (MQ=255)                                                                                                                                             tGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGAtt        <  2:1375709/150‑1 (MQ=255)                                                                                                                                                aCGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTcga     <  1:1286567/150‑1 (MQ=255)                                                                                                                                                   cTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTcgacga  <  2:428002/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GATAGCGCACATCAGCGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGA  >  NC_000913/3738388‑3738683

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