BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,316,701:1	+T	7.2%	coding (937/1209 nt)	aspC →	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,316,701	1	.	T	7.2%	926.8 / 66.5	222	coding (937/1209 nt)	aspC	aminotransferase
Reads supporting (aligned to +/- strand): ref base . (77/129); new base T (5/11); total (82/140)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

GCAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGT  >  NC_002947/4316552‑4316849

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

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

      tCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTgg                                                                                                                                                >  2:114504/1‑150 (MQ=255)

       cAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGa                                                                                                                                               <  1:2962418/150‑1 (MQ=255)

       cAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGa                                                                                                                                               <  1:3690414/150‑1 (MQ=255)

       cAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGa                                                                                                                                               <  1:1827329/150‑1 (MQ=255)

       cAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGa                                                                                                                                               <  2:3148468/150‑1 (MQ=255)

       cAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGC‑CATGTGCCTGGa                                                                                                                                               <  2:2906332/150‑1 (MQ=255)

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

          ccTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGgt                                                                                                                                            >  1:3964904/1‑150 (MQ=255)

          ccTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGgt                                                                                                                                            <  1:347752/150‑1 (MQ=255)

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

            tgctgcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGgtgt                                                                                                                                          <  1:196851/150‑1 (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:4327020/1‑150 (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)

                gcATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTaaa                                                                                                                                      <  2:2141736/150‑1 (MQ=255)

                 cATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAACGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAAc                                                                                                                                     <  1:3945167/150‑1 (MQ=255)

                   tGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGc                                                                                                                                   <  1:2718103/150‑1 (MQ=255)

                   tGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGc                                                                                                                                   >  2:3978346/1‑150 (MQ=255)

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

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

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

                                 cgcgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCt                                                                                                                     <  1:811328/150‑1 (MQ=255)

                                   cgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCGGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCttt                                                                                                                   <  1:732162/150‑1 (MQ=255)

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

                                   cgTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCttt                                                                                                                   >  1:3557674/1‑150 (MQ=255)

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

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

                                     tATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGt                                                                                                                 >  1:565549/1‑150 (MQ=255)

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

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

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

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

                                             tCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGc                                                                                                         >  2:1860129/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                       gACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGc                                                                                               <  1:3227411/150‑1 (MQ=255)

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

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

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

                                                           tGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCt                                                                                           <  1:2853227/150‑1 (MQ=255)

                                                           tGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCt                                                                                           <  1:2936255/150‑1 (MQ=255)

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

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

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

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

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

                                                                       gacCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGc                                                                               <  2:3602179/150‑1 (MQ=255)

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

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

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

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

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

                                                                                   gCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACAGTCGcc                                                                   >  1:3340385/1‑150 (MQ=255)

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

                                                                                     cGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCAc                                                                 >  1:346191/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                          aGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATa                                            <  2:3212327/150‑1 (MQ=255)

                                                                                                          aGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATa                                            >  2:1393372/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                gTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATc                      >  1:2029589/1‑150 (MQ=255)

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

                                                                                                                                     acacctcCAGGC‑CATGTGCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGAc                 >  1:3611259/5‑150 (MQ=255)

                                                                                                                                       ttCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGa               <  1:3869055/150‑1 (MQ=255)

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

                                                                                                                                         ctcCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGAtt             >  2:1328739/1‑150 (MQ=255)

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

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

                                                                                                                                          tcCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTc            <  1:953807/150‑1 (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)

                                                                                                                                               gCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTAt       <  1:2303061/150‑1 (MQ=255)

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

                                                                                                                                                    tG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGt  <  2:2682520/150‑1 (MQ=255)

                                                                                                                                                    tG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGt  <  2:2578949/150‑1 (MQ=255)

GCAAGCTCAGCCTGCTGCATGAAAGTGGGCAGCCGCGTATTTTACTCGATGCCGTGACCTGCATTCTCGACGACCGGGGCTTCGCCGAACGCCTTTGCCAGTGCATAGCTCAGCGACGTGCGGTAATGGTCGATATTCTCCAGGCACATG‑GCCTGGAGGTGTTAAACGCTGACGGCACATTCTTTGTTTACGTGCGATGTCCTGCGGCTGTGAGCAACGGCAGGACATTCGCCACTGCCAGGGACTTTTCGCAATACCTGGCCACTGAGCTGGGTATCCTGACGATTCCCTATCAGGT  >  NC_002947/4316552‑4316849

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