Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,704,880	A→C	9.7%	Y144D (TAC→GAC) ‡	lepB ←	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,704,880	0	A	C	9.7%	171.2 / 11.7	62	Y144D (TAC→GAC) ‡	lepB	signal peptidase I
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base C (2/4);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTCG  >  NC_000913/2704732‑2705029                                                                                                                                                     |                                                                                                                                                      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:652774/1‑150 (MQ=255)                                                                                           tttGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGa                                                            <  1:176326/150‑1 (MQ=255)                                                                                           tttGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGa                                                            >  1:96208/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:943794/150‑1 (MQ=255)                                                                                                        cGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCAt                                               <  1:939101/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                                    <  1:1343224/142‑1 (MQ=255)                                                                                                                           cGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTtc                                    >  2:1343224/1‑142 (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)                                                                                                                                                     aaTCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGTCCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTCg  >  2:494062/1‑150 (MQ=255)                                                                                                                                                     aaTCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTCg  >  1:407807/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AAATCGCTCGGTTCCACGTTTGAGTAGGTGACCGGCAGCGCGTTTTCACACGCCTGGCCGGAACTGCATCCCGGTTGAATCGTCAGCTCTTTTGAGACCGGATCGTAAGTGACTTTATCGCCCGGTAAACCCACCGCGCGCTTGATGTAATCAAGCTTTGGATCTTCCGGATATTTAAAGACCACGATATCGCCGCGTTTCGGATGACCGGTTTCGATCAGCGTTTTCTGGTAGATAGGATCTTTAATGCCATAAGCAAACTTCTCTACCAGAATAAAATCACCAATTAACAGAGTCG  >  NC_000913/2704732‑2705029

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