BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,211,571	A→G	R44R (CGT→CGC)	yjaB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,211,571	0	A	G	99.7%	1113.5 / ‑6.1	359	R44R (CGT→CGC)	yjaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (183/175); total (183/176)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 423 total aligned reads displayed.

ATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTC  >  NC_000913/4211481‑4211660

aTCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCTGCAGGAAGgagcg                                                                                         >  2:6096882/1‑93 (MQ=255)

aTCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcg                                                                                         >  1:2715345/1‑93 (MQ=255)

aTCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcg                                                                                         >  1:1752164/1‑93 (MQ=255)

 tCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcga                                                                                        >  2:6425038/1‑93 (MQ=255)

 tCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcga                                                                                        >  2:6720274/1‑93 (MQ=255)

 tCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcga                                                                                        >  2:1291996/1‑93 (MQ=255)

 tCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcga                                                                                        <  2:1835717/93‑1 (MQ=255)

 accATATGCTGCCTACTAAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcga                                                                                        <  1:6324722/92‑1 (MQ=255)

   cATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcgaac                                                                                      >  1:7782791/1‑93 (MQ=255)

   cATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcgaac                                                                                      <  2:3326468/93‑1 (MQ=255)

   cATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcgaac                                                                                      <  2:901643/93‑1 (MQ=255)

   cATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCAACCCACAACGGCGCTTCCGGCAGGAAGgagcgaac                                                                                      <  1:2227055/93‑1 (MQ=255)

     tatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTa                                                                                    >  2:4203326/1‑93 (MQ=255)

     tatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTa                                                                                    >  1:3377990/1‑93 (MQ=255)

      atGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAg                                                                                   <  1:7292589/93‑1 (MQ=255)

       tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAgg                                                                                  <  2:1465682/93‑1 (MQ=255)

       tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAgg                                                                                  >  2:4534189/1‑93 (MQ=255)

       tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAgg                                                                                  <  2:489314/93‑1 (MQ=255)

       tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAgg                                                                                  >  2:5577617/1‑93 (MQ=255)

       tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAgg                                                                                  >  2:3044380/1‑93 (MQ=255)

        gctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGt                                                                                 <  2:6501893/93‑1 (MQ=255)

          tgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcc                                                                               >  1:6201195/1‑93 (MQ=255)

              cACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCa                                                                           >  1:4377740/1‑93 (MQ=255)

               aCTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAg                                                                          >  2:435118/1‑93 (MQ=255)

               aCTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAg                                                                          >  2:586955/1‑93 (MQ=255)

                cTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGc                                                                         >  1:2900129/1‑93 (MQ=255)

                cTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGc                                                                         <  2:1194265/93‑1 (MQ=255)

                 ttAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCt                                                                        >  1:3868804/1‑93 (MQ=255)

                 ttAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCt                                                                        >  1:7485376/1‑93 (MQ=255)

                 ttAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGACCTACAGCt                                                                        >  1:4686863/1‑93 (MQ=255)

                  tAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTc                                                                       <  2:7366523/93‑1 (MQ=255)

                  tAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTc                                                                       >  1:6681867/1‑93 (MQ=255)

                   aGCAAAATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGATTCCGGCAGGAAGGAGCGAACTAGGTCCTTCAGCTCg                                                                      >  2:6324230/1‑93 (MQ=255)

                     cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGt                                                                    <  1:5370851/93‑1 (MQ=255)

                      aaCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTc                                                                   >  1:20252/1‑93 (MQ=255)

                      aaCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTc                                                                   <  1:655962/93‑1 (MQ=255)

                      aaCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTACTCCAGCTCGGTc                                                                   >  2:7208528/1‑93 (MQ=255)

                       aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTcc                                                                  >  1:485599/1‑93 (MQ=255)

                       aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGT‑CTCCAGCTCGGTTCg                                                                 >  2:1314657/1‑93 (MQ=255)

                          tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGAt                                                               >  1:7122914/1‑93 (MQ=255)

                            aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATAc                                                             >  1:6603828/1‑93 (MQ=255)

                            aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATAc                                                             >  1:6902803/1‑93 (MQ=255)

                                 aaCCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTctgc                                                        <  1:3552927/93‑1 (MQ=255)

                                      gCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATa                                                   >  2:7581102/1‑93 (MQ=255)

                                      gCTCATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATa                                                   <  1:5228006/93‑1 (MQ=255)

                                        tGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGa                                                 >  1:1523336/1‑93 (MQ=255)

                                         gATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaa                                                >  1:3414318/1‑93 (MQ=255)

                                          aTCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaaa                                               >  1:7908755/1‑93 (MQ=255)

                                           acccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaaaa                                              <  1:6817233/92‑1 (MQ=255)

                                            cccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAAt                                             >  2:5020967/1‑93 (MQ=255)

                                            cccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAAt                                             <  1:6313469/93‑1 (MQ=255)

                                               gCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGt                                          >  2:5400071/1‑93 (MQ=255)

                                                 tCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGa                                        >  1:793649/1‑93 (MQ=255)

                                                 tCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGa                                        <  2:7191177/93‑1 (MQ=255)

                                                 tCATAAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGa                                        <  2:7113854/93‑1 (MQ=255)

                                                  cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAg                                       <  1:461234/93‑1 (MQ=255)

                                                  cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAg                                       <  2:3388924/93‑1 (MQ=255)

                                                   aTTAACCGCGACCCACAACGGCGCTTCCGGCAGTAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGt                                      <  1:5544264/93‑1 (MQ=255)

                                                    ttACCCGCGGCCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTg                                     <  2:3226537/93‑1 (MQ=255)

                                                     tAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTgg                                    <  1:1404156/93‑1 (MQ=255)

                                                       aCCGCGTCCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCa                                  <  1:3216098/93‑1 (MQ=255)

                                                       aCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCa                                  <  2:3798943/93‑1 (MQ=255)

                                                       aCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCa                                  <  1:976342/93‑1 (MQ=255)

                                                       aCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCa                                  >  1:4550706/1‑93 (MQ=255)

                                                       aCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCa                                  >  1:557883/1‑93 (MQ=255)

                                                         cgcgACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATc                                <  2:2005478/93‑1 (MQ=255)

                                                          gcgACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCg                               >  1:2570090/1‑93 (MQ=255)

                                                          gcgACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCg                               >  2:6021901/1‑93 (MQ=255)

                                                            gACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGAc                             <  2:7655845/93‑1 (MQ=255)

                                                             aCCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACa                            <  2:7846898/93‑1 (MQ=255)

                                                             aCCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACa                            >  2:2120713/1‑93 (MQ=255)

                                                             aCCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACa                            >  2:158623/1‑93 (MQ=255)

                                                              cccACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAg                           <  2:258667/93‑1 (MQ=255)

                                                               ccACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGa                          >  1:3985288/1‑93 (MQ=255)

                                                                 acaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAAc                        >  1:3332295/1‑93 (MQ=255)

                                                                   aaCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGa                      <  1:1584034/93‑1 (MQ=255)

                                                                   aaCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGa                      <  1:4141421/93‑1 (MQ=255)

                                                                   aaCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGa                      <  2:6243722/93‑1 (MQ=255)

                                                                    aCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGac                     >  2:3889371/1‑93 (MQ=255)

                                                                     cGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAAAGacc                    >  1:6015275/1‑92 (MQ=255)

                                                                     cGGCGCGTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCAGCGACAGAACGaca                    <  2:476678/93‑1 (MQ=255)

                                                                      ggCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGacac                   <  1:7566617/93‑1 (MQ=255)

                                                                       gcgcTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACAcc                  >  2:3530384/1‑93 (MQ=255)

                                                                       gcgcTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACAcc                  >  2:6593843/1‑93 (MQ=255)

                                                                        cgcTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCa                 >  2:1054935/1‑93 (MQ=255)

                                                                        cgcTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCa                 >  1:2808096/1‑93 (MQ=255)

                                                                          cTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCaaa               <  2:2074407/93‑1 (MQ=255)

                                                                           ttCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAAt              >  1:2401725/1‑93 (MQ=255)

                                                                              cGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGc           >  1:1212395/1‑93 (MQ=255)

                                                                                gCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCaa         >  2:5117780/1‑93 (MQ=255)

                                                                                gCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCaa         <  1:640227/93‑1 (MQ=255)

                                                                                 cAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAAc        <  1:5811090/93‑1 (MQ=255)

                                                                                   ggaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATCCTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGa      <  2:4414910/93‑1 (MQ=255)

                                                                                   ggaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGa      >  2:5686131/1‑93 (MQ=255)

                                                                                   ggaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGa      <  1:1025978/93‑1 (MQ=255)

                                                                                    gaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAg     >  2:6213847/1‑93 (MQ=255)

                                                                                     aaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGt    <  1:3038366/93‑1 (MQ=255)

                                                                                     aaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGt    <  2:2761292/93‑1 (MQ=255)

                                                                                      aggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGtt   >  1:5436864/1‑93 (MQ=255)

                                                                                       ggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTtc  >  2:6936912/1‑93 (MQ=255)

ATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTC  >  NC_000913/4211481‑4211660

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