BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,781,636	A→T	5.0%	L828Q (CTG→CAG)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,781,636	0	A	T	5.0%	332.8 / 6.0	99	L828Q (CTG→CAG)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand): ref base A (41/53); new base T (2/3); total (43/56)
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.73e-01

Only 100 of 103 total aligned reads displayed.

GGCGGCTGCACAGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGATGATGGTCGGGTAGTCCTCGGG  >  NC_002947/1781490‑1781781

ggcggcTGCACAGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCcagca                                                                                                                                                <  2:568413/150‑1 (MQ=255)

 gcggcTGCACAGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCcagcag                                                                                                                                               <  1:1809387/150‑1 (MQ=255)

   ggcTGCACAGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGc                                                                                                                                             <  2:2047382/150‑1 (MQ=255)

    gcTGCACAGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGcc                                                                                                                                            >  2:2398372/1‑150 (MQ=255)

      tGCACAGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCg                                                                                                                                          >  1:954677/1‑150 (MQ=255)

      tGCACAGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCg                                                                                                                                          >  2:1896113/1‑150 (MQ=255)

       gCACAGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCgg                                                                                                                                         <  1:1189746/150‑1 (MQ=255)

         acaGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCTGCTGGCCCggct                                                                                                                                       <  2:1098587/150‑3 (MQ=255)

         acaGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGc                                                                                                                                       >  2:221729/1‑150 (MQ=255)

             cTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTc                                                                                                                                   <  1:368179/150‑1 (MQ=255)

                      tGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGt                                                                                                                          <  2:1036563/150‑1 (MQ=255)

                        aCAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGa                                                                                                                        <  1:1755763/150‑1 (MQ=255)

                           gCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGAttt                                                                                                                     <  2:715624/150‑1 (MQ=255)

                              gCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGa                                                                                                                  <  1:1075658/150‑1 (MQ=255)

                              gCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGAAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGa                                                                                                                  <  1:162802/150‑1 (MQ=255)

                                   ttGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATg                                                                                                             <  1:2388510/150‑1 (MQ=255)

                                   ttGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATg                                                                                                             <  1:1426295/150‑1 (MQ=255)

                                     gTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGt                                                                                                           <  1:1248315/150‑1 (MQ=255)

                                     gTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGt                                                                                                           >  2:1248315/1‑150 (MQ=255)

                                         gCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCcagcag                                                                                                                                               >  1:2183462/1‑110 (MQ=255)

                                         gCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCcagcag                                                                                                                                               <  2:2183462/110‑1 (MQ=255)

                                           aTCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCTGCTGGCCCggctg                                                                                                                                      <  2:554051/117‑4 (MQ=255)

                                           aTCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCTGCTGGCCCggctg                                                                                                                                      >  1:554051/1‑114 (MQ=255)

                                              ggTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCg                                                                                                  >  2:2259525/1‑150 (MQ=255)

                                              ggTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCg                                                                                                  >  1:2310163/1‑150 (MQ=255)

                                               gTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGc                                                                                                 >  1:2346407/1‑150 (MQ=255)

                                                 aaTGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGGCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTg                                                                                               <  2:1452135/150‑1 (MQ=255)

                                                  aTGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTgg                                                                                              <  2:1573108/150‑1 (MQ=255)

                                                    gaagaaCACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGt                                                                                                                          >  2:1007391/1‑120 (MQ=255)

                                                    gaagaaCACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGt                                                                                                                          <  1:1007391/120‑1 (MQ=255)

                                                       gaaCACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGa                                                                                                                  <  2:447426/125‑1 (MQ=255)

                                                       gaaCACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGa                                                                                                                  >  1:447426/1‑125 (MQ=255)

                                                           acGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTgg                                                                                              >  2:1109099/1‑141 (MQ=255)

                                                           acGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTgg                                                                                              <  1:1109099/141‑1 (MQ=255)

                                                            cGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGa                                                                                    <  1:2398372/150‑1 (MQ=255)

                                                            cGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGa                                                                                    >  1:1780842/1‑150 (MQ=255)

                                                             gTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGt                                                                                                                          <  1:1960149/111‑1 (MQ=255)

                                                             gTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGt                                                                                                                          >  2:1960149/1‑111 (MQ=255)

                                                               cttctACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCAGGGCCAGCAGGCCCGGGCGGTGTAACTCGGTGATTTCGATGATGGTCACCGGCCGCTGGGCATCGTTGAGGa                                                                                 <  1:1622894/150‑1 (MQ=255)

                                                                ttctACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGAt                                                                                >  1:1707517/1‑150 (MQ=255)

                                                                ttctACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGAt                                                                                >  2:313913/1‑150 (MQ=255)

                                                                 tctACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATg                                                                               >  1:1228079/1‑150 (MQ=255)

                                                                    aCCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGTAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTc                                                                            <  1:1084038/150‑1 (MQ=255)

                                                                     cccGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCa                                                                           >  2:2380938/1‑150 (MQ=255)

                                                                        gCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACtt                                                                        <  1:2051724/150‑1 (MQ=255)

                                                                        gCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACtt                                                                        >  2:2056844/1‑150 (MQ=255)

                                                                             ccGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTgcggc                                                                   >  2:1376402/1‑150 (MQ=255)

                                                                             ccGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTgcggc                                                                   <  2:2310163/150‑1 (MQ=255)

                                                                               gAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGcggcgg                                                                 <  1:221729/150‑1 (MQ=255)

                                                                               gAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGcggcgg                                                                 <  2:2311595/150‑1 (MQ=255)

                                                                                  ggTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGaa                                                              <  2:1707517/150‑1 (MQ=255)

                                                                                  ggTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGaa                                                              <  2:2346407/150‑1 (MQ=255)

                                                                                   gTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAg                                                             >  2:2074228/1‑150 (MQ=255)

                                                                                        gATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGaa                                                        <  2:669760/150‑1 (MQ=255)

                                                                                        gATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGaa                                                        <  1:2380938/150‑1 (MQ=255)

                                                                                            ttGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGTGGGAAGTCGAAGTGc                                                    <  1:1896113/150‑1 (MQ=255)

                                                                                            ttGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGc                                                    <  1:2074228/150‑1 (MQ=255)

                                                                                            ttGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGc                                                    <  1:1669999/150‑1 (MQ=255)

                                                                                            ttGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGc                                                    >  1:1969497/1‑150 (MQ=255)

                                                                                             tGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCt                                                   <  2:800175/150‑1 (MQ=255)

                                                                                              ggCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCtt                                                  >  2:54202/1‑150 (MQ=255)

                                                                                               gCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTg                                                 >  1:108627/1‑150 (MQ=255)

                                                                                                   tttGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCTAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCt                                             >  1:1951319/1‑150 (MQ=255)

                                                                                                        aGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCcgc                                        >  1:2037901/1‑150 (MQ=255)

                                                                                                          tGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGc                                                    >  1:1272245/1‑136 (MQ=255)

                                                                                                          tGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGc                                                    <  2:1272245/136‑1 (MQ=255)

                                                                                                             cAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTAc                                   <  2:182576/150‑1 (MQ=255)

                                                                                                                gTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCt                                                   <  2:1800465/131‑1 (MQ=255)

                                                                                                                gTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCt                                                   >  1:1800465/1‑131 (MQ=255)

                                                                                                                 tCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCgg                               >  1:908522/1‑150 (MQ=255)

                                                                                                                   gAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGgcg                             <  2:1228079/150‑1 (MQ=255)

                                                                                                                   gAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGgcg                             <  2:1184261/150‑1 (MQ=255)

                                                                                                                       tCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCtt                                                  >  1:651829/1‑125 (MQ=255)

                                                                                                                       tCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCtt                                                  <  2:651829/125‑1 (MQ=255)

                                                                                                                       tCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTgg                         >  1:1820622/1‑150 (MQ=255)

                                                                                                                       tCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTgg                         >  1:233388/1‑150 (MQ=255)

                                                                                                                       tCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTgg                         >  1:1958778/1‑150 (MQ=255)

                                                                                                                       tCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGCAATTTCGGGGATGGTCACCTTTCGCTGGTCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTgg                         >  1:295788/1‑150 (MQ=255)

                                                                                                                        ccAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGga                        <  2:108627/150‑1 (MQ=255)

                                                                                                                        ccAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGga                        <  1:313913/150‑1 (MQ=255)

                                                                                                                         cAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGgat                       >  2:483692/1‑150 (MQ=255)

                                                                                                                         cAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGgat                       >  2:2297786/1‑150 (MQ=255)

                                                                                                                            gAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCtt                                                  >  2:1713873/1‑120 (MQ=255)

                                                                                                                            gAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCtt                                                  <  1:1713873/120‑1 (MQ=255)

                                                                                                                            gAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGgatgat                    >  2:2261137/1‑150 (MQ=255)

                                                                                                                             aaGATCCGCCCCAGCCGGGCGAGCAGGCCCGGGCGGTCTCGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGatgatg                   <  2:219996/150‑1 (MQ=255)

                                                                                                                             aaGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCGGGTTCCGCTGGGCATCGATGAGGATGGTCACTTGCGGCGGGCAGTCGAAGTCCGTGAGCCGCCGCGGTACCCGGCGCTGGatgatg                   <  1:1376402/150‑1 (MQ=255)

                                                                                                                             aaGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGatgatg                   <  2:2037901/150‑1 (MQ=255)

                                                                                                                                 tCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGaa                                                        >  2:615810/1‑109 (MQ=255)

                                                                                                                                 tCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGaa                                                        <  1:615810/109‑1 (MQ=255)

                                                                                                                                  ccGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGGGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCCGCGCGGGATGTTGGTctg              <  2:531966/150‑3 (MQ=255)

                                                                                                                                  ccGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGATGATGGTCgg              <  2:1969497/150‑1 (MQ=255)

                                                                                                                                   cGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGATGATGGTCggg             <  1:698233/150‑1 (MQ=255)

                                                                                                                                    gCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGATGATGGTCGGGt            >  1:1842729/1‑150 (MQ=255)

                                                                                                                                     ccccAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGATGATGGTCGGGTa           >  1:1063517/1‑150 (MQ=255)

                                                                                                                                      cccAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGaa                                                              <  2:2285884/98‑1 (MQ=255)

                                                                                                                                       ccAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGATGATGGTCGGGTAGt         <  2:1958778/150‑1 (MQ=255)

                                                                                                                                        cgcccGGGCCTGCTGGCCC‑GGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTgcggcg                                                                  >  1:803858/4‑91 (MQ=255)

                                                                                                                                        cgcccGGGCCTGCTGGCCC‑GGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTgcggcg                                                                  <  2:803858/88‑1 (MQ=255)

                                                                                                                                              ggCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGTTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGTGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGATGATGGTCGGGTAGTCCTCggg  <  1:1828691/150‑1 (MQ=255)

GGCGGCTGCACAGCTGCGGGTCTGACAGCGGCTGGTTGTCGGCATCGGTAATGAAGAACACGTCTTCTACCCGCTCGCCGAGGGTGGCGATCTTGGCGTTTTGCAGTGACAGGTCGAACTCCAGGAAGATCCGCCCCAGCCGGGCCAGCAGGCCCGGGCGGTCTGGCGCGGTGATTTCGAGGATGGTCACCGGCCGCTGGGCATCGTTGAGGATGGTCACTTGCGGCGGGAAGTCGAAGTGCTTGAGCTGCCGCGGTACCCGGCGCTGGATGATGGTCGGGTAGTCCTCGGG  >  NC_002947/1781490‑1781781

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