BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,564,973	A→T	8.5%	L133Q (CTG→CAG)	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,564,973	0	A	T	8.5%	256.2 / 10.6	92	L133Q (CTG→CAG)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand): ref base A (39/45); new base T (5/3); total (44/48)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 101 total aligned reads displayed.

GCGGGCGCCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATT  >  NC_002947/1564824‑1565112

gCGGGCGCCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTggcca                                                                                                                                             <  2:868322/150‑1 (MQ=255)

gCGGGCGCCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTggcca                                                                                                                                             <  2:1269132/150‑1 (MQ=255)

 cGGGCGCCACCATTTCAGGCCGATCAGCACCAGTGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTggccag                                                                                                                                            <  1:2051087/150‑1 (MQ=255)

 cGGGCGCCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTggccag                                                                                                                                            <  1:1391818/150‑1 (MQ=255)

  gggCGCCACCATTTCAGGCCGATCAGCACCAGGGCAGCAGCGGCGAAGCCGGCCATTGGCGAGACCAGCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTggccagg                                                                                                                                           >  2:2156175/1‑150 (MQ=255)

    gcgcCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTggccaggcc                                                                                                                                         >  2:679118/1‑150 (MQ=255)

    gcgcCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTggccaggcc                                                                                                                                         >  2:1183004/1‑150 (MQ=255)

      gcCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCaa                                                                                                                                       <  1:847547/150‑1 (MQ=255)

       ccaccaTTTCAGGCCGATCAGAACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAac                                                                                                                                      <  2:1061745/150‑1 (MQ=255)

        caccaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGGGAGCACCAGCGACATGGCGATATCGATCGCTTTCTGACAGTTGACGCCCTCGCCCAGCGGAATGTCGTTGATCATGGCGTTGGCCAGGCCAaca                                                                                                                                     >  2:2304580/1‑150 (MQ=255)

        caccaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGTGTTGGCCTGGCCAacg                                                                                                                                     >  2:171470/1‑149 (MQ=255)

        caccaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAaca                                                                                                                                     >  2:1237489/1‑150 (MQ=255)

        caccaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAaca                                                                                                                                     >  1:1085814/1‑150 (MQ=255)

         accaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAacac                                                                                                                                    >  2:1155609/1‑150 (MQ=255)

         accaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAacac                                                                                                                                    >  1:2242016/1‑150 (MQ=255)

             tttCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                <  2:1327067/150‑1 (MQ=255)

               tCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGATCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGa                                                                                                                              <  2:1430454/150‑1 (MQ=255)

               tCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGa                                                                                                                              <  1:1183004/150‑1 (MQ=255)

                 aGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATg                                                                                                                            >  2:2091278/1‑150 (MQ=255)

                  ggCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGGGACCACCAGCGACATGGCGATATCGATCGCTTTATGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCTAGGCCAACATCGAGGATgg                                                                                                                           >  1:514067/1‑150 (MQ=255)

                  ggCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATgg                                                                                                                           >  2:599180/1‑150 (MQ=255)

                   gCCGATCAGCACCAGGGCAGCAAAGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGa                                                                                                                          >  1:642707/1‑150 (MQ=255)

                      gATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGACGGTATCGATCGCTTTCTGCGAGTTGACGCCGTCGCCCAGCGGAATATCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGCAAcc                                                                                                                       >  2:444289/1‑150 (MQ=255)

                       aTCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGTGTTGGCCTGGCCAACGcc                                                                                                                                   >  2:2254998/1‑137 (MQ=255)

                       aTCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGTGTTGGCCTGGCCAACGcc                                                                                                                                   <  1:2254998/137‑1 (MQ=255)

                        tCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGa                                                                                                                     <  1:1155609/150‑1 (MQ=255)

                          aGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATc                                                                                                                   >  2:1325291/1‑150 (MQ=255)

                           gCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCa                                                                                                                  <  2:2028765/150‑1 (MQ=255)

                             aCCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGc                                                                                                                <  1:679118/150‑1 (MQ=255)

                             aCCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGc                                                                                                                >  1:2081709/1‑150 (MQ=255)

                             aCCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAATTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGc                                                                                                                >  1:883038/1‑150 (MQ=255)

                                aGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCgtg                                                                                                             >  2:244498/1‑150 (MQ=255)

                                 gggCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCgtgt                                                                                                            <  1:2202070/150‑1 (MQ=255)

                                 gggCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAATTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCgtgt                                                                                                            <  2:883038/150‑1 (MQ=255)

                                       caACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTTGCCAGGCCAACACCGAGGATGGAAc                                                                                                                        <  1:1251973/132‑1 (MQ=255)

                                       caACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTTGCCAGGCCAACACCGAGGATGGAAc                                                                                                                        >  2:1251973/1‑132 (MQ=255)

                                           ggCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGa                                                                                                  >  2:2025827/1‑150 (MQ=255)

                                                aGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCgg                                                                                             <  2:2333794/150‑1 (MQ=255)

                                                         ttGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGAAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCaa                                                                                    >  1:1654609/1‑150 (MQ=255)

                                                         ttGGCGAGACCACCAGCGACATGGCCCGAGCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCaa                                                                                    >  1:1230039/1‑150 (MQ=255)

                                                          tGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCaaa                                                                                   <  2:642707/150‑1 (MQ=255)

                                                               agaCCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTAcc                                                                              <  2:2151967/150‑1 (MQ=255)

                                                                 accaccAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAg                                                                            >  1:1669218/1‑150 (MQ=255)

                                                                  ccaccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTAGCCCGGCTGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAgg                                                                           <  2:1230039/150‑1 (MQ=255)

                                                                   caccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                >  1:864870/1‑96 (MQ=255)

                                                                   caccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                <  2:864870/96‑1 (MQ=255)

                                                                   caccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAAc                                                                                                                        >  1:749650/1‑104 (MQ=255)

                                                                   caccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAAc                                                                                                                        <  2:749650/104‑1 (MQ=255)

                                                                    accaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTg                                                                         >  1:903762/1‑150 (MQ=255)

                                                                          gACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAgg                                                                   <  1:1325291/150‑1 (MQ=255)

                                                                           aCATGGCGATATCGATCGCTTTCTGCCAGTTCACGTCGTCGCCCGGCGGAATGGCGTTGATCGGGGCGTTGGCCAGGCAAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGt                                                                  <  1:2304580/150‑1 (MQ=255)

                                                                             aTGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGCTGCCACCGTACTCGCTACCCTGGTTc                                                                >  2:501015/1‑150 (MQ=255)

                                                                              tGGCGATATCGCCCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGCCGTTCTCCAGGCTCACAAAGAGGATGGAACCGATCAGCGTGTGAGAACTGTAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTcc                                                               <  1:1333055/150‑1 (MQ=255)

                                                                              tGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTcc                                                               >  1:897277/1‑150 (MQ=255)

                                                                                gCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAATTACCAGGTGCCCAGGTTCCAg                                                             <  1:1237489/150‑1 (MQ=255)

                                                                                  gATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGt                                                           <  1:2091278/150‑1 (MQ=255)

                                                                                  atcttcgatcgCCCTCTGAGAGGTGACCCCCCCGCCCAGCTGATTGTCTTTGAGCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGt                                                           <  2:1668121/146‑1 (MQ=255)

                                                                                      tcgatcgCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATg                                                       <  1:501015/150‑1 (MQ=255)

                                                                                       cgatcgCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGagat                                                      >  1:1668121/1‑147 (MQ=255)

                                                                                        gatcgCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGc                                                     >  2:1442298/1‑150 (MQ=255)

                                                                                              tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCa                                                              >  2:2147808/1‑135 (MQ=255)

                                                                                              tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCa                                                              <  1:2147808/135‑1 (MQ=255)

                                                                                                tCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGcagcca                                             <  2:897277/150‑1 (MQ=255)

                                                                                                  tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                >  2:481753/1‑65 (MQ=255)

                                                                                                  tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                <  1:481753/65‑1 (MQ=255)

                                                                                                  tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTAcc                                                                              <  2:1050851/115‑1 (MQ=255)

                                                                                                  tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTAcc                                                                              >  1:1050851/1‑115 (MQ=255)

                                                                                                  tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGcagccagc                                           >  2:1934528/1‑150 (MQ=255)

                                                                                                       gTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACga                                      <  1:244498/150‑1 (MQ=255)

                                                                                                          gACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTg                                                                         <  2:977697/112‑1 (MQ=255)

                                                                                                          gACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTg                                                                         >  1:977697/1‑112 (MQ=255)

                                                                                                          gACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGagaa                                   <  2:546334/150‑1 (MQ=255)

                                                                                                            cGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGagaaga                                 <  1:2025827/150‑1 (MQ=255)

                                                                                                             gCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGAGCAGCGTGTGCGAACTGGAGGCGGGCATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGAc                                <  1:315114/150‑1 (MQ=255)

                                                                                                               cGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGc                                                                                       <  2:668709/93‑1 (MQ=255)

                                                                                                               cGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGc                                                                                       >  1:668709/1‑93 (MQ=255)

                                                                                                                gTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCAt                             >  1:1643889/1‑150 (MQ=255)

                                                                                                                 tCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATg                            <  2:2081709/150‑1 (MQ=255)

                                                                                                                 tCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGCATGGAACCGCTTACCGCGCACGATCTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATg                            <  1:2156175/150‑1 (MQ=255)

                                                                                                                   gCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGc                          <  1:599180/150‑1 (MQ=255)

                                                                                                                   gCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGc                          >  2:1000245/1‑150 (MQ=255)

                                                                                                                            aaTGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTg                 >  2:2164167/1‑150 (MQ=255)

                                                                                                                            aaTGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCACGGCCTGGCCATg                 <  1:1749556/150‑1 (MQ=255)

                                                                                                                                tCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTACGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATg                            >  1:1232949/1‑135 (MQ=255)

                                                                                                                                tCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTACGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATg                            <  2:1232949/135‑1 (MQ=255)

                                                                                                                                ctgTTGATCAGATCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCg             <  2:722235/148‑1 (MQ=255)

                                                                                                                                  gTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACTATGGCCAGGCCGTGGccatg           <  1:968341/150‑4 (MQ=255)

                                                                                                                                  gTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGgt           >  2:1862083/1‑150 (MQ=255)

                                                                                                                                   ttGGTCAGGGCGGTGGCAAGGCCAACACCGAGGATGGAACGGATCAGCGTGTGGGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCTAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGgtg          <  2:2242016/150‑1 (MQ=255)

                                                                                                                                      aTCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGAt                                                        >  1:90365/1‑101 (MQ=255)

                                                                                                                                      aTCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTc       <  1:1000245/150‑1 (MQ=255)

                                                                                                                                      aTCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGACCAGGTTCCAGGTGAt                                                        <  2:90365/101‑1 (MQ=255)

                                                                                                                                       tCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGATCCAAGGTACCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTca      >  1:1795643/1‑150 (MQ=255)

                                                                                                                                        ctcggtgTTGGCCTGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCACGTGCCCAGGTTCGAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTcac     >  2:665309/7‑150 (MQ=255)

                                                                                                                                         tcggtgTTGGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACga                                      >  1:2081266/6‑116 (MQ=255)

                                                                                                                                         tcggtgTTGGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACga                                      <  2:2081266/111‑1 (MQ=255)

                                                                                                                                           ggtgTTGGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAg                                                             >  2:736946/4‑91 (MQ=255)

                                                                                                                                           ggtgTTCGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAg                                                             <  1:736946/88‑1 (MQ=255)

                                                                                                                                           ggCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACAtt  <  1:1934528/150‑1 (MQ=255)

                                                                                                                                                  gccaggccaACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGagaa                                   <  1:999903/110‑1 (MQ=255)

GCGGGCGCCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATT  >  NC_002947/1564824‑1565112

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