BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,772,254	T→G	5.8%	S351S (TCT→TCG)	phoB →	alkaline phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,772,254	0	T	G	5.8%	240.3 / 8.9	87	S351S (TCT→TCG)	phoB	alkaline phosphatase
Reads supporting (aligned to +/- strand): ref base T (47/34); new base G (2/3); total (50/37)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

TCGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTAA  >  CP000730/2772124‑2772390

tCGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACTAATGGAATCCTCAACCGATTAAATCTATGAAACAct                                                                                                                                >  2:1084983/1‑141 (MQ=255)

tCGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGAAAAGATAAAGGCTACGAATGAAATCCTAAACCGATTAAATCAAAGAAACAct                                                                                                                                >  1:209495/1‑141 (MQ=255)

tCGATGACGCCATTCAAGATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACAct                                                                                                                                <  1:604440/141‑1 (MQ=255)

 cGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAATGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACActc                                                                                                                               >  2:274408/1‑141 (MQ=255)

 cGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAAAAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAGAGCTATGAAACActc                                                                                                                               >  1:227222/1‑141 (MQ=255)

  gATGACGCCATTCAATATTCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACActct                                                                                                                              <  2:364363/141‑1 (MQ=255)

  gATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACActct                                                                                                                              <  2:531000/141‑1 (MQ=255)

    tGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTgg                                                                                                                            <  2:31606/141‑1 (MQ=255)

    tGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTgg                                                                                                                            <  1:726608/141‑1 (MQ=255)

      aCGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGAt                                                                                                                          >  2:134972/1‑141 (MQ=255)

       cGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATc                                                                                                                         >  2:546371/1‑141 (MQ=255)

       cGCCATTCAATATGCTAAGAAACATAAAGATACACGTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGATACACTCTGGATc                                                                                                                         <  2:116967/141‑1 (MQ=255)

        gCCATTCAATATGCTAAGAAACATAAAGATACACTTTTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCa                                                                                                                        <  1:341958/141‑1 (MQ=255)

        gCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCa                                                                                                                        <  2:444084/141‑1 (MQ=255)

        gCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCa                                                                                                                        <  1:410429/141‑1 (MQ=255)

        gCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCa                                                                                                                        <  1:1053644/141‑1 (MQ=255)

             tCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCatatat                                                                                                                   <  2:1080039/141‑1 (MQ=255)

                    gCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGaa                                                                                                            >  2:262231/1‑141 (MQ=255)

                        aGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGaaaaaa                                                                                                        >  1:889253/1‑141 (MQ=255)

                         gAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAAt                                                                                                       >  1:356679/1‑141 (MQ=255)

                           aaCATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATgg                                                                                                     <  1:160428/141‑1 (MQ=255)

                                aaaGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTaaa                                                                                                <  1:864083/141‑1 (MQ=255)

                                 aaGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAATGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATACATGACAGAAAAAATGGCTAAAg                                                                                               >  2:22367/1‑141 (MQ=255)

                                  aGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAgg                                                                                              >  2:57041/1‑141 (MQ=255)

                                  aGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAgg                                                                                              >  2:883223/1‑141 (MQ=255)

                                     tACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTGTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGa                                                                                           >  2:934716/1‑141 (MQ=255)

                                        acTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGa                                                                                        >  2:511993/1‑141 (MQ=255)

                                         cTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGAc                                                                                       <  2:65860/141‑1 (MQ=255)

                                           tGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGAccc                                                                                     <  1:211690/141‑1 (MQ=255)

                                            gTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCt                                                                                    <  1:134972/141‑1 (MQ=255)

                                              cGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGa                                                                                  <  2:113939/141‑1 (MQ=255)

                                                ttGCAACAGCAGATCACTCTACTGGTGGTCTAACAATGGGTAAAGATAAAGGCTACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaa                                                                                >  1:23110/1‑141 (MQ=255)

                                                 tGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaa                                                                               >  2:410680/1‑141 (MQ=255)

                                                 tGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaa                                                                               >  1:245484/1‑141 (MQ=255)

                                                 tGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaa                                                                               >  1:271064/1‑141 (MQ=255)

                                                 tGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaa                                                                               <  2:971337/141‑1 (MQ=255)

                                                  gCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGaaaaa                                                                              >  2:1024897/1‑141 (MQ=255)

                                                   cAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAg                                                                             <  1:380754/141‑1 (MQ=255)

                                                    aaCAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGt                                                                            <  2:566682/141‑1 (MQ=255)

                                                     aCAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGtt                                                                           >  2:107883/1‑141 (MQ=255)

                                                     aCAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATAGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGAAAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGtt                                                                           <  2:980633/141‑1 (MQ=255)

                                                             tCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGa                                                                   >  1:1027791/1‑141 (MQ=255)

                                                             tCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGa                                                                   >  1:91724/1‑141 (MQ=255)

                                                             tCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGa                                                                   >  2:520232/1‑141 (MQ=255)

                                                              cACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGaa                                                                  >  2:937563/1‑141 (MQ=255)

                                                                ctctACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAgg                                                                >  2:682370/1‑141 (MQ=255)

                                                                 tctACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGt                                                               <  1:870956/141‑1 (MQ=255)

                                                                   tACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTa                                                             <  1:692556/141‑1 (MQ=255)

                                                                     cTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATg                                                           >  2:871346/1‑141 (MQ=255)

                                                                      tggaggTGTAACAATAGCAAAAGATAAAGAAGACGAAGGGAATCCTCAACCGATTAAATCTATGAAACACTCTGAATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATgg                                                          <  2:378186/141‑1 (MQ=255)

                                                                       ggtggtCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCAGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGTAGGTTATGGa                                                         <  2:895075/141‑1 (MQ=255)

                                                                         tggtCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGAtt                                                       <  1:160803/141‑1 (MQ=255)

                                                                          ggtCTAACAATTGGTAAAGATAACGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTACAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTCTGGAttt                                                      >  1:765685/1‑141 (MQ=255)

                                                                          ggtCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGAttt                                                      >  2:139021/1‑141 (MQ=255)

                                                                          ggtCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTAAGAAACACTCGGGACCATATAGGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGAttc                                                      >  1:647773/1‑140 (MQ=255)

                                                                            tCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGa                                                    <  1:581345/141‑1 (MQ=255)

                                                                               aaCAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAAtt                                                 >  2:822839/1‑141 (MQ=255)

                                                                               aaCAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAAtt                                                 >  1:705723/1‑141 (MQ=255)

                                                                                 cAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAAtttt                                               >  1:218742/1‑141 (MQ=255)

                                                                                 cAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCGATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGCTTTGAAtttt                                               >  2:119942/1‑141 (MQ=255)

                                                                                     tGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTcaca                                           >  1:147080/1‑141 (MQ=255)

                                                                                       gTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACACTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAg                                         >  2:436124/1‑141 (MQ=255)

                                                                                           aGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatga                                     <  1:100944/141‑1 (MQ=255)

                                                                                           aGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatga                                     <  2:808396/141‑1 (MQ=255)

                                                                                           aGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatga                                     <  1:1435/141‑1 (MQ=255)

                                                                                           aGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatga                                     <  1:713554/141‑1 (MQ=255)

                                                                                            gATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAgatgat                                    >  1:226416/1‑141 (MQ=255)

                                                                                                aaGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGa                                >  1:115619/1‑141 (MQ=255)

                                                                                                  ggATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGTACTCGCGGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACGAGATGATATGaaa                              >  1:19797/1‑141 (MQ=255)

                                                                                                   gATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGaaaa                             >  1:582889/1‑141 (MQ=255)

                                                                                                   gATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGaaaa                             >  2:150802/1‑141 (MQ=255)

                                                                                                     tACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAg                           >  1:532939/1‑141 (MQ=255)

                                                                                                       cGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGt                         <  2:80760/141‑1 (MQ=255)

                                                                                                       cGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGt                         >  2:1030650/1‑141 (MQ=255)

                                                                                                        gAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGtt                        >  1:164992/1‑141 (MQ=255)

                                                                                                            ggAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTaaaa                    >  1:696970/1‑141 (MQ=255)

                                                                                                               aTCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCAAAAGGTGAAGACCCTGAAAAAGATATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAag                 >  2:607996/1‑141 (MQ=255)

                                                                                                               aTCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATCACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAag                 >  1:572806/1‑141 (MQ=255)

                                                                                                                 ccTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAagag               <  1:678967/141‑1 (MQ=255)

                                                                                                                      aCCGATTAAATCTATGAAACACTCTGGATAATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATaa          <  2:824889/141‑1 (MQ=255)

                                                                                                                           ttAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGc     <  1:57041/141‑1 (MQ=255)

                                                                                                                           ttAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGc     >  2:802091/1‑141 (MQ=255)

                                                                                                                             aaaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTa   >  1:994129/1‑141 (MQ=255)

                                                                                                                              aaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTaa  <  2:889253/141‑1 (MQ=255)

                                                                                                                              aaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTaa  <  2:387531/141‑1 (MQ=255)

                                                                                                                              aaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTaa  >  1:492218/1‑141 (MQ=255)

                                                                                                                              aaTCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAACGACACTTATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTaa  >  1:565585/1‑141 (MQ=255)

TCGATGACGCCATTCAATATGCTAAGAAACATAAAGATACACTTGTCGTTGCAACAGCAGATCACTCTACTGGTGGTCTAACAATTGGTAAAGATAAAGGATACGAATGGAATCCTCAACCGATTAAATCTATGAAACACTCTGGATCATATATGACAGAAAAAATGGCTAAAGGTGAAGACCCTGAAAAAGTTATTAATGAAGGTTATGGATTTGAATTTTCACAAGATGATATGAAAAAGGTTAAAAAAGAGAATAAAAAGCTAA  >  CP000730/2772124‑2772390

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