Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,910,295	A→T	12.8%	D7V (GAC→GTC)	yebQ →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,910,295	0	A	T	12.8%	189.1 / 12.3	83	D7V (GAC→GTC)	yebQ	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (32/39);  new base T (4/6);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGG  >  NC_000913/1910146‑1910441                                                                                                                                                      |                                                                                                                                                   attatAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGa                                                                                                                                                    >  1:136279/1‑150 (MQ=255)  ttatAATTTATCACAATGGTGGTAGGAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGAc                                                                                                                                                   <  1:346778/150‑1 (MQ=255)    ataATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGTCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACgg                                                                                                                                                 >  1:603244/1‑150 (MQ=255)    ataATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACgg                                                                                                                                                 >  2:1142090/1‑150 (MQ=255)    ataATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACgg                                                                                                                                                 <  1:200921/150‑1 (MQ=255)    ataATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACgg                                                                                                                                                 >  1:1597239/1‑150 (MQ=255)      aaTTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGcc                                                                                                                                               >  2:439500/1‑150 (MQ=255)      aaTTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATAATACCCCTGTTGAGTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGcc                                                                                                                                               >  2:82967/1‑150 (MQ=255)          tCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaa                                                                                                                                           <  2:748760/150‑3 (MQ=255)          tCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaa                                                                                                                                           <  1:1685081/150‑3 (MQ=255)          tCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGGGAATCTCAGTGCTAGAATCATACCCCTGTTGATTATTCTCCAAAGATATAAAATTCCTATGCCAAAAGATCAGGCCGGCGGCCTgaa                                                                                                                                           <  1:30702/150‑3 (MQ=255)            gCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCAt                                                                                                                                         <  1:1231159/150‑1 (MQ=255)             cACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaactt                                                                                                                                        <  1:1317879/150‑6 (MQ=255)              aCCAGGGTGGTCGCAATCCATCTTTTG‑CTGTTAGTTACAATTCTGCGACATCCACCGTGAATTTCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaactttt                                                                                                                                      <  1:1129271/150‑8 (MQ=255)                cAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGcc                                                                                                                                     >  1:593981/1‑150 (MQ=255)                 aGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGccc                                                                                                                                    >  1:55758/1‑150 (MQ=255)                   ggtggtCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCa                                                                                                                                  <  1:227227/150‑1 (MQ=255)                         cGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATa                                                                                                                            <  1:571420/150‑1 (MQ=255)                          gCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATAc                                                                                                                           <  2:1205315/150‑1 (MQ=255)                             aTCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCGCTGTTGATTATTCACCAAAGATATAAAATTCCTATGACAAAAGTTCAGGCCGACGGCCTGCCCTTGCCCCAGCGATACGGt                                                                                                                        >  1:887411/1‑150 (MQ=255)                              tCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTg                                                                                                                       <  2:326299/150‑1 (MQ=255)                                 aTCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGa                                                                                                                    >  1:1070089/1‑150 (MQ=255)                                     tttGCCGGTTAGTTACAATTCTGCGACATCCACCCTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGATGGCCTGCCATTGCCCCAGCGATACGGTGCGATAtt                                                                                                                <  1:72854/150‑1 (MQ=255)                                          cGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCa                                                                                                           >  2:1249042/1‑150 (MQ=255)                                             ttagttaCAATTCTGCGACATCCACCGTGAATATCTTTGCTAGAATGATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAATAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATtg                                                                                                        <  2:1105995/150‑1 (MQ=255)                                                gttaCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGa                                                                                                     >  1:1284945/1‑150 (MQ=255)                                                 ttaCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGAt                                                                                                    >  2:428792/1‑150 (MQ=255)                                                 ttaCAATTCTGCGACATCCACAGTTAATATCATTGCTAGAATCATAACCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGAt                                                                                                    <  2:908450/150‑1 (MQ=255)                                                   aCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTg                                                                                                  <  1:483391/150‑1 (MQ=255)                                                   aCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTg                                                                                                  <  2:647368/150‑1 (MQ=255)                                                    cAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTgg                                                                                                 <  2:1070089/150‑1 (MQ=255)                                                      aTTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACTAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTa                                                                                               >  2:683385/1‑150 (MQ=255)                                                         cttcGACATCAACCGTGAATATCAGAGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAAATCCTATGCCAAAAGTTCAGTCCGACGGCCTGCCATTGCCCCAGCGATACGTTGCGATATTAACCATTGTGATTGGTAttt                                                                                            <  1:749468/147‑1 (MQ=255)                                                           gCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCg                                                                                          >  2:1334446/1‑150 (MQ=255)                                                             gACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGAt                                                                                        >  2:499520/1‑150 (MQ=255)                                                               cATCCCCCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATgg                                                                                      <  2:1046512/150‑1 (MQ=255)                                                               cATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGCCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATgg                                                                                      >  2:55135/1‑150 (MQ=255)                                                                  ccACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGCCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCg                                                                                   <  1:1489461/150‑1 (MQ=255)                                                                    aCCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTc                                                                                 >  1:1714383/1‑150 (MQ=255)                                                                     ccGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGAGATATTAACCATTGTGATTGGTATTTCGATgtccgtcc                                                                                >  2:1162853/1‑150 (MQ=255)                                                                          aaTATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGAc                                                                           <  1:545216/150‑1 (MQ=255)                                                                              tCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgcg                                                                       <  2:646163/150‑1 (MQ=255)                                                                              tCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgcg                                                                       <  1:387958/150‑1 (MQ=255)                                                                                aGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCa                                                                     >  1:1254505/1‑150 (MQ=255)                                                                                 gTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCAAGGACGgccatc                                                                    <  2:724878/150‑5 (MQ=255)                                                                                   gCTAGAATCATACCCCTGTTTATTATTCACCAAAGATATAAAATTCATATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCCCTCCGGTGCGATAATAACCATTGAGCGTGCTATTTACATCGCCGTCCTTGACGGCGCAATc                                                                  >  2:1035558/1‑150 (MQ=255)                                                                                   gCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATc                                                                  <  1:1263766/150‑1 (MQ=255)                                                                                    cTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCg                                                                 <  2:136279/150‑1 (MQ=255)                                                                                    cTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAACTGCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCTTCCTTGACGGCGCAATCg                                                                 <  1:823462/150‑1 (MQ=255)                                                                                      aGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGcc                                                               <  1:1578530/150‑1 (MQ=255)                                                                                           cATACCCCTGTTGATTAATCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGt                                                          <  1:737324/150‑1 (MQ=255)                                                                                               ccccTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGcc                                                      >  2:1389873/1‑150 (MQ=255)                                                                                                    gTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGcc                                                 >  2:955891/1‑150 (MQ=255)                                                                                                    gTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGCTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGcc                                                 <  1:1175119/150‑1 (MQ=255)                                                                                                    gTTGATTATTCACCAAAGATATAAAATTCCTAAGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGcc                                                 <  2:522765/150‑1 (MQ=255)                                                                                                     ttGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCTATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCca                                                <  1:937996/150‑1 (MQ=255)                                                                                                     ttGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCca                                                <  1:255031/150‑1 (MQ=255)                                                                                                      tGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGCCGGCGCAATCGCCAACGTCGCCCTGCcaa                                               >  2:931914/1‑150 (MQ=255)                                                                                                       gATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCcaac                                              <  2:198040/150‑1 (MQ=255)                                                                                                        attattCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCcaaca                                             <  1:1142090/150‑1 (MQ=255)                                                                                                           attCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATc                                          <  1:499520/150‑1 (MQ=255)                                                                                                             tCACTAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGc                                        <  1:683385/150‑1 (MQ=255)                                                                                                             tCACCAAAGATATAAACTTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCGCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGc                                        <  1:265923/150‑1 (MQ=255)                                                                                                               aCCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGCCGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGccc                                      >  2:35064/1‑149 (MQ=255)                                                                                                                ccAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCAc                                     <  1:52137/150‑1 (MQ=255)                                                                                                                           aaaTTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATg                          <  2:683503/150‑1 (MQ=255)                                                                                                                            aaTTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCCTTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCAAACGTCGCCATGACAACAATAGCCACGGACCGTCAGGc                         >  2:869212/1‑150 (MQ=255)                                                                                                                            aaTTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCTTTGTGATTGGTATTTCGATGGCCGTCCTTGTCGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGc                         <  1:82967/150‑1 (MQ=255)                                                                                                                                ccTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGTGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACg                     >  1:1239935/1‑150 (MQ=255)                                                                                                                                ccTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACg                     >  1:958922/1‑150 (MQ=255)                                                                                                                                     gCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCATGCCAACAATCGCCAGGGACCTTCATGCCACgcccgc                >  2:148706/1‑150 (MQ=255)                                                                                                                                      ccAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgccagcc               >  1:1171845/1‑150 (MQ=255)                                                                                                                                       cAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACgccagcca              <  1:490112/150‑1 (MQ=255)                                                                                                                                         aaaGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGt            >  1:478513/1‑150 (MQ=255)                                                                                                                                          aatggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCGTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGtt           >  1:412458/6‑150 (MQ=255)                                                                                                                                           atggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTc          >  1:329526/5‑150 (MQ=255)                                                                                                                                           atggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTc          >  1:999254/5‑150 (MQ=255)                                                                                                                                           aGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTc          <  2:593981/150‑1 (MQ=255)                                                                                                                                           aGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTc          >  1:1306637/1‑150 (MQ=255)                                                                                                                                             ggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCTGCCAGTTCCa        >  2:797091/3‑150 (MQ=255)                                                                                                                                             ggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCTGCCAGTTCCa        <  1:797091/148‑1 (MQ=255)                                                                                                                                                aGGCCGACGGCCTGCCATTGCCCCAGCGATACGGAGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCt     >  2:317293/1‑150 (MQ=255)                                                                                                                                                 ggCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTg    >  1:139706/1‑150 (MQ=255)                                                                                                                                                   ccGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTggg  <  2:329526/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   ATTATAATTTCGCACCAGGGTGGTCGCAATCCATCTTTTGCCGGTTAGTTACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCCAACAATCGCCACGGACCTTCATGCCACGCCAGCCAGTTCCATCTGGG  >  NC_000913/1910146‑1910441

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