BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,338,677	T→C	7.5%	noncoding (353/395 nt)	PP_mr53 →	CrcZ ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,338,677	0	T	C	7.5%	336.2 / 13.0	120	noncoding (353/395 nt)	PP_mr53	CrcZ ncRNA
Reads supporting (aligned to +/- strand): ref base T (58/53); new base C (5/4); total (63/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-01

Only 100 of 120 total aligned reads displayed.

ACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATGCT  >  NC_002947/5338531‑5338825

aCATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGTTTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAg                                                                                                                                                   >  1:1294930/1‑150 (MQ=255)

aCATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAg                                                                                                                                                   >  1:2255785/1‑150 (MQ=255)

  aTCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGc                                                                                                                                                 >  2:1940098/1‑150 (MQ=255)

     aGCGGCCAAAGCAATCCGGTTGCTCTTTACCCCCGGTTTGGGGGTCGTAAACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTcc                                                                                                                                              <  2:15043/150‑1 (MQ=255)

      gCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAACCGAACTAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTccc                                                                                                                                             >  1:23575/1‑150 (MQ=255)

      gCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTccc                                                                                                                                             <  2:1950318/150‑1 (MQ=255)

         gCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCtt                                                                                                                                          >  1:1605019/1‑150 (MQ=255)

          ccAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTAGGGGAGCCTAAGCTccccaag                                                                                                                                         <  2:1029004/150‑4 (MQ=255)

           cAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGt                                                                                                                                        <  1:57077/150‑1 (MQ=255)

           cAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGt                                                                                                                                        >  2:2188857/1‑150 (MQ=255)

            aaaGTAATAAGTTTGCTCTTTAACCCCGTTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTc                                                                                                                                       <  2:1961977/150‑1 (MQ=255)

            aaaGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAg                                                                                                                                                   <  1:1140192/138‑1 (MQ=255)

            aaaGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTc                                                                                                                                       <  1:575073/150‑1 (MQ=255)

            aaaGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTc                                                                                                                                       >  2:8132/1‑150 (MQ=255)

              aGCAATCCGTTTGCTCTTTACCCCCGGTTGGGGGGTCGTCCACAAGCTTCGAGAGTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCCCGCAACGACTTTTTGTGGAGCTTAGGCTCCCCTcgtcgt                                                                                                                                     >  1:1084697/1‑150 (MQ=255)

                  aTCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACACAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCt                                                                                                                                 >  2:161379/1‑150 (MQ=255)

                     cGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTAcc                                                                                                                              >  2:906476/1‑150 (MQ=255)

                       tttGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGAGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCt                                                                                                                            >  1:2169945/1‑150 (MQ=255)

                          gCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCt                                                                                                                         <  2:760114/150‑1 (MQ=255)

                          gCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCt                                                                                                                         <  1:1514948/150‑1 (MQ=255)

                           ctctTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTg                                                                                                                        <  1:231444/150‑1 (MQ=255)

                           ctctTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTg                                                                                                                        <  1:509615/150‑1 (MQ=255)

                            tctTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTgg                                                                                                                       >  2:859114/1‑150 (MQ=255)

                             ctTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCCTAAGCTccccaag                                                                                                                                         >  2:326146/1‑128 (MQ=255)

                             ctTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCCTAAGCTccccaag                                                                                                                                         <  1:326146/131‑4 (MQ=255)

                               ttACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTggcag                                                                                                                    <  2:2255890/150‑1 (MQ=255)

                                   ccccGGTTTGGGGGTCGTCCGCAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGAAATTTTCTACTTCTTGGGGAACTTAGGCTCCCCTTGTCGGTTCTACCCTTCTGGCAGGCAt                                                                                                                >  1:1207098/1‑150 (MQ=255)

                                   ccccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACGGGGCACTCAACAAAAACAAGAAGCCCAGCAAAGAACCAATAAGAGCACGCAACGACTTCTTGGGGAGCTCAGGCTCCCCTTGTCGTTTCTACCCTTCTCGCCGGCAt                                                                                                                >  2:43859/1‑150 (MQ=255)

                                   ccccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCAt                                                                                                                >  1:314002/1‑150 (MQ=255)

                                   ccccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCAAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCAt                                                                                                                >  2:1371105/1‑150 (MQ=255)

                                   ccccGGTTTGGGGGTCGTCCACAAGATTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCAt                                                                                                                >  1:2157820/1‑150 (MQ=255)

                                    cccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATc                                                                                                               >  2:1526495/1‑150 (MQ=255)

                                     ccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATcc                                                                                                              >  2:393863/1‑150 (MQ=255)

                                     ccGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAACAGGTGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATcc                                                                                                              >  1:1333503/1‑150 (MQ=255)

                                      cGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCt                                                                                                             <  1:2188857/150‑1 (MQ=255)

                                         tttGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCt                                                                                                          <  1:321202/150‑1 (MQ=255)

                                               ggTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGc                                                                                                    <  1:906476/150‑1 (MQ=255)

                                               ggTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGc                                                                                                    <  2:1907184/150‑1 (MQ=255)

                                                gtcgtcCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGcc                                                                                                   <  2:130769/150‑1 (MQ=255)

                                                    tcCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTaca                                                                                               <  1:8132/150‑1 (MQ=255)

                                                    tcCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTaca                                                                                               <  1:859114/150‑1 (MQ=255)

                                                    tcCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTaca                                                                                               <  1:1940098/150‑1 (MQ=255)

                                                      cacaAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACAcc                                                                                             >  2:736861/1‑150 (MQ=255)

                                                      cacaAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTGCGGGGGGTCATCCTGCTACCCCCCTACAcc                                                                                             >  2:300608/1‑150 (MQ=255)

                                                       acaAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCACCCTGCTACCCGCCTACACCa                                                                                            >  1:2230005/1‑150 (MQ=255)

                                                           gCTTCGAGATTTGTGGACAGGGCGCTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTcc                                                                                        >  2:1185812/1‑150 (MQ=255)

                                                              tCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCg                                                                                     >  1:2003759/1‑150 (MQ=255)

                                                               cGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACACGCCTACACCATTCCCCGa                                                                                    <  1:393863/150‑1 (MQ=255)

                                                                       gtgGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCCTAAGCTccccaa                                                                                                                                          <  2:1435609/88‑3 (MQ=255)

                                                                               ggCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATcc                                                                    >  2:2338970/1‑150 (MQ=255)

                                                                               ggCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATcc                                                                    <  2:1333503/150‑1 (MQ=255)

                                                                               agCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATcc                                                                    <  2:2169945/149‑1 (MQ=255)

                                                                                  aCTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATcc                                                                                                              <  1:1443222/105‑1 (MQ=255)

                                                                                  aCTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATcc                                                                                                              >  2:1443222/1‑105 (MQ=255)

                                                                                    tCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGcc                                                               >  2:1777583/1‑150 (MQ=255)

                                                                                      aacaaAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTaca                                                                                               >  1:605649/1‑116 (MQ=255)

                                                                                      aacaaAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTaca                                                                                               <  2:605649/116‑1 (MQ=255)

                                                                                      aacaaAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCAAGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCca                                                             >  1:1006752/1‑150 (MQ=255)

                                                                                        caaAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCcatc                                                           <  2:314002/150‑1 (MQ=255)

                                                                                        caaAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGCGAGCTTAGGCTCCGCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCcatc                                                           <  2:1294930/150‑1 (MQ=255)

                                                                                         aaaaaCAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCcatca                                                          >  2:1010668/1‑150 (MQ=255)

                                                                                               aagaagCCCAGCAAAAAAACAATAAGAGCACGCAACGCCTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGc                                                    >  1:1496629/1‑150 (MQ=255)

                                                                                                agaagCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGcc                                                   >  1:1061626/1‑150 (MQ=255)

                                                                                                 gaagCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGcc                                                                                                   <  2:1156495/101‑1 (MQ=255)

                                                                                                 gaagCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGcc                                                                                                   >  1:1156495/1‑101 (MQ=255)

                                                                                                      ccAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATtctc                                             <  2:2254652/150‑1 (MQ=255)

                                                                                                       cAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCATGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTcc                                            <  1:161379/150‑1 (MQ=255)

                                                                                                        aGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCt                                           >  1:1221606/1‑150 (MQ=255)

                                                                                                        aGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCt                                           >  2:493317/1‑150 (MQ=255)

                                                                                                           aaaaaaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATtctc                                             >  2:909579/1‑145 (MQ=255)

                                                                                                           aaaaaaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATtctc                                             <  1:909579/145‑1 (MQ=255)

                                                                                                             aaaaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTc                                                                                                                          <  2:2395431/66‑1 (MQ=255)

                                                                                                             aaaaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTc                                                                                                                          >  1:2395431/1‑66 (MQ=255)

                                                                                                             aaaaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCt                                      >  2:2370837/1‑150 (MQ=255)

                                                                                                               aaaCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTAtt                                        >  2:2047536/1‑146 (MQ=255)

                                                                                                                 aCAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGc                                  <  2:601794/150‑1 (MQ=255)

                                                                                                                   aaTAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCg                                >  2:845281/1‑150 (MQ=255)

                                                                                                                    aTAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGa                               >  2:1028904/1‑150 (MQ=255)

                                                                                                                      atgtgcgcGCCACGACTTCTTGGGGAGCTTCGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAAc                             <  2:1084697/143‑1 (MQ=255)

                                                                                                                        gagCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACAt                           >  1:2445456/1‑150 (MQ=255)

                                                                                                                           cACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTCCACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTcc                        >  1:1779088/1‑150 (MQ=255)

                                                                                                                             cGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCtt                      <  2:155706/150‑1 (MQ=255)

                                                                                                                                aaCGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCa                                                  >  2:2118645/1‑119 (MQ=255)

                                                                                                                                aaCGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCa                                                  <  1:2118645/119‑1 (MQ=255)

                                                                                                                                     cttcttGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt              <  2:1815221/150‑1 (MQ=255)

                                                                                                                                     cttcttGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGt              <  1:999377/150‑1 (MQ=255)

                                                                                                                                      ttcttGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTg             >  1:1772254/1‑150 (MQ=255)

                                                                                                                                      ttcttGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTg             >  1:1283597/1‑150 (MQ=255)

                                                                                                                                       acaaggggAGCCTAAGCTCCCCTTGTCGTTTATAAACTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGc            >  1:1774041/5‑150 (MQ=255)

                                                                                                                                         ttGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCAt          <  2:1221606/150‑1 (MQ=255)

                                                                                                                                          tGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATc         <  2:1605019/150‑1 (MQ=255)

                                                                                                                                          tGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATc         <  2:1571316/150‑1 (MQ=255)

                                                                                                                                           ggggAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATcc        <  2:71267/150‑1 (MQ=255)

                                                                                                                                            gggAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATccc       >  2:497079/1‑150 (MQ=255)

                                                                                                                                             ggAGCCTAAGCTCCCCTTGTCGTTTCTACCCTTCCGGCATGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCa      <  2:1774041/150‑1 (MQ=255)

                                                                                                                                               aGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATg    <  1:50497/150‑1 (MQ=255)

                                                                                                                                                gCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATGc   <  1:2338970/150‑1 (MQ=255)

                                                                                                                                                gCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATGc   <  2:2255785/150‑1 (MQ=255)

                                                                                                                                                gCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATGc   <  2:1779088/150‑1 (MQ=255)

                                                                                                                                                 cTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATGCt  >  1:86237/1‑150 (MQ=255)

ACATCAGCGGCCAAAGCAATCCGTTTGCTCTTTACCCCCGGTTTGGGGGTCGTCCACAAGCTTCGAGATTTGTGGACAGGGCACTCAACAAAAACAAGAAGCCCAGCAAAAAAACAATAAGAGCACGCAACGACTTCTTGGGGAGCTTAGGCTCCCCTTGTCGTTTCTACCCTTCTGGCAGGCATCCTGCTACCCGCCTACACCATTCCCCGAGTAAATGCTAGAATCCCCGCCCATCATGCGGCCATTCTCCTATTCTTGGCCGAACATTCCTTCACACAGTGCATCCCATGCT  >  NC_002947/5338531‑5338825

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