BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,857,969	A→C	7.6%	T4T (ACT→ACG) ‡	PP_4272 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,857,969	0	A	C	7.6%	285.9 / 9.5	100	T4T (ACT→ACG) ‡	PP_4272	membrane protein
Reads supporting (aligned to +/- strand): ref base A (52/40); new base C (3/5); total (55/45)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

TCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGA  >  NC_002947/4857827‑4858116

tCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAGGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACt                                                                                                                                              >  2:1615856/1‑150 (MQ=255)

 cagcagCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCACTAAACGTCATTTGAACAAGCCCTCGAAGGTTAACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACtt                                                                                                                                             >  1:773020/1‑150 (MQ=255)

   gcagcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGACTTCAACTtct                                                                                                                                           <  2:10591/150‑2 (MQ=255)

   gcagcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGACTTCAACTtct                                                                                                                                           <  2:10179/150‑2 (MQ=255)

    cagcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTtcaa                                                                                                                                          >  1:1524153/1‑150 (MQ=255)

    cagcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACCGGCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTtcaa                                                                                                                                          >  1:519171/1‑150 (MQ=255)

     agcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGACTTCAACTtctat                                                                                                                                         <  2:1626515/150‑4 (MQ=255)

     agcaACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTtcaa                                                                                                                                          <  2:1953154/149‑1 (MQ=255)

       caACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTtcaagta                                                                                                                                       >  2:816701/1‑149 (MQ=255)

         aCCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTAAACAAGCCCTCGAAGGTCTACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTtcaagtcaa                                                                                                                                     <  1:1978515/150‑1 (MQ=255)

           ccccACAGGTCCATCAATACAGTGGCACCAATGCAAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAgg                                                                                                                                   >  2:1481243/1‑150 (MQ=255)

               acaGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAg                                                                                                                                                    <  1:1632722/129‑1 (MQ=255)

               acaGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAg                                                                                                                                                    >  2:1632722/1‑129 (MQ=255)

                  ggTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGaa                                                                                                                            >  1:477868/1‑150 (MQ=255)

                   gTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGaaa                                                                                                                           <  1:2363898/150‑1 (MQ=255)

                    tCCATAACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAAt                                                                                                                          <  1:103271/150‑1 (MQ=255)

                        tCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGAc                                                                                                                      >  1:1958540/1‑150 (MQ=255)

                         cACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCACCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACa                                                                                                                     <  1:2359778/150‑1 (MQ=255)

                            tACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTg                                                                                                                  >  1:1937673/1‑150 (MQ=255)

                              ccgTGGCACGAATGCCAATCCTCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGTTTGACGGTCATTACCATGCCTCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGcc                                                                                                                <  2:519171/148‑1 (MQ=255)

                                gTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGa                                                                                                              <  1:1615856/150‑1 (MQ=255)

                                 tGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGaa                                                                                                                                                     <  2:1068118/110‑1 (MQ=255)

                                 tGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGaa                                                                                                                                                     >  1:1068118/1‑110 (MQ=255)

                                   gCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGt                                                                                                           >  1:1187933/1‑150 (MQ=255)

                                     aCCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCg                                                                                                         >  2:1369113/1‑150 (MQ=255)

                                      ccAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGc                                                                                                        <  2:477868/150‑1 (MQ=255)

                                              aaTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTtcaagtcaag                                                                                                                                    <  1:1714950/114‑1 (MQ=255)

                                              aaTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTtcaagtcaag                                                                                                                                    >  2:1714950/1‑114 (MQ=255)

                                              aaTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCa                                                                                                >  2:250226/1‑150 (MQ=255)

                                               aTCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCAc                                                                                               >  2:1370411/1‑150 (MQ=255)

                                                   gCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGaa                                                                                                                                                     >  1:1660921/1‑92 (MQ=255)

                                                     aaTGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGt                                                                                         >  1:1525410/1‑150 (MQ=255)

                                                     aaTGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGt                                                                                         >  2:1683221/1‑150 (MQ=255)

                                                        gCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTAcc                                                                                      <  2:1187933/150‑1 (MQ=255)

                                                         cAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCt                                                                                     <  1:719830/150‑1 (MQ=255)

                                                            cTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGGCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCa                                                                                  >  2:654444/1‑150 (MQ=255)

                                                              aGAAACATCGCAGTAAACGTCATTCGAACAAGCCCTCGAAGGTTGACGGTCACTCCCCCG‑ATCTCCCCACGCGTTATTGGCGTCAACTGCAAGTCACGGGCCCAACCTGGACCTTGCCGACGTCGCCAAACTCACTGGCGTACCTGCATc                                                                                <  1:551762/150‑1 (MQ=255)

                                                               gAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCg                                                                               >  2:1444211/1‑150 (MQ=255)

                                                                aaaCATCGTAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGa                                                                              <  2:2306357/150‑1 (MQ=255)

                                                                aaaCATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGa                                                                              >  1:1485370/1‑150 (MQ=255)

                                                                 aaCATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGac                                                                             >  1:1276158/1‑150 (MQ=255)

                                                                  aCATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCa                                                                                                <  1:709867/130‑1 (MQ=255)

                                                                  aCATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCa                                                                                                >  2:709867/1‑130 (MQ=255)

                                                                    aTCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCG‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAGACGAACTGGCGTACTTGCATCGACACt                                                                          >  2:1811490/1‑150 (MQ=255)

                                                                    aTCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGACAACTTCAAGTGAAGGGCCCGAAATGCACATTGCCGACGTCGCCAAACGCACTCGCATACCTCCATAAacacc                                                                          >  2:1209775/1‑149 (MQ=255)

                                                                      cGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGc                                                                        <  1:1369113/150‑1 (MQ=255)

                                                                         aGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGtt                                                                     <  1:87630/150‑1 (MQ=255)

                                                                            aaaCGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCg                                                                                                         >  1:639769/1‑111 (MQ=255)

                                                                            aaaCGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCg                                                                                                         <  2:639769/111‑1 (MQ=255)

                                                                              aCGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTtactac                                                                >  1:2330272/1‑150 (MQ=255)

                                                                                gTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACga                                                              >  1:1385488/1‑150 (MQ=255)

                                                                                   aTTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGagaa                                                           >  1:1422511/1‑150 (MQ=255)

                                                                                    tttGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGc                                                                                           >  2:251780/1‑117 (MQ=255)

                                                                                    tttGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGc                                                                                           <  1:251780/117‑1 (MQ=255)

                                                                                    tttGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGagaag                                                          >  1:594953/1‑150 (MQ=255)

                                                                                    tttGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGagaag                                                          <  1:2018424/150‑1 (MQ=255)

                                                                                    tttGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGagaag                                                          >  1:500982/1‑150 (MQ=255)

                                                                                      tGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAgaagaa                                                        <  2:1485370/150‑1 (MQ=255)

                                                                                       gAACAAGCCCTCGAAGGTTGACGGTCATTACAATGCCACGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCACGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGaaa                                                       <  2:2296338/150‑1 (MQ=255)

                                                                                        aaCAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAg                                                      >  2:1529146/1‑150 (MQ=255)

                                                                                        aaCAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAg                                                      >  1:683383/1‑150 (MQ=255)

                                                                                         aCAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAgg                                                     >  2:670003/1‑150 (MQ=255)

                                                                                         aCAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAgg                                                     <  1:1481243/150‑1 (MQ=255)

                                                                                         aCAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAg                                                      <  1:995407/149‑1 (MQ=255)

                                                                                         aCAAGCCCGCGAAGGTTGACGGGCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATCGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAgg                                                     <  2:1446677/150‑1 (MQ=255)

                                                                                           aaGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCt                                                   >  1:1181659/1‑150 (MQ=255)

                                                                                             gCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTg                                                 >  1:1808542/1‑150 (MQ=255)

                                                                                                cTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTc                                              >  2:1381669/1‑150 (MQ=255)

                                                                                                 tCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGCCGTCCCCAGGGAATCAATGACCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCa                                             >  2:1791808/1‑150 (MQ=255)

                                                                                                 tCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAAc                                                                                                   >  1:1913953/1‑96 (MQ=255)

                                                                                                 tCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAAc                                                                                                   <  2:1913953/96‑1 (MQ=255)

                                                                                                 tCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCa                                             >  2:1603547/1‑150 (MQ=255)

                                                                                                  cGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCaa                                            <  2:1276158/150‑1 (MQ=255)

                                                                                                        ttGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCtt                                      <  1:1683221/150‑1 (MQ=255)

                                                                                                        ttGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCtt                                      >  1:333728/1‑150 (MQ=255)

                                                                                                        ttGACGGTCATTACAATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGAAAAAAAAAAAAAAAAAAGGCTTGCTCAAGTCCtt                                      <  2:2407000/150‑1 (MQ=255)

                                                                                                         tGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTg                                     >  2:454966/1‑150 (MQ=255)

                                                                                                           aCGGTCGTTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGc                                   >  2:1199002/1‑150 (MQ=255)

                                                                                                              gTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCgg                                >  1:1387118/1‑150 (MQ=255)

                                                                                                               tCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCCGCATCGACACTTCGTTACTACGAGAATAAAGCCGTGCTCAAGTCCTTGGCCgcg                               >  2:1330887/1‑148 (MQ=255)

                                                                                                                 aTTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGa                                                                              >  1:2320575/1‑101 (MQ=255)

                                                                                                                 aTTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGa                                                                              <  2:2320575/101‑1 (MQ=255)

                                                                                                                    aCCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCgg                          <  1:250226/150‑1 (MQ=255)

                                                                                                                      cATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGagaaga                                                         >  2:1236988/1‑117 (MQ=255)

                                                                                                                      cATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGagaaga                                                         <  1:1236988/117‑1 (MQ=255)

                                                                                                                       aTGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCa                       <  2:500982/150‑1 (MQ=255)

                                                                                                                       aTGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCa                       <  2:1385488/150‑1 (MQ=255)

                                                                                                                              gcagca‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGa                <  2:1808542/150‑1 (MQ=255)

                                                                                                                              gcagca‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGa                <  1:670003/150‑1 (MQ=255)

                                                                                                                                    cttgAAGTTGACTTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGtt           >  2:1355475/4‑150 (MQ=255)

                                                                                                                                     ttgAAGTTGACTTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTgg                                    >  1:1201342/3‑124 (MQ=255)

                                                                                                                                     ttgAAGTTGACTTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTgg                                    <  2:1201342/122‑1 (MQ=255)

                                                                                                                                        aaGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCa       <  2:683383/150‑1 (MQ=255)

                                                                                                                                          gTTGAAGTCAACTTCAAGTCACGGGCACGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGc     >  1:1089653/1‑150 (MQ=255)

                                                                                                                                          gTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGc     <  1:816701/150‑1 (MQ=255)

                                                                                                                                           ttGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCg                                                                       >  2:1991512/1‑83 (MQ=255)

                                                                                                                                           ttGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCg                                                                       <  1:1991512/83‑1 (MQ=255)

                                                                                                                                           ttGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCtt                                      <  2:353805/116‑1 (MQ=255)

                                                                                                                                           ttGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCtt                                      >  1:353805/1‑116 (MQ=255)

                                                                                                                                             gAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGa  <  1:1199002/150‑1 (MQ=255)

TCAGCAGCAACCCCCACAGGTCCATCACTACAGTGGCACCAATGCCAATCGGCAATGCAGCTAGAAACATCGCAGTAAACGTCATTTGAACAAGCCCTCGAAGGTTGACGGTCATTACCATGCCGCGCAGCA‑TAGAAGTTGAAGTCAACTTCAAGTCAAGGGCCCGAAATGGACATTGCCGACGTCGCCAAACGCACTGGCGTACCTGCATCGACACTGCGTTACTACGAGAAGAAAGGCTTGCTCAAGTCCTTGGCCGGGCGCGGCCAGCGTCGACAGTTTGCAGCCGA  >  NC_002947/4857827‑4858116

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