Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	691,321	A→T	L155Q (CTG→CAG)	ybeX ←	putative ion transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	691,321	0	A	T	98.5%	241.3 / ‑3.3	65	L155Q (CTG→CAG)	ybeX	putative ion transport
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (28/36);  total (28/37)
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.25e-01

CGATATCATCTTCTTCGTCATACTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTA  >  NC_000913/691178‑691467                                                                                                                                                |                                                                                                                                                   cGATATCATCTTCTTCGTCATACTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATcc                                                                                                                                              <  1:991156/150‑1 (MQ=255) cGATATCATCTTCTTCGTCATACTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCg                                                                                                                                             >  2:69160/1‑151 (MQ=255)       catcTTCTTCGTCATACTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTa                                                                                                                                        <  2:665049/150‑1 (MQ=255)            tcttcGTCATACTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCt                                                                                                                                  >  1:706781/1‑151 (MQ=255)               tcGTCATACTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCt                                                                                                                                  <  1:67702/148‑1 (MQ=255)                cGTCATACTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACt                                                                                                                              >  2:255627/1‑151 (MQ=255)                    aTACTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCa                                                                                                                          <  2:1068208/151‑1 (MQ=255)                      aCTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAg                                                                                                                         <  1:799901/150‑1 (MQ=255)                         cGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAAc                                                                                                                     >  2:265470/1‑151 (MQ=255)                          gTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACg                                                                                                                    >  2:656866/1‑151 (MQ=255)                             ttCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACg                                                                                                                    >  2:727169/1‑148 (MQ=255)                              tCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACa                                                                                                                 <  2:686126/150‑1 (MQ=255)                              tCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACa                                                                                                                 <  1:1163642/150‑1 (MQ=255)                                 aTCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGAAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAAcc                                                                                                              >  1:954845/1‑150 (MQ=255)                                     cACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCg                                                                                                             <  1:796090/147‑1 (MQ=255)                                       ccAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGAc                                                                                                       >  1:1062069/1‑151 (MQ=255)                                       ccAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGAc                                                                                                       >  2:267089/1‑151 (MQ=255)                                       ccAACAACCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTAAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTg                                                                                                         <  2:500864/149‑1 (MQ=255)                                        caacaaTCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACg                                                                                                      <  2:1071421/151‑1 (MQ=255)                                          acaaACAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTa                                                                                                    <  2:92743/151‑1 (MQ=255)                                             aTCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAac                                                                                                  >  2:855378/1‑150 (MQ=255)                                                   tCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTg                                                                                            <  1:518532/150‑1 (MQ=255)                                                     cAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGATAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTcc                                                                                         >  1:1181327/1‑151 (MQ=255)                                                        gATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATg                                                                                      >  1:154669/1‑151 (MQ=255)                                                              ttCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGa                                                                                  >  2:974911/1‑149 (MQ=255)                                                               tCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAg                                                                                <  1:981142/150‑1 (MQ=255)                                                                cAATGGTCCCCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGaa                                                                                 <  2:122216/148‑1 (MQ=255)                                                                   tGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTc                                                                           >  1:558839/1‑151 (MQ=255)                                                                       cACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGc                                                                        <  2:347092/150‑1 (MQ=255)                                                                         ccAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCAt                                                                      <  1:960226/150‑1 (MQ=255)                                                                                    aCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAGCTCTTTCTGCATCCGTTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATaa                                                          <  2:511434/151‑1 (MQ=255)                                                                                    aCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATa                                                           >  2:360413/1‑150 (MQ=255)                                                                                    aCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATa                                                           >  1:1173922/1‑150 (MQ=255)                                                                                       ccACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACg                                                       <  1:595457/151‑1 (MQ=255)                                                                                        cACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACgg                                                      <  1:686920/151‑1 (MQ=255)                                                                                         aCCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgc                                                     >  1:603720/1‑151 (MQ=255)                                                                                          ccGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAATCTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACGCTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgc                                                     <  2:189742/150‑1 (MQ=255)                                                                                            gAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgcagc                                                  >  2:416618/1‑151 (MQ=255)                                                                                                tCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAg                                               <  1:217381/150‑1 (MQ=255)                                                                                                  gTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCt                                            >  1:1184266/1‑151 (MQ=255)                                                                                                   tCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCt                                            <  1:134782/150‑1 (MQ=255)                                                                                                    cAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCttt                                          >  2:701246/1‑151 (MQ=255)                                                                                                       tAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCg                                        >  2:482675/1‑150 (MQ=255)                                                                                                        aaCGATCGCCATGGGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGc                                       <  2:118632/150‑1 (MQ=255)                                                                                                         aCGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCa                                     <  2:444985/151‑1 (MQ=255)                                                                                                          cGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCa                                     <  2:730519/150‑1 (MQ=255)                                                                                                           gATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATc                                   >  1:868337/1‑151 (MQ=255)                                                                                                            aTCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATc                                   >  1:854833/1‑150 (MQ=255)                                                                                                             tCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCa                                  <  2:95741/150‑1 (MQ=255)                                                                                                              cGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAg                                 <  2:1107100/150‑1 (MQ=255)                                                                                                               gCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGa                                >  1:1128348/1‑150 (MQ=255)                                                                                                               gCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGa                                >  2:278272/1‑150 (MQ=255)                                                                                                               gCCATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGa                                >  1:964615/1‑150 (MQ=255)                                                                                                                 cATGTGGTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATc                             <  2:553618/151‑1 (MQ=255)                                                                                                                      ggTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACGCTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATcccgt                         <  2:779248/150‑3 (MQ=255)                                                                                                                      ggTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCtt                         <  1:531/150‑1 (MQ=255)                                                                                                                      ggTAACGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTc                        <  1:512945/151‑1 (MQ=255)                                                                                                                         aaCGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCaa                      <  1:992952/150‑1 (MQ=255)                                                                                                                           cGCTGAGAGCGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtgt                   <  1:570978/151‑1 (MQ=255)                                                                                                                                    cGAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg           <  2:737402/150‑1 (MQ=255)                                                                                                                                     gAAACTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTc         >  1:1048017/1‑151 (MQ=255)                                                                                                                                      aaaCTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCt        >  1:1158451/1‑151 (MQ=255)                                                                                                                                       aaCTCTTTCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCtt       <  1:1045535/151‑1 (MQ=255)                                                                                                                                            tttCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCg     <  2:1138819/148‑1 (MQ=255)                                                                                                                                             ttCTGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTa  <  2:1089937/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CGATATCATCTTCTTCGTCATACTCGTCTTCAATCTCACCAACAATCAGTTCCAGGATGTCTTCAATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTA  >  NC_000913/691178‑691467

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