BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,316,697	Δ1 bp	7.4%	coding (933/1209 nt)	aspC →	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,316,697	0	A	.	7.4%	750.0 / 50.4	217	coding (933/1209 nt)	aspC	aminotransferase
Reads supporting (aligned to +/- strand): ref base A (77/124); new base . (5/11); total (82/135)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 223 total aligned reads displayed.

CTTGGCAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCA  >  NC_002947/4316548‑4316846

cTTGGCAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGca                                                                                                                                                        >  2:3227411/1‑150 (MQ=255)

cTTGGCAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGca                                                                                                                                                        >  1:1742184/1‑150 (MQ=255)

    gCAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATg                                                                                                                                                    <  2:876268/150‑1 (MQ=255)

         cTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTg                                                                                                                                              <  2:377796/150‑1 (MQ=255)

             tccTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCTAGTGCATAGCTCAGCGACGTGCGGTAATGGTGGATATTCTCCAGGC‑CATGTGCCTGGAga                                                                                                                                          <  2:789695/149‑2 (MQ=255)

              ccTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGC‑CATGTGCCTGGAgaa                                                                                                                                         <  1:3205969/150‑3 (MQ=255)

               ctgctgCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAggat                                                                                                                                        <  2:1046192/150‑3 (MQ=255)

               ctgctgCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGgtg                                                                                                                                        <  1:3305455/150‑1 (MQ=255)

                tgctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGgtgt                                                                                                                                       >  2:35043/1‑150 (MQ=255)

                tgctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGgtgt                                                                                                                                       <  2:2441109/150‑1 (MQ=255)

                tgctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGC‑CATGTGCCTGGAgaata                                                                                                                                       <  2:2963977/150‑5 (MQ=255)

                 gctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGtt                                                                                                                                      <  1:1555024/150‑1 (MQ=255)

                  ctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTa                                                                                                                                     >  2:602872/1‑150 (MQ=255)

                  ctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGC‑CATGTGCCTGGAgaatatc                                                                                                                                     <  1:1545670/150‑7 (MQ=255)

                  ctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATCCTCCAGGC‑CATGTGCCTGGAgaatatc                                                                                                                                     <  1:168090/150‑7 (MQ=255)

                   tgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATCCTCCAGGC‑CATGTGCCTGGAgaatatcg                                                                                                                                    <  1:1132348/150‑8 (MQ=255)

                     cATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAAc                                                                                                                                  <  2:2602354/150‑1 (MQ=255)

                      aTGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACg                                                                                                                                 >  2:3293582/1‑150 (MQ=255)

                        gAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCt                                                                                                                               <  2:2204701/150‑1 (MQ=255)

                         aaaGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTg                                                                                                                              >  1:2879243/1‑150 (MQ=255)

                           aGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGAc                                                                                                                            <  2:1407720/150‑1 (MQ=255)

                           aGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTACATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGAc                                                                                                                            >  2:125090/1‑150 (MQ=255)

                               ggCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGca                                                                                                                        >  1:3648651/1‑150 (MQ=255)

                                gCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGcac                                                                                                                       <  1:867302/150‑1 (MQ=255)

                                       cgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCttt                                                                                                                <  2:4317946/150‑1 (MQ=255)

                                       cgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCttt                                                                                                                <  2:3129898/150‑1 (MQ=255)

                                        gTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTg                                                                                                               <  2:3980019/150‑1 (MQ=255)

                                         tATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGt                                                                                                              <  2:3648651/150‑1 (MQ=255)

                                          aTTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGtt                                                                                                             <  1:2366840/150‑1 (MQ=255)

                                           ttttACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGttt                                                                                                            <  2:3112372/150‑1 (MQ=255)

                                            tttACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTa                                                                                                           >  1:1150060/1‑150 (MQ=255)

                                               aCTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGt                                                                                                        >  2:1470524/1‑150 (MQ=255)

                                               aCTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGt                                                                                                        >  1:4118297/1‑150 (MQ=255)

                                                 tCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGc                                                                                                      <  2:4380775/150‑1 (MQ=255)

                                                 tCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGc                                                                                                      <  2:1150377/150‑1 (MQ=255)

                                                 tCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGc                                                                                                      <  2:3427152/150‑1 (MQ=255)

                                                 tCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGCTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGc                                                                                                      >  2:3764677/1‑150 (MQ=255)

                                                  cGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCg                                                                                                     <  2:562653/150‑1 (MQ=255)

                                                   gATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGa                                                                                                    <  1:4519471/150‑1 (MQ=255)

                                                   gATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGa                                                                                                    >  1:680318/1‑150 (MQ=255)

                                                    aTGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGAt                                                                                                   <  2:2051195/150‑1 (MQ=255)

                                                    aTGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGAt                                                                                                   >  2:347621/1‑150 (MQ=255)

                                                          tGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTACTGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATGCTTTGTTTACGTGCGATGTCCTg                                                                                             >  1:559900/1‑150 (MQ=255)

                                                              cTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGc                                                                                         >  1:245486/1‑150 (MQ=255)

                                                              cTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGc                                                                                         >  2:704430/1‑150 (MQ=255)

                                                              cTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGc                                                                                         >  2:2114897/1‑150 (MQ=255)

                                                              cTGCATTCTCGACGACCGGGGCTTCGCCGAACGCATTTGCCAGTGCATAGATCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGc                                                                                         >  2:4116868/1‑150 (MQ=255)

                                                               tGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCt                                                                                        >  1:1860274/1‑150 (MQ=255)

                                                               tGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCt                                                                                        >  1:3906005/1‑150 (MQ=255)

                                                                gCATTCTCGACGTCCGGTGCTTCTCCGAACGCCTTTGACAGTGCATAGCTCTGCGACGTGCTGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCAGATTCTTTGTTTACGTGCGATGTCCTGCGGCtg                                                                                       <  1:933055/150‑1 (MQ=255)

                                                                gCAGTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtg                                                                                       <  2:1329389/150‑1 (MQ=255)

                                                                 cATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCtgt                                                                                      <  2:2879183/150‑1 (MQ=255)

                                                                   ttCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGa                                                                                    <  2:209188/150‑1 (MQ=255)

                                                                   ttCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGa                                                                                    <  1:1389144/150‑1 (MQ=255)

                                                                   ttCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCAGAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGa                                                                                    <  1:2717042/150‑1 (MQ=255)

                                                                      tcGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCa                                                                                 >  1:1839682/1‑150 (MQ=255)

                                                                       cgacgaCCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCaa                                                                                <  2:1456866/150‑1 (MQ=255)

                                                                            acCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCa                                                                           <  1:3282703/150‑1 (MQ=255)

                                                                            acCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCa                                                                           <  1:3370858/150‑1 (MQ=255)

                                                                             ccGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAg                                                                          <  2:1069389/150‑1 (MQ=255)

                                                                              cGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAgg                                                                         <  1:3721339/150‑1 (MQ=255)

                                                                                gggCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGAc                                                                       >  2:2335821/1‑150 (MQ=255)

                                                                                    ttCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTc                                                                   <  2:1589346/150‑1 (MQ=255)

                                                                                    ttCGCCGAACGCCATTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTc                                                                   <  2:4258026/150‑1 (MQ=255)

                                                                                       gCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGcc                                                                >  1:3950105/1‑150 (MQ=255)

                                                                                        ccGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCa                                                               <  2:2299169/150‑1 (MQ=255)

                                                                                          gAACGCCATTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGATAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACt                                                             <  2:2263505/150‑1 (MQ=255)

                                                                                              gCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCa                                                         >  1:2179862/1‑150 (MQ=255)

                                                                                                 tttGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAggg                                                      >  2:3686120/1‑150 (MQ=255)

                                                                                                 tttGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAggg                                                      >  2:1311626/1‑150 (MQ=255)

                                                                                                 tttGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAggg                                                      >  1:3258350/1‑150 (MQ=255)

                                                                                                   tGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGAc                                                    <  2:2628245/150‑1 (MQ=255)

                                                                                                      cAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACttt                                                 >  1:1499781/1‑150 (MQ=255)

                                                                                                        gTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTc                                               <  1:276503/150‑1 (MQ=255)

                                                                                                         tGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCg                                              <  1:347621/150‑1 (MQ=255)

                                                                                                         tGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCg                                              <  2:1543244/150‑1 (MQ=255)

                                                                                                          gCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGc                                             >  1:594406/1‑150 (MQ=255)

                                                                                                               gCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATAc                                        <  1:7443/150‑1 (MQ=255)

                                                                                                                  cAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTg                                     <  2:3964904/150‑1 (MQ=255)

                                                                                                                  cAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTg                                     <  1:3263473/150‑1 (MQ=255)

                                                                                                                   aGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTgg                                    <  1:2686619/150‑1 (MQ=255)

                                                                                                                    gCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGc                                   <  1:244090/150‑1 (MQ=255)

                                                                                                                      gACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCa                                 <  2:319383/150‑1 (MQ=255)

                                                                                                                         gTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTg                              >  1:1998616/1‑150 (MQ=255)

                                                                                                                              gTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTg                         <  2:725691/150‑1 (MQ=255)

                                                                                                                              gTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTg                         <  1:35043/150‑1 (MQ=255)

                                                                                                                               tAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                        <  2:1856126/150‑1 (MQ=255)

                                                                                                                               tAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTgg                        <  1:903099/150‑1 (MQ=255)

                                                                                                                                        atatTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGa               <  2:4496059/150‑1 (MQ=255)

                                                                                                                                           acctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACt                                                             >  1:1219324/3‑101 (MQ=255)

                                                                                                                                           acctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGa            >  1:356733/3‑150 (MQ=255)

                                                                                                                                            tctcCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAt           >  2:2136827/1‑150 (MQ=255)

                                                                                                                                             ctcCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAtt          <  1:1677596/150‑1 (MQ=255)

                                                                                                                                             ctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAtt          >  1:2958924/1‑150 (MQ=255)

                                                                                                                                             ctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAtt          >  1:1173710/1‑150 (MQ=255)

                                                                                                                                              tcCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTc         >  2:1658672/1‑150 (MQ=255)

                                                                                                                                               ccAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTcc        <  2:198782/150‑1 (MQ=255)

                                                                                                                                                cAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTccc       >  1:4342254/1‑150 (MQ=255)

                                                                                                                                                 aGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCt      >  1:1027624/1‑150 (MQ=255)

                                                                                                                                                     acaTG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCa  <  1:4296196/150‑1 (MQ=255)

CTTGGCAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCA  >  NC_002947/4316548‑4316846

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