Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,991,335	T→A	7.3%	I61I (ATT→ATA)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,335	0	T	A	7.3%	188.8 / 2.6	73	I61I (ATT→ATA)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (35/29);  new base A (3/3);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

ATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGG  >  NC_000913/3991195‑3991476                                                                                                                                             |                                                                                                                                              aTAAATCAATTGCGTGTGGATCGCGCTCTTGCTGCTATCATGACTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTg                                                                                                                                             >  2:374594/1‑143 (MQ=255)  tAAATTAATTGCGTGTGTATCGCGCGCGTGCGGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCTCTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGc                                                                                                                                            <  1:733793/143‑1 (MQ=255)    aaTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCAGGGTTACCTTGATGGTAACGTTCCCCAAGGCATTTGCCt                                                                                                                                          >  2:775870/1‑143 (MQ=255)    aaTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCAGCCATCGGTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGcgt                                                                                                                                          >  1:613410/1‑141 (MQ=255)     aTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtt                                                                                                                                         >  1:206070/1‑143 (MQ=255)       cgaTTGCGTGTGGATCGCGCGCTAGCTGCGATGGGGCCTGCGTTCCAACAGGTCTACAGTCTACTGCCGACAGTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                       <  2:422017/141‑2 (MQ=255)       cAATTGCGTGTGGATCGCGCGCTTGCTGCTACGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                                                                                                       <  1:538615/143‑2 (MQ=255)       cAATTGCGTGTGGATCGCGCGCGTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGATACCTTGATGGTCACGCTCCCAAAGGCATTTGCCTTTa                                                                                                                                       >  2:718641/1‑143 (MQ=255)         aTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttggg                                                                                                                                     <  2:206070/143‑4 (MQ=255)                gtgGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                              >  1:550611/1‑143 (MQ=255)                   gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgatga                                                                                                                           >  2:576314/1‑143 (MQ=255)                   gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTCCCGACATTGTTGCCCTATCACCATCCGCTAATGCCCGGTTACCTTGATGGTACCGTTCCCAAAGGCACTTGCCTTCACTCGcattctaa                                                                                                                           >  2:1072503/1‑136 (MQ=255)                   gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACACTCTACTGCCGACATCGTTGCCCTATCACCATCCGCCAACGCCCGGTTACCTCGATGGTAACGCTCCCAAAGGCAGCCGCCTTCACACGCCtgatgc                                                                                                                           >  2:56024/1‑142 (MQ=255)                     tCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaag                                                                                                                         <  2:528361/143‑2 (MQ=255)                     tCGCGCGCTGGCTGCTATGGGGCCTTCATCTCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaaa                                                                                                                         <  2:967141/143‑1 (MQ=255)                     tCGCCCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCACCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaaa                                                                                                                         <  1:149534/143‑1 (MQ=255)                       gcgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACt                                                                                                                       <  1:1023787/143‑1 (MQ=255)                        cgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTCACCTTGATGGTAACGTTCCCAAAGGCATATGCCTTTACACGCCTGATGAACCTc                                                                                                                      >  2:8831/1‑143 (MQ=255)                             ttGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGc                                                                                                                 >  1:575043/1‑143 (MQ=255)                               gctgctATGGGGCCTGCATTCCAACAGGTCTCCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCa                                                                                                               <  1:1044681/143‑1 (MQ=255)                               gctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCCTCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAACGCCTTTGCCTTTACACGCCGGATGATACTCACCGccc                                                                                                               >  1:181149/1‑142 (MQ=255)                                  gctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                                                                                                       >  1:800118/1‑116 (MQ=255)                                  gctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                                                                                                       <  2:800118/116‑1 (MQ=255)                                    tATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTAcc                                                                                                          <  2:347789/143‑1 (MQ=255)                                     attggtgcTGTATTCCAACCGGACAACAGTCTACTGCCGACACTGGTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCt                                                                                                         <  1:520409/136‑1 (MQ=255)                                     aTGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGGTACCTTGATGGTAACGTTCCCAAACGCATTTGCCTTTACACGCCGGATGAAACACAACGCCAcgagct                                                                                                         >  2:791415/1‑138 (MQ=255)                                      tGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACCCGCCTGATGAAACCCAACGCCGCTACCTg                                                                                                        >  2:141561/1‑143 (MQ=255)                                        gggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGATAGTGTTTCACTATCACCATCCGCTTAGCCCGGGGTACCTTGATGGTAACGTTCCCAAAGGCATTCGCCTTTACACGCCTGATGACACTCATCGTCACTACTTGaa                                                                                                      >  1:252320/1‑143 (MQ=255)                                            cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCGCCAGCAGATAATGCCCGGTTCCCGGGATGGTCACGTTCGCACAGGCATGTGCCGTTACACGCCTGCTGAAACACAGCGACACTCCCTgcgggcg                                                                                                  >  2:178977/1‑137 (MQ=255)                                            cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAg                                                                                                  >  1:45341/1‑143 (MQ=255)                                            cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCGTTACACGCCTGATGAAACTCAACACCACTACCTCACCGac                                                                                                  >  1:368113/1‑142 (MQ=255)                                             tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                                                                                                 >  1:860196/1‑143 (MQ=255)                                             tGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                                                                                                 <  2:603454/143‑1 (MQ=255)                                             cgCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGc                                                                                                 <  1:550949/142‑1 (MQ=255)                                               cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTCGCCTTTACACGCCTGATCAAACTCAACGCCACTACCTGAACCAGCtt                                                                                               >  2:960222/1‑143 (MQ=255)                                                  ttcAACACGTCTCCAGTCTACGGCCGACATTGTTGCGCTATCACCCTCCGCTAATTCCGGGTTATCTTGATGGTAACGTTCCCAAATCGATTTGCCTTTACACGTCTGACGTAGCTCAACGCCAATATCTGAACGAGCTTGaa                                                                                            <  2:181553/141‑1 (MQ=255)                                                  tCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                            <  2:550611/143‑1 (MQ=255)                                                  gccATCAGGTCTACGGACGACTGCCGCCGTTGTGGCAGTATCGCGATCCGCGAGTGCCGGGTTACCTGGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCAGTACCTGAACGAGCTTGaa                                                                                            <  2:1037654/142‑1 (MQ=255)                                                    cAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACt                                                                                          >  2:736696/1‑143 (MQ=255)                                                     aaCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAACTTGAACTg                                                                                         <  1:603853/143‑1 (MQ=255)                                                       cAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGGTCCCGATGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTa                                                                                       <  1:449303/143‑1 (MQ=255)                                                        aGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTACTGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAACCTCAACGCACCTACCTGACCGAGCTTGAACTGTAt                                                                                      >  2:829974/1‑143 (MQ=255)                                                            cgaCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCGTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTg                                                                                  <  1:352841/141‑1 (MQ=255)                                                            cTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATCAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTg                                                                                  >  1:195419/1‑143 (MQ=255)                                                                  tCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGt                                                                            <  1:576314/143‑1 (MQ=255)                                                                   gtaatGCCGACATGGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                           <  1:43454/139‑1 (MQ=255)                                                                   cTACTGCCGACATTGTGGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                           <  1:142199/143‑1 (MQ=255)                                                                    tgcTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCa                                                                          <  2:1068360/141‑1 (MQ=255)                                                                       tGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGTCTTTACATGTTTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTa                                                                       <  1:527193/143‑1 (MQ=255)                                                                       tGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGCTCCCTTGATGGTAACGTTCCCAAAGGCATTTAACATTAAACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGACTGTCAGTa                                                                       >  1:29442/1‑143 (MQ=255)                                                                           gACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGACTTTACACGCCTGATGAAACTTAACGGCACTACCTTAACGAACTTGAAATGTATAATCGAATGTCAGCACAga                                                                   >  2:102235/1‑142 (MQ=255)                                                                                 gttgCACTATCACCATCCGCTAATGCCGAGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTAACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCCTGAACTGCATCGTGGAATGCCAGTACAGGATccgc                                                             >  2:432597/1‑143 (MQ=255)                                                                                     cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAGTTGCCGTTACACGCCTGATGAAACTCAACGCCACTACCTGACCGCGCTTGAACTGTATCGTGGAACGTCAGTACAGGAACCGCCGaa                                                         >  2:238085/1‑143 (MQ=255)                                                                                         aTCACCATCCGCTAATGCCGGGTTACCTTGAGGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGCGCTTGGGCCGTAGCGTCGCATGGCAGCACAGCACCCGCCGAAAGGt                                                     >  2:679829/1‑143 (MQ=255)                                                                                              cATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                >  2:29123/1‑143 (MQ=255)                                                                                               aTCCGCTAATGCCGGGTTACCGTGATGGTAACGTTCCCAACGGCATATGCCTTTACACGACTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGGGGAATGTCAGTACAGGATCCGCCGAAAGGTGGGCtt                                               >  2:57975/1‑143 (MQ=255)                                                                                                    cTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCGTCCact                                          >  1:468517/1‑141 (MQ=255)                                                                                                     tAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTCCACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTACCTGCTCATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAAtt                                         >  1:667703/1‑143 (MQ=255)                                                                                                      aaTGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACACGATCCGCCGAAAGGTGACCTTCCaccta                                        >  1:415156/1‑139 (MQ=255)                                                                                                       aTGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTAc                                       <  1:736696/143‑1 (MQ=255)                                                                                                         gCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                     >  1:322536/1‑143 (MQ=255)                                                                                                                tgccTTGATGGTAACGTTCCCAAAGGCATGTGCCTTTACACGCCTGAGGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGGACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGtata                              <  1:1072503/141‑1 (MQ=255)                                                                                                                tACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGtata                              <  1:374594/143‑1 (MQ=255)                                                                                                                      gATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCAt                        <  1:633679/143‑1 (MQ=255)                                                                                                                         ggTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATggg                     <  2:45341/143‑1 (MQ=255)                                                                                                                          gTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAACGGTTAGCTTCCAATTACTGGTGTATACCCCCTGGgc                    >  1:652781/1‑143 (MQ=255)                                                                                                                            aaCGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgcag                  <  1:738304/143‑1 (MQ=255)                                                                                                                               gTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAc               >  1:963171/1‑143 (MQ=255)                                                                                                                                  cccAAAGGCATTTGCCTCTACACGACCGATTAAAATTCGCGTAATTCCCTTGACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTc            <  1:775870/143‑1 (MQ=255)                                                                                                                                   ccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGCGCTTCCAATTACTAGTGTATACACCCTGGGCAGCACCTCg           >  2:89223/1‑143 (MQ=255)                                                                                                                                     aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAGTTACTGGTGGATACACCCTGGGGAGCGCCTCGTc         >  2:1022610/1‑143 (MQ=255)                                                                                                                                          cATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAg    >  1:943794/1‑143 (MQ=255)                                                                                                                                           aTTTGCCTTTACACGCCTGATGAACCTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAgg   <  2:734052/143‑1 (MQ=255)                                                                                                                                            tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCATGTACTGGTGTATCCACCAGGGAGAGCACCTCGTCCGAAggg  >  2:6249/1‑143 (MQ=255)                                                                                                                                            tttGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAggg  <  1:214329/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                              ATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGG  >  NC_000913/3991195‑3991476

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