BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,135,356	T→C	5.2%	Q339R (CAG→CGG) ‡	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,135,356	0	T	C	5.2%	406.6 / 8.9	134	Q339R (CAG→CGG) ‡	pitB	phosphate transporter
Reads supporting (aligned to +/- strand): ref base T (70/57); new base C (2/5); total (72/62)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

Only 100 of 149 total aligned reads displayed.

CCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGA  >  NC_000913/3135207‑3135505

 cATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctg                                                                                                                                                      <  2:1067009/150‑1 (MQ=255)

 cATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctg                                                                                                                                                      <  2:1129574/150‑1 (MQ=255)

   tgatgaTCCACACCGGAGTGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                                                                    >  2:88654/1‑150 (MQ=255)

   tgatgaTCCACACCGGAGTGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                                                                    >  2:1337816/1‑150 (MQ=255)

   tgatgaTCCACACCGGAGGGTACTCAATGGTGCTTAACATATCGCTGCGAAGGTTTTTCAGCAGGTTCTGGTCTTCTTTTCTGACGACTGGCAGTTTCGCTAGCTGCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                                                                    <  2:1064115/150‑1 (MQ=255)

   tgatgaTCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGCCGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                                                                    >  1:1752796/1‑150 (MQ=255)

   tgatgaTCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                                                                    >  2:1317497/1‑150 (MQ=255)

   tgatgaTCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAgctggc                                                                                                                                                    >  2:181101/1‑150 (MQ=255)

         tCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCt                                                                                                                                              <  2:931868/150‑1 (MQ=255)

         tCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgcc                                                                                                                                              <  2:579112/150‑2 (MQ=255)

          ccACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccg                                                                                                                                             <  1:768659/150‑3 (MQ=255)

           cacaCCGGAGCGTACTCAATTGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCCAGCTgcgccgc                                                                                                                                            <  2:2011846/150‑4 (MQ=255)

           cacaCCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGa                                                                                                                                            >  1:676279/1‑150 (MQ=255)

           cacaCCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGa                                                                                                                                            >  1:1817049/1‑150 (MQ=255)

            acacCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGAc                                                                                                                                           <  1:1725498/150‑1 (MQ=255)

                   aGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCt                                                                                                                                    <  1:1890189/150‑1 (MQ=255)

                   aGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCt                                                                                                                                    <  2:678176/150‑1 (MQ=255)

                    gCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCtt                                                                                                                                   >  1:1952611/1‑150 (MQ=255)

                        aCTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTATTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACg                                                                                                                               >  2:1156787/1‑150 (MQ=255)

                         cTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACgg                                                                                                                              <  1:1615406/150‑1 (MQ=255)

                           cAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCt                                                                                                                            >  2:1857969/1‑150 (MQ=255)

                                 tGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAc                                                                                                                      <  1:1379887/150‑1 (MQ=255)

                                  gCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACt                                                                                                                     >  1:1084851/1‑150 (MQ=255)

                                      aaCATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTcc                                                                                                                 >  2:926545/1‑150 (MQ=255)

                                        cATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCat                                                                                                               >  1:183554/1‑150 (MQ=255)

                                        cATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTACat                                                                                                               >  2:1818825/1‑150 (MQ=255)

                                         atagCGCTGCGAAGTTTTTTCTGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCata                                                                                                              <  1:128674/150‑1 (MQ=255)

                                           atCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATAtt                                                                                                            >  2:229452/1‑150 (MQ=255)

                                            tCGCTGCAAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTg                                                                                                           >  2:380075/1‑150 (MQ=255)

                                             cGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGc                                                                                                          <  2:1760711/150‑1 (MQ=255)

                                               cTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCt                                                                                                        >  2:1648542/1‑150 (MQ=255)

                                                tGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTg                                                                                                       <  2:1717138/150‑1 (MQ=255)

                                                ggCGAAGTTTATACAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTg                                                                                                       <  1:83314/149‑1 (MQ=255)

                                                 gCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTgg                                                                                                      >  1:152388/1‑150 (MQ=255)

                                                  cGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgc                                                                                                     <  1:2017094/150‑1 (MQ=255)

                                                     aGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgcagc                                                                                                  >  1:501285/1‑150 (MQ=255)

                                                     aGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGgcagc                                                                                                  <  1:1378852/150‑1 (MQ=255)

                                                      gTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCAGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGTCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACATTCCATATTGCCTGgcagca                                                                                                 >  1:84602/1‑150 (MQ=255)

                                                       gtttttCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACTGCTCGTAACTTTCCATATTGCCTGgcagcct                                                                                                >  1:124018/2‑148 (MQ=255)

                                                            tccgcagGTTCTGGGCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGttt                                                                                           <  2:986416/147‑1 (MQ=255)

                                                             cagcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGtttt                                                                                          <  2:1766921/150‑1 (MQ=255)

                                                              tgcagGTGCTGGTCTTCTTTACTGACGCCTGGCTGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCGCGTAACTTTCCATATTGCCTGGCAGAATCGTTTTa                                                                                         <  2:1090759/149‑1 (MQ=255)

                                                              agcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTa                                                                                         >  1:664640/1‑150 (MQ=255)

                                                              agcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCT‑GCAGCATCgactgat                                                                                        >  2:1670916/1‑144 (MQ=255)

                                                               gcagGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCTCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCAGCAACGTTTTaa                                                                                        >  1:417825/1‑150 (MQ=255)

                                                                  ggTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgc                                                                                     >  1:199420/1‑150 (MQ=255)

                                                                     tCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgcgcg                                                                                  >  1:1514523/1‑150 (MQ=255)

                                                                     tCTGGTCTTCTTTACTGACGCCGGGCAGTTTCGCTAGCTTCGCGGATGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAcgcgcg                                                                                  <  2:1957485/150‑1 (MQ=255)

                                                                       tGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCa                                                                                <  2:777346/150‑1 (MQ=255)

                                                                       tGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCa                                                                                <  2:132052/150‑1 (MQ=255)

                                                                        ggTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCCCGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCaa                                                                               >  2:1267744/1‑150 (MQ=255)

                                                                           cttcttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAGCGCGCGCAATAg                                                                            >  2:1350209/1‑150 (MQ=255)

                                                                            ttcttTACTGCCGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGc                                                                           <  1:1775524/150‑1 (MQ=255)

                                                                                ttACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCa                                                                       <  2:14695/150‑1 (MQ=255)

                                                                                   cTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAg                                                                    >  1:220912/1‑150 (MQ=255)

                                                                                   cTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAg                                                                    <  1:497252/150‑1 (MQ=255)

                                                                                    tGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGGTTAACGGCTCGTAACTTTACATAATGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAgg                                                                   >  2:38834/1‑150 (MQ=255)

                                                                                    tGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAgg                                                                   >  2:1523407/1‑150 (MQ=255)

                                                                                     gACGCCTGGCAGTATCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATAGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGt                                                                  <  2:1212888/150‑1 (MQ=255)

                                                                                     gACGCCTCGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTTACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCATTAGCATCAAAGGt                                                                  >  1:42401/1‑150 (MQ=255)

                                                                                            ggCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCGTATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGcc                                                           <  2:1565041/150‑1 (MQ=255)

                                                                                            ggCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGcc                                                           <  1:751037/150‑1 (MQ=255)

                                                                                             gCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCg                                                          >  1:840910/1‑150 (MQ=255)

                                                                                              cAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCgg                                                         <  1:1317497/150‑1 (MQ=255)

                                                                                              cAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCgg                                                         >  1:1672664/1‑150 (MQ=255)

                                                                                                gTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGAt                                                       >  2:524898/1‑150 (MQ=255)

                                                                                                gTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGAt                                                       <  1:102603/150‑1 (MQ=255)

                                                                                                    cGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACa                                                   <  1:1016105/150‑1 (MQ=255)

                                                                                                     gCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAg                                                  >  2:1814085/1‑150 (MQ=255)

                                                                                                      cTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGt                                                 <  1:1337816/150‑1 (MQ=255)

                                                                                                      cTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGt                                                 <  1:88654/150‑1 (MQ=255)

                                                                                                           ttCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAAc                                            <  1:793896/150‑1 (MQ=255)

                                                                                                           ttCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAAc                                            <  1:911537/150‑1 (MQ=255)

                                                                                                              gcgGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAActct                                         >  2:17086/1‑150 (MQ=255)

                                                                                                                ggAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGt                                       <  2:894032/150‑1 (MQ=255)

                                                                                                                  aGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTa                                     <  2:14482/150‑1 (MQ=255)

                                                                                                                   ggTATCAGAGATGCACAGCATATTGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTAc                                    <  2:1022680/150‑1 (MQ=255)

                                                                                                                        cAGAGATGCACTGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCGCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTACCGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCg                               >  2:495924/1‑150 (MQ=255)

                                                                                                                           agaTGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGc                            <  2:968024/150‑1 (MQ=255)

                                                                                                                             aTGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCa                          <  2:579663/150‑1 (MQ=255)

                                                                                                                              tGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCCCGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCAt                         >  1:1667517/1‑150 (MQ=255)

                                                                                                                              tGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCAt                         <  2:1416558/150‑1 (MQ=255)

                                                                                                                                cacaGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCa                       <  1:229452/150‑1 (MQ=255)

                                                                                                                                    gCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCg                   <  2:676279/150‑1 (MQ=255)

                                                                                                                                     cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGa                  <  2:1955382/150‑1 (MQ=255)

                                                                                                                                      aTAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAt                 >  2:1386100/1‑150 (MQ=255)

                                                                                                                                      aTAATGAGGCGCAGCTGGCTGCGCTGACTCTCGCTTAACGGGTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAt                 <  2:850165/150‑1 (MQ=255)

                                                                                                                                        aaTGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAtgc               >  1:1211121/1‑150 (MQ=255)

                                                                                                                                            cggcgCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCACTAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCGGCa           >  1:1871155/1‑150 (MQ=255)

                                                                                                                                            cggcgCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCa           <  1:983106/150‑1 (MQ=255)

                                                                                                                                            cggcgCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGGTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAtgctgcc           >  1:1983275/1‑149 (MQ=255)

                                                                                                                                            cagcgcAGCCAGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTGTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCa           >  1:1467165/3‑150 (MQ=255)

                                                                                                                                             ggcgCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTGGAACGCGCGCAGTAGCGTCAAAGGTAGTTGCCGGATGAGAGGGAAACTCGGTTACTGGCGTGCCAACAGTCGATGCTGCAg          >  1:369567/1‑150 (MQ=255)

                                                                                                                                              gcgcAGCTGGCTGCGCTGACTCACGCTTACCGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCTAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCGATCAGTCGATGCTGCCgg         >  2:933744/1‑150 (MQ=255)

                                                                                                                                              gcgcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAgg         <  1:181101/150‑1 (MQ=255)

                                                                                                                                                   gctggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATg    <  2:1869540/150‑1 (MQ=255)

                                                                                                                                                    ctggccgCGCTGACTCACGCTTAACGGCTCGTAACTTACCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATgg   <  2:769163/150‑1 (MQ=255)

                                                                                                                                                     tggttgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAATGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGa  <  1:1212913/150‑1 (MQ=255)

                                                                                                                                                     tggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGa  <  1:380075/150‑1 (MQ=255)

                                                                                                                                                     tggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGa  <  2:183554/150‑1 (MQ=255)

CCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGA  >  NC_000913/3135207‑3135505

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