Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,365,448	T→A	6.2%	L406Q (CTG→CAG)	puuB →	gamma‑glutamylputrescine oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,365,448	0	T	A	6.2%	190.5 / 3.3	65	L406Q (CTG→CAG)	puuB	gamma‑glutamylputrescine oxidase
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base A (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

TTACATGCAGGGCTACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGCCACTCCG  >  NC_000913/1365305‑1365594                                                                                                                                                |                                                                                                                                                   ttACATGCAGGGCTACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCgtgt                                                                                                                                              <  2:946367/150‑1 (MQ=255)    cATGCAGGGCTACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGcc                                                                                                                                           >  2:232876/1‑150 (MQ=255)    cATGCAGGGCTACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGcc                                                                                                                                           >  1:258171/1‑150 (MQ=255)      tGCAGGGCTACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGt                                                                                                                                         <  2:282004/150‑1 (MQ=255)       gCAGGGCTACCGCGGCCACGGCGTGACCTGTACTCATCAAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGTTGCGTGTGCCGtt                                                                                                                                        <  2:244336/150‑1 (MQ=255)       gCAGGGCTACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGtt                                                                                                                                        <  1:1323208/150‑1 (MQ=255)        cAGGGCTACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGttt                                                                                                                                       >  2:719195/1‑150 (MQ=255)          gggCTACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTAc                                                                                                                                     >  1:23373/1‑150 (MQ=255)              tACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGtgtg                                                                                                                                             >  2:712747/1‑138 (MQ=255)              tACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGtgtg                                                                                                                                             <  1:712747/138‑1 (MQ=255)                    ggCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGgc                                                                                                                           >  1:912936/1‑150 (MQ=255)                       cACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGcgcg                                                                                                                        <  2:93085/150‑1 (MQ=255)                         cGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCgcggc                                                                                                                      <  2:1258315/150‑1 (MQ=255)                              tGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTAtt                                                                                                                 <  2:60999/150‑1 (MQ=255)                                  cTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAg                                                                                                             <  2:253733/150‑1 (MQ=255)                                    gTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGcc                                                                                                           <  1:816215/150‑1 (MQ=255)                                    gTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGcc                                                                                                           <  1:737286/150‑1 (MQ=255)                                         cATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTgcgc                                                                                                      >  1:698918/1‑150 (MQ=255)                                           tCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGa                                                                                                    >  1:390894/1‑150 (MQ=255)                                           tCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGa                                                                                                    >  1:1113336/1‑150 (MQ=255)                                            cTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCTATTCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGTCGTATTACAGCCTGCGCGAt                                                                                                   <  2:1067971/150‑1 (MQ=255)                                              aGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGTCGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCg                                                                                                 <  1:1316127/150‑1 (MQ=255)                                              aGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCg                                                                                                 <  1:312602/150‑1 (MQ=255)                                              aGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCg                                                                                                 <  1:232876/150‑1 (MQ=255)                                                 cGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCt                                                                                              >  2:183953/1‑150 (MQ=255)                                                      gTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCg                                                                                         >  2:338100/1‑150 (MQ=255)                                                       tttGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGt                                                                                        >  1:62486/1‑150 (MQ=255)                                                       tttGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGt                                                                                        >  2:5068/1‑150 (MQ=255)                                                       tttGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGt                                                                                        >  2:347494/1‑150 (MQ=255)                                                            ttGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAAt                                                                                   >  1:996861/1‑150 (MQ=255)                                                                        tgcGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACcgg                                                                                                                                  <  1:208522/91‑2 (MQ=255)                                                                        tgcGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACcgc                                                                                                                                  >  2:208522/1‑91 (MQ=255)                                                                              gcgACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAag                                                                 <  2:825011/150‑1 (MQ=255)                                                                                   gCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCa                                                            >  2:219765/1‑150 (MQ=255)                                                                                     cGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCaaa                                                          >  1:239196/1‑150 (MQ=255)                                                                                      gAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCaaaa                                                         <  1:183953/150‑1 (MQ=255)                                                                                        aCGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAgg                                                       <  2:23373/150‑1 (MQ=255)                                                                                              cGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGaa                                                 <  2:698918/150‑1 (MQ=255)                                                                                               gATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGaaa                                                <  1:256530/150‑1 (MQ=255)                                                                                                   ccTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATg                                            <  2:964888/150‑1 (MQ=255)                                                                                                       cGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCg                                                                              >  2:229365/1‑112 (MQ=255)                                                                                                       cGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCg                                                                              <  1:229365/112‑1 (MQ=255)                                                                                                       cGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGca                                        >  2:538458/1‑150 (MQ=255)                                                                                                          cAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTaa                                                                         <  2:833902/114‑1 (MQ=255)                                                                                                          cAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTaa                                                                         >  1:833902/1‑114 (MQ=255)                                                                                                          cAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGcaaca                                     >  2:1230362/1‑150 (MQ=255)                                                                                                          cAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGcaaca                                     >  1:1337470/1‑150 (MQ=255)                                                                                                             tCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATg                                  <  2:390894/150‑1 (MQ=255)                                                                                                              cTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGa                                 <  1:986538/150‑1 (MQ=255)                                                                                                                gCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAAt                               >  2:1404661/1‑150 (MQ=255)                                                                                                                  cGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTc                             <  1:888816/150‑1 (MQ=255)                                                                                                                     aTTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCAt                          >  1:943713/1‑150 (MQ=255)                                                                                                                        aCCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAg                       >  2:418421/1‑150 (MQ=255)                                                                                                                          ccGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGcg                     <  2:236417/150‑1 (MQ=255)                                                                                                                                    aacggcaCACGCAGCGTGCGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCt           >  1:1034174/4‑150 (MQ=255)                                                                                                                                    aacggcaCACGCAGCGTGCGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCt           >  1:1039174/4‑150 (MQ=255)                                                                                                                                     cgggcgCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTg          <  1:502643/150‑1 (MQ=255)                                                                                                                                       ggcaCACGCAGCGTGCGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGcc        <  2:1034174/150‑1 (MQ=255)                                                                                                                                       ggcaCACGCAGCGTGCGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGcc        <  2:1039174/150‑1 (MQ=255)                                                                                                                                         cgcACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGCCAc      <  2:912936/150‑1 (MQ=255)                                                                                                                                          gcACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGCCACt     <  2:239196/150‑1 (MQ=255)                                                                                                                                           cACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGCCACTc    >  1:499227/1‑150 (MQ=255)                                                                                                                                             cGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGCCACTCcg  <  1:5068/150‑1 (MQ=255)                                                                                                                                             cGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGCCACTCcg  >  2:430698/1‑150 (MQ=255)                                                                                                                                             cGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGCCACTCcg  >  2:1247020/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                   TTACATGCAGGGCTACAGCGGCCACGGCGTGACCTGTACTCATCTAGCCGGACGTTTGATTGCCGAACTGCTGCGCGGCGACGCCGAACGTTTCGATGCCTTCGCCAATCTGCCGCATTACCCGTTCCCCGGCGGGCGCACGCTGCGTGTGCCGTTTACCGCGATGGGCGCGGCGTATTACAGCCTGCGCGATCGTCTGGGCGTTTAATTTCCGATTAACCGTGAAGAGTCAAAAGGTGTGAAACATGAGCAACAATGAATTCCATCAGCGTCGTCTTTCTGCCACTCCG  >  NC_000913/1365305‑1365594

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