Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	961,154	A→T	7.1%	intergenic (+126/‑47)	ycaL → / → cmk	periplasmic protease YcaL/cytidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	961,154	0	A	T	7.1%	220.6 / 2.5	73	intergenic (+126/‑47)	ycaL/cmk	periplasmic protease YcaL/cytidylate kinase
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base T (2/4);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGTATCGCCTCTGGTAAGTAAATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTG  >  NC_000913/961005‑961302                                                                                                                                                      |                                                                                                                                                     gATCGTATCGCCTCTGGTAAGTAAATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGtata                                                                                                                                                      <  1:1205505/150‑1 (MQ=255)   tCGTATCGCCTCTGGTAAGTAAATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAAt                                                                                                                                                    >  2:572209/1‑150 (MQ=255)    cGTATCGCCTCTGGTAAGTAAATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATg                                                                                                                                                   >  1:1542923/1‑150 (MQ=255)             gctGGTAAGTAAATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGCATTATacgcgccgac                                                                                                                                          <  1:1351219/149‑5 (MQ=255)                ggTAAGTAAATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATg                                                                                                                                       <  2:544380/150‑1 (MQ=255)                ggTAAGTAAATCATTGTCATCTTTCGGGCGGGGCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACACACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGGCGGCGCGCATTATacgcgccgacgcc                                                                                                                                       <  2:1342573/150‑8 (MQ=255)                   aagtaaATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTa                                                                                                                                    >  2:1617191/1‑150 (MQ=255)                    agtaaATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTaa                                                                                                                                   <  1:256272/150‑1 (MQ=255)                      taaATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACg                                                                                                                                 <  2:764489/150‑1 (MQ=255)                      taaATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACg                                                                                                                                 <  2:1525551/150‑1 (MQ=255)                        aaTCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGt                                                                                                                               >  2:1171548/1‑150 (MQ=255)                          tCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTAc                                                                                                                             <  1:1534857/150‑1 (MQ=255)                               gTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTa                                                                                                                                    >  1:822927/1‑138 (MQ=255)                               gTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTa                                                                                                                                    <  2:822927/138‑1 (MQ=255)                                tCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGt                                                                                                                       <  1:631593/150‑1 (MQ=255)                                 cATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGtt                                                                                                                      <  2:101418/150‑1 (MQ=255)                                 cATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGAAAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGtt                                                                                                                      >  2:1597074/1‑150 (MQ=255)                                      ttCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAg                                                                                                                 <  2:644537/150‑1 (MQ=255)                                      ttCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAg                                                                                                                 <  2:1542923/150‑1 (MQ=255)                                       tCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAgg                                                                                                                <  1:1321833/150‑1 (MQ=255)                                       tCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAgg                                                                                                                >  1:442248/1‑150 (MQ=255)                                          ggCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGgaga                                                                                                             >  1:978618/1‑150 (MQ=255)                                               gTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGGCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAg                                                                                                        <  2:610134/150‑1 (MQ=255)                                                       ccAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGc                                                                                                <  1:779113/150‑1 (MQ=255)                                                         aGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCaa                                                                                              >  1:1228030/1‑150 (MQ=255)                                                          gCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAAt                                                                                             >  1:411684/1‑150 (MQ=255)                                                              ggtttaATTTCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGcc                                                                                         <  1:9339/146‑1 (MQ=255)                                                                       gcgcAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGttatt                                                                                >  2:943312/1‑150 (MQ=255)                                                                        cgcAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGttatta                                                                               >  2:1620901/1‑150 (MQ=255)                                                                         gcAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTAc                                                                              >  2:251872/1‑150 (MQ=255)                                                                          cAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTAcc                                                                             <  2:473494/150‑1 (MQ=255)                                                                          cAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCAGGTTATTAcc                                                                             >  1:961675/1‑150 (MQ=255)                                                                            aTAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCAt                                                                           <  2:1490826/150‑1 (MQ=255)                                                                             tAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCAtt                                                                          >  1:1090954/1‑150 (MQ=255)                                                                             tAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCAtt                                                                          <  1:242381/150‑1 (MQ=255)                                                                                 tCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATg                                                                      >  2:1532442/1‑150 (MQ=255)                                                                                           aaTAGAGACAATTATGGTTTTTGGACACAACTTTAACGATTTTTTGCGTCGGCGCGTATTATGCGCGAGGTTATGTTAACGGTACGCCTTTTTTAATGATATAAAGATGACGGGAATTGCCCCGGTTATTACCATTGATGGCCGAAGCgg                                                            <  2:384865/150‑1 (MQ=255)                                                                                               cAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCa                                                        <  1:1120268/150‑1 (MQ=255)                                                                                                aGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAg                                                       >  2:1396839/1‑150 (MQ=255)                                                                                                 gACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAgg                                                      >  2:1331633/1‑150 (MQ=255)                                                                                                 gACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAgg                                                      >  1:1129148/1‑150 (MQ=255)                                                                                                  aCAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAggg                                                     <  2:257762/150‑1 (MQ=255)                                                                                                   cAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGa                                                    <  1:1617191/150‑1 (MQ=255)                                                                                                      aaCTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGAGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGAAATTGCCCCGGTTATTACAATTGATGGCCCAAGCGGTGCAGGGAAAg                                                 >  2:191303/1‑150 (MQ=255)                                                                                                       aCTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAgg                                                >  1:1592647/1‑150 (MQ=255)                                                                                                          ggTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGc                                               <  1:479463/148‑1 (MQ=255)                                                                                                          ggTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGc                                               >  2:479463/1‑148 (MQ=255)                                                                                                          ggTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCAc                                             >  2:323022/1‑150 (MQ=255)                                                                                                           gTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCAcc                                            <  2:725913/150‑1 (MQ=255)                                                                                                             ttttGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCtt                                          >  2:540280/1‑150 (MQ=255)                                                                                                              tttGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTtg                                         <  2:978618/150‑1 (MQ=255)                                                                                                              tttGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTtg                                         <  1:23920/150‑1 (MQ=255)                                                                                                                tGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTtgtg                                       >  1:768361/1‑150 (MQ=255)                                                                                                                tGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTtgtg                                       <  2:719831/150‑1 (MQ=255)                                                                                                                  cacacaACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTa                                     <  2:442248/150‑1 (MQ=255)                                                                                                                      caACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGc                                 >  1:1305256/1‑150 (MQ=255)                                                                                                                        aCGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTa                               >  1:507551/1‑150 (MQ=255)                                                                                                                         cGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTAt                              >  2:1080167/1‑150 (MQ=255)                                                                                                                          gTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGGAAGGCTATg                             >  2:611997/1‑150 (MQ=255)                                                                                                                            tAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGc                           >  1:1163694/1‑150 (MQ=255)                                                                                                                               cGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGa                        >  1:518049/1‑150 (MQ=255)                                                                                                                               cGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGa                        <  1:251872/150‑1 (MQ=255)                                                                                                                                 aTTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAg                      <  2:1228030/150‑1 (MQ=255)                                                                                                                                  tttGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGc                     >  1:1595425/1‑150 (MQ=255)                                                                                                                                       ggCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGc                <  2:692545/150‑1 (MQ=255)                                                                                                                                       ggCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGc                <  1:129352/150‑1 (MQ=255)                                                                                                                                       ggCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGc                >  2:962657/1‑150 (MQ=255)                                                                                                                                         taaccgcgcgcATTATACGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCaa              >  1:1587831/6‑150 (MQ=255)                                                                                                                                           tCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATg            <  2:1595425/150‑1 (MQ=255)                                                                                                                                           accgcgcgcATTATACGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATg                             <  2:103718/130‑1 (MQ=255)                                                                                                                                           accgcgcgcATTATACGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATg                             >  1:103718/4‑133 (MQ=255)                                                                                                                                              gcgcGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCa         <  1:1171548/150‑1 (MQ=255)                                                                                                                                               cgcgTATAATGCGCGCGGTTATGTTAACGGTACGCATGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTCTAAGGCTATGGCGGAAGCGTAGCAATGGCAt        >  2:613640/1‑150 (MQ=255)                                                                                                                                                  gTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATctg     >  1:951059/1‑150 (MQ=255)                                                                                                                                                     taATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATctgctg  >  1:702065/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GATCGTATCGCCTCTGGTAAGTAAATCATTGTCATCTTTCGGGCTGGTCTTCTGCCAGCCCGCTATAATTGCGCAATAAATCCCCATCTGAATACAGACAAAACTGGTTTTTGCACACAACGTTAACGATTTGTGGCGTCGGCGCGTATAATGCGCGCGGTTATGTTAACGGTACGCCTGTTTTAAGGAGATAAAGATGACGGCAATTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTG  >  NC_000913/961005‑961302

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: ‑