Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,704,879	T→A	9.5%	Y144F (TAC→TTC) ‡	lepB ←	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,704,879	0	T	A	9.5%	183.5 / 8.9	63	Y144F (TAC→TTC) ‡	lepB	signal peptidase I
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base A (2/4);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

AAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGT  >  NC_000913/2704732‑2705027                                                                                                                                                    |                                                                                                                                                     aaaTCGCTCGGTTCCACGTTTGGGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTaa                                                                                                                                                    >  2:262995/1‑150 (MQ=255)          ggTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcgcg                                                                                                                                           <  1:966552/150‑4 (MQ=255)          ggTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcgcg                                                                                                                                           <  2:1318872/150‑4 (MQ=255)              ccACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCGCGGTTGAATCGTCAGCTCTTTTGAGACCCGATCGTAAGTGACTTTATCGCCCGGTAAACACACCGCGCCCTTGAGGTAATCAAGCTTTGGAt                                                                                                                                       >  1:934293/1‑150 (MQ=255)               cACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATc                                                                                                                                      <  2:1068335/150‑1 (MQ=255)                     tGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCgg                                                                                                                                >  1:803455/1‑150 (MQ=255)                          aGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcgc                                                                                                                                            >  1:1339952/1‑131 (MQ=255)                          aGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcgc                                                                                                                                            <  2:1339952/133‑3 (MQ=255)                          aGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATAtt                                                                                                                           <  1:1252770/150‑1 (MQ=255)                           ggTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATAttt                                                                                                                          >  1:245532/1‑150 (MQ=255)                           ggTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATAttt                                                                                                                          >  2:748188/1‑150 (MQ=255)                           ggTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATATTCCGGATAttt                                                                                                                          >  2:760256/1‑150 (MQ=255)                            gTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTa                                                                                                                         >  2:762835/1‑150 (MQ=255)                            gTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCGTCCGGATATTTa                                                                                                                         >  2:557039/1‑150 (MQ=255)                                  ggCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAcc                                                                                                                   >  2:758626/1‑150 (MQ=255)                                  ggCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGAcc                                                                                                                   <  1:5679/150‑1 (MQ=255)                                    cAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCAc                                                                                                                 <  1:1366904/150‑1 (MQ=255)                                         cgTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGatat                                                                                                            <  1:758626/150‑1 (MQ=255)                                         cgTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGatat                                                                                                            >  1:1087794/1‑150 (MQ=255)                                                 cacGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcg                                                                                                                                             <  2:764474/109‑2 (MQ=255)                                                 cacGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATTACATCAAGcg                                                                                                                                             >  1:764474/1‑108 (MQ=255)                                                  acGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGCCCACGATATCGCCGCGtt                                                                                                   >  1:760593/1‑150 (MQ=255)                                                    gCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAATCCAACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTc                                                                                                 >  2:964457/1‑150 (MQ=255)                                                     ccTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCg                                                                                                >  1:1112968/1‑150 (MQ=255)                                                      cTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCAACCGCGCGCTTGATGTAATCAAGCTTTGGATATTCCGGATATTTAAATACAACGATATCGCCGCGTTTCgg                                                                                               >  2:1384945/1‑150 (MQ=255)                                                         gCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATg                                                                                            <  2:245532/150‑1 (MQ=255)                                                         gCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATg                                                                                            >  2:939101/1‑150 (MQ=255)                                                         gCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATg                                                                                            >  2:943794/1‑150 (MQ=255)                                                                  gCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTc                                                                                   <  2:1087794/150‑1 (MQ=255)                                                                    aTCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGa                                                                                 <  2:795381/150‑1 (MQ=255)                                                                      cccGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATc                                                                               >  2:1410810/1‑150 (MQ=255)                                                                                    cAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGCCCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTa                                                                 >  2:931517/1‑150 (MQ=255)                                                                                    cAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTa                                                                 >  2:529135/1‑150 (MQ=255)                                                                                    cAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTa                                                                 <  1:557039/150‑1 (MQ=255)                                                                                     aGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAg                                                                >  1:1038877/1‑150 (MQ=255)                                                                                     aGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAg                                                                <  2:492682/150‑1 (MQ=255)                                                                                           tttGAGACCGGATCGTAAGTGACTTTATCGCCCTGTAATCCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAACACCACGATATGGCCGCGTTTCGGATGACCGGTTTTGATCAGCGTTTTGTGGTAGATTggt                                                          <  2:508089/150‑2 (MQ=255)                                                                                           tttGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGa                                                          <  1:176326/150‑1 (MQ=255)                                                                                           tttGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGa                                                          >  1:96208/1‑150 (MQ=255)                                                                                           tttGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGa                                                          >  1:652774/1‑150 (MQ=255)                                                                                             tGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATc                                                        <  1:161612/150‑1 (MQ=255)                                                                                                tacgggATCGTAAGTGACTTTATCGCCCGGTAAACCGTCCGCGGGCTTGTTGTAATCAAGGTTTGGATCTTCCGGATATTTAAAGACCACGTTATCGGCTCGTTTCGGATGACTGGTTTCGTTCAGCGTTTTCTGGTAGATAGGATCttt                                                     <  2:934293/146‑1 (MQ=255)                                                                                                   cGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAAt                                                  <  2:491585/150‑1 (MQ=255)                                                                                                    ggATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATg                                                 <  2:3084/150‑1 (MQ=255)                                                                                                       tCGTAAGTGACTTTATCGCCCGGTAAACCCTCCTCGCGCTTGATTTTATCAAGCTTTGTATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCa                                              <  1:1384945/150‑1 (MQ=255)                                                                                                       tCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCa                                              <  2:803455/150‑1 (MQ=255)                                                                                                        cGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCAt                                             <  1:939101/150‑1 (MQ=255)                                                                                                        cGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCAt                                             <  1:943794/150‑1 (MQ=255)                                                                                                           aaGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAg                                          <  2:1355763/150‑1 (MQ=255)                                                                                                                     aTCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtctc                                <  2:1112968/150‑1 (MQ=255)                                                                                                                         cccGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACTATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGTATCTTTAATGCCATAAGCAAACTTCTCTAcc                            <  1:762835/150‑1 (MQ=255)                                                                                                                           cGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtc                                  >  2:1343224/1‑142 (MQ=255)                                                                                                                           cGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtc                                  <  1:1343224/142‑1 (MQ=255)                                                                                                                            ggTAAACCCACCGCGCGCTGGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGa                         <  1:166344/150‑1 (MQ=255)                                                                                                                             gTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGaa                        >  2:734975/1‑150 (MQ=255)                                                                                                                                 aCCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATaaa                    >  2:270400/1‑150 (MQ=255)                                                                                                                                        gcgcgcTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCa             <  1:1410810/150‑1 (MQ=255)                                                                                                                                            gcTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTa         >  1:1019617/1‑150 (MQ=255)                                                                                                                                             cTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTaa        >  1:527231/1‑150 (MQ=255)                                                                                                                                              ttGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAAc       <  1:748188/150‑1 (MQ=255)                                                                                                                                               tGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACa      >  1:1112606/1‑150 (MQ=255)                                                                                                                                                gATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAACCTTCTCTACCAGAATAAAATCACCAATTAACag     >  2:788771/1‑150 (MQ=255)                                                                                                                                                   gTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGt  >  2:959792/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     AAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGT  >  NC_000913/2704732‑2705027

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