BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,109,062	A→T	5.2%	I104N (ATT→AAT)	PP_1881 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,109,062	0	A	T	5.2%	354.8 / 4.6	114	I104N (ATT→AAT)	PP_1881	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (55/52); new base T (3/3); total (59/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 118 total aligned reads displayed.

TGTGAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGTGC  >  NC_002947/2108915‑2109209

tgtgAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATg                                                                                                                                                   >  1:610092/1‑150 (MQ=255)

  tgAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGc                                                                                                                                                 >  2:633725/1‑150 (MQ=255)

   gAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTGGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATTGcc                                                                                                                                                <  2:1440615/150‑1 (MQ=255)

    aaaCATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGccc                                                                                                                                               >  1:367647/1‑150 (MQ=255)

    aaaCATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGccc                                                                                                                                               >  2:619365/1‑150 (MQ=255)

    aaaCATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGccc                                                                                                                                               >  2:1855659/1‑150 (MQ=255)

       cATATTTGTCGCCCTTGTTGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGCCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCATTGGCCCTgt                                                                                                                                            <  1:1043596/150‑2 (MQ=255)

        atatTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGaa                                                                                                                                           >  1:1101454/1‑150 (MQ=255)

        atatTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGaa                                                                                                                                           <  2:1636500/150‑1 (MQ=255)

         tatTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAg                                                                                                                                          <  2:829858/150‑1 (MQ=255)

         tatTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAg                                                                                                                                          <  1:1161360/150‑1 (MQ=255)

              gTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGtt                                                                                                                                     >  2:223846/1‑150 (MQ=255)

                  cccTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGgtgt                                                                                                                                      <  2:480731/145‑1 (MQ=255)

                  cccTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGgtgt                                                                                                                                      >  1:480731/1‑145 (MQ=255)

                        tgGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCa                                                                                                                           <  1:2296026/150‑1 (MQ=255)

                                cAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGACCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAggg                                                                                                                   >  2:948429/1‑149 (MQ=255)

                                  gTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAg                                                                                                                     <  2:1712493/146‑1 (MQ=255)

                                  gTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAg                                                                                                                     >  1:1712493/1‑146 (MQ=255)

                                  gTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAAATTTTCGACGGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGa                                                                                                                 <  2:327989/150‑1 (MQ=255)

                                   ttCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGaa                                                                                                                >  2:1387992/1‑150 (MQ=255)

                                       tAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGCAGGTGTTTTCGGCGCCAGTCAAGGt                                                                                                                   >  2:295221/1‑143 (MQ=255)

                                       tAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGt                                                                                                                   <  1:295221/143‑1 (MQ=255)

                                          cGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGt                                                                                                         <  2:625213/150‑1 (MQ=255)

                                           gAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCTCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGtt                                                                                                        <  1:2223403/150‑1 (MQ=255)

                                           aaaTCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGtt                                                                                                        <  1:1959733/149‑1 (MQ=255)

                                              tCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTa                                                                                                     >  1:132918/1‑150 (MQ=255)

                                                   aGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTc                                                                                                <  1:782028/150‑1 (MQ=255)

                                                     tttGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCaa                                                                                                                                                     >  2:1321496/1‑95 (MQ=255)

                                                     tttGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCaa                                                                                                                                                     <  1:1321496/95‑1 (MQ=255)

                                                          cATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGagtagt                                                                                         >  2:106052/1‑150 (MQ=255)

                                                          cATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGagtagt                                                                                         >  2:2105235/1‑150 (MQ=255)

                                                           aTTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTc                                                                                        >  1:1523316/1‑150 (MQ=255)

                                                            ttGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTc                                                                                                <  2:1013985/141‑1 (MQ=255)

                                                            ttGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTc                                                                                                >  1:1013985/1‑141 (MQ=255)

                                                              gCGTGGAATGGACGCCTGGGGAATCTGAAGTTCTTACTTAAAGTCTCCTGGTAAGAGGGGGGGGGGCAAGACGGCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCAtt                                                                                     <  1:633725/150‑1 (MQ=255)

                                                               ggTGGAATGGACGCCTGGGGAAAAGCGAGTGCTAACGGAAAGTCTCCTGGCAAGAGGGGCGGGAGAAACTTTTCGACAGGGACAATGGCCCTGGAGGTGTTTTCCCCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCAttt                                                                                    <  1:996667/149‑1 (MQ=255)

                                                               cGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGgt                                                                                                                                        <  2:134194/98‑1 (MQ=255)

                                                                gTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTg                                                                                   >  2:1326068/1‑150 (MQ=255)

                                                                gTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTg                                                                                   >  1:1743416/1‑150 (MQ=255)

                                                                 tGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCAt                                                                                      >  2:10632/1‑146 (MQ=255)

                                                                        gACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTg                                                                           >  2:1159273/1‑150 (MQ=255)

                                                                           gCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAg                                                                                                                     <  1:1102381/105‑1 (MQ=255)

                                                                           gCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGccc                                                                        >  2:1235666/1‑150 (MQ=255)

                                                                           gCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGccc                                                                        >  1:693470/1‑150 (MQ=255)

                                                                           gCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCAGGCCCGACTGAGCGCAGAACTGCTGAGACAGGACCAGGGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGccc                                                                        >  1:635858/1‑150 (MQ=255)

                                                                            ccTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCg                                                                       <  1:106052/150‑1 (MQ=255)

                                                                             cTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAg                                                                                                                     <  2:1021920/103‑1 (MQ=255)

                                                                             cTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAg                                                                                                                     >  1:1021920/1‑103 (MQ=255)

                                                                             cTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCgg                                                                      >  2:2116891/1‑150 (MQ=255)

                                                                              tGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGt                                                                     <  1:2105235/150‑1 (MQ=255)

                                                                              tGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGCTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGt                                                                     >  1:2236346/1‑150 (MQ=255)

                                                                               ggggAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTaa                                                                                                    <  1:715841/118‑1 (MQ=255)

                                                                               ggggAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTg                                                                    <  1:1855659/150‑1 (MQ=255)

                                                                               ggggAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTg                                                                    <  2:1101454/150‑1 (MQ=255)

                                                                                gggAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTgg                                                                   <  2:610092/150‑1 (MQ=255)

                                                                                 ggAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAGCTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGag                                                                                             >  2:396508/1‑123 (MQ=255)

                                                                                 ggAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAGCTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGag                                                                                             <  1:396508/123‑1 (MQ=255)

                                                                                 ggAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGc                                                                  >  2:770539/1‑150 (MQ=255)

                                                                                  gAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCa                                                                 >  2:1304708/1‑150 (MQ=255)

                                                                                   aaaaGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTGTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAGCTGGAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAg                                                                <  1:1590033/150‑1 (MQ=255)

                                                                                     aaGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGt                                                              <  1:1304708/150‑1 (MQ=255)

                                                                                       gCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACGTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCa                                                            <  1:765276/150‑1 (MQ=255)

                                                                                            ggcgaaCGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCa                                                       >  2:34127/5‑150 (MQ=255)

                                                                                                aaCGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCGCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCg                                                   >  2:1510130/1‑150 (MQ=255)

                                                                                                 aCGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCgg                                                  <  1:1387992/150‑1 (MQ=255)

                                                                                                   ggAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCgggg                                                >  2:2219397/1‑150 (MQ=255)

                                                                                                    gAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGa                                               >  2:1221117/1‑150 (MQ=255)

                                                                                                     aaaGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGAt                                              >  1:291045/1‑150 (MQ=255)

                                                                                                              tGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGAtc                                             <  2:1632804/142‑1 (MQ=255)

                                                                                                                gCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCgg                                                  >  1:1636032/1‑135 (MQ=255)

                                                                                                                    gAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGAt                                                                                >  1:1959659/1‑101 (MQ=255)

                                                                                                                       gggCGGGAGAAACTTTTCGACAGGGACAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGt                            <  2:1003730/150‑1 (MQ=255)

                                                                                                                        ggCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGtt                           >  1:1569701/1‑150 (MQ=255)

                                                                                                                            ggAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACTGTCAGGTCGGTTTTCa                       <  1:1221117/150‑1 (MQ=255)

                                                                                                                            ggAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCa                       >  2:930794/1‑150 (MQ=255)

                                                                                                                            ggAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCa                       >  2:2164520/1‑150 (MQ=255)

                                                                                                                             gagaAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAg                      <  1:1159273/150‑1 (MQ=255)

                                                                                                                              agaACCTTTTCGACACGGCCAAGGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTCCAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAgg                     <  1:1096015/150‑1 (MQ=255)

                                                                                                                                  aCTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGTATCTCGACGGTCAGGTCGGTTTTCAGGCGGt                 <  1:1510130/150‑1 (MQ=255)

                                                                                                                                   cTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTc                >  2:2132948/1‑150 (MQ=255)

                                                                                                                                   cTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTc                >  2:1468667/1‑150 (MQ=255)

                                                                                                                                   cTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTc                >  2:1468680/1‑150 (MQ=255)

                                                                                                                                    ttttCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCgg                             >  2:1206122/1‑136 (MQ=255)

                                                                                                                                    ttttCGAAAGTGCTAAGGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCt               <  2:132918/150‑1 (MQ=255)

                                                                                                                                      ttCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCttt             <  1:1007288/150‑1 (MQ=255)

                                                                                                                                        cttcAGGGCCATTGGCCCTGTAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTtgt           >  2:849712/4‑150 (MQ=255)

                                                                                                                                        cttcAGGGCCATTGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGc                                                                  <  2:1204419/92‑1 (MQ=255)

                                                                                                                                        cttcAGGGCCATTGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGc                                                                  >  1:1204419/4‑95 (MQ=255)

                                                                                                                                         ttcAGGGCCATTGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTtgtg          >  1:1307593/3‑150 (MQ=255)

                                                                                                                                         gACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTtttg          <  2:1743416/150‑1 (MQ=255)

                                                                                                                                             gggCCATTGGCCCTGTAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGc      <  1:849712/150‑1 (MQ=255)

                                                                                                                                              ggCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGAt                                              <  2:1618327/109‑1 (MQ=255)

                                                                                                                                              ggCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGAt                                              >  1:1618327/1‑109 (MQ=255)

                                                                                                                                              ggCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCg     >  2:2326165/1‑150 (MQ=255)

                                                                                                                                               gCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGAc                                                         <  1:1678871/97‑1 (MQ=255)

                                                                                                                                               gCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGAc                                                         >  2:1678871/1‑97 (MQ=255)

                                                                                                                                                 cAATGGCCTTGACGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGTGc  <  1:948429/150‑1 (MQ=255)

                                                                                                                                                 cAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATTTGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCg                   <  1:1033472/133‑1 (MQ=255)

                                                                                                                                                 cAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCg                   >  2:1033472/1‑133 (MQ=255)

                                                                                                                                                   aTGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTtg            <  2:2338154/138‑1 (MQ=255)

TGTGAAACATATTGGTCGCCCTTGTGGTCGACCAGTTCGTAGCGAATCCACAGTTTGCCATTGCGTGGAATGGACGCCTGGGGAAAAGCCAGTGCTAACGGAAAGTCTCCTGGCACGAGGGGCGGGAGAAACTTTTCGACAGGGCCAATGGCCCTGAAGGTGTTTTCGGCGCCAGTCAAGGTGAACTGCAGTTGTAACGTCGAGTAGTCATTTGGATCACTGCCCGGTGGCAGGTCAGACCATGCGGGGATCTCGACGGTCAGGTCGGTTTTCAGGCGGTCTTTGTGCAGCGTGC  >  NC_002947/2108915‑2109209

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