BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	894,155	T→G	100%	D124E (GAT→GAG)	potF →	putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	894,155	0	T	G	100.0%	362.9 / ‑6.7	105	D124E (GAT→GAG)	potF	putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base T (0/1); major base G (55/48); minor base A (1/0); total (56/49)

Only 100 of 105 total aligned reads displayed.

AGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGCCC  >  NC_000913/894058‑894251

aGTTGACTGTGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGcc                                                                                                <  1:535616/100‑1 (MQ=255)

aGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGcc                                                                                                <  1:210152/100‑1 (MQ=255)

aGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGcc                                                                                                <  2:1559942/100‑1 (MQ=255)

  ttGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCg                                                                                              <  2:2427859/100‑1 (MQ=255)

   tGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGa                                                                                             >  1:363987/1‑100 (MQ=255)

   tGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGa                                                                                             <  1:488427/100‑1 (MQ=255)

    gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                            <  2:795580/100‑1 (MQ=255)

    gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                            >  1:1210018/1‑100 (MQ=255)

    gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                            <  1:817381/100‑1 (MQ=255)

    gACTGCGGGAGGTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                            <  2:759615/100‑1 (MQ=255)

     aCTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                            <  2:1397120/99‑1 (MQ=255)

      cTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACaa                                                                                          >  1:202698/1‑100 (MQ=255)

        gCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATa                                                                                        <  1:1305245/100‑1 (MQ=255)

               ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                                 <  1:2333276/100‑1 (MQ=255)

               ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                                 >  1:2224020/1‑100 (MQ=255)

               ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                                 >  1:2423344/1‑100 (MQ=255)

               ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                                 <  2:1210010/100‑1 (MQ=255)

                tttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGGAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCt                                                                                <  2:1580057/100‑1 (MQ=255)

                tttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCt                                                                                >  1:2004205/1‑100 (MQ=255)

                tttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCt                                                                                <  2:383762/100‑1 (MQ=255)

                 tgccAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGGCAATAAATTTGCt                                                                                <  2:1781890/97‑1 (MQ=255)

                   ccAGCCGCTGGACGAAAGCAAATTGCCGGAGTGGAAGAGTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATg                                                                             >  2:1579624/1‑100 (MQ=255)

                    cAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGCACCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATGTACTATg                                                                             >  2:1782996/1‑99 (MQ=255)

                      gCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGccc                                                                          <  1:188377/100‑1 (MQ=255)

                        cGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCta                                                                        <  2:314761/100‑1 (MQ=255)

                        cGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCta                                                                        <  1:1140012/100‑1 (MQ=255)

                         gCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCtat                                                                       >  1:2061949/1‑100 (MQ=255)

                           tGGACAAAAGCAAATTGCCGGAGTGGAAGACTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTatat                                                                     <  1:1078409/100‑1 (MQ=255)

                           tGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTatat                                                                     <  2:288124/100‑1 (MQ=255)

                             gACAAAAGCAAATTGCCGGAATGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtg                                                                    >  2:1569687/1‑99 (MQ=255)

                              aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                                  >  2:1816637/1‑100 (MQ=255)

                              aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                                  >  1:1001435/1‑100 (MQ=255)

                                aaaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTggg                                                                <  1:1516374/100‑1 (MQ=255)

                                aaaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACAGCTCAAGCTGGTCGCCAAACAAGAGCCCGACAATAAATTTGCTATGCCCTATATGTgg                                                                 >  2:293691/1‑99 (MQ=255)

                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                               <  2:1461212/100‑1 (MQ=255)

                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                               >  1:800380/1‑100 (MQ=255)

                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                               >  2:1817391/1‑100 (MQ=255)

                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                               >  2:1324822/1‑100 (MQ=255)

                                 aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                               >  1:1721542/1‑100 (MQ=255)

                                  aaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcg                                                              <  2:194291/100‑1 (MQ=255)

                                  aaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcg                                                              >  2:1783510/1‑100 (MQ=255)

                                    gCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcgac                                                            >  2:1196395/1‑100 (MQ=255)

                                    gCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcgac                                                            >  1:161762/1‑100 (MQ=255)

                                    gCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcgac                                                            >  2:763797/1‑100 (MQ=255)

                                     cAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcgacg                                                           <  1:2240684/100‑1 (MQ=255)

                                        aTTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCAGAAGCGGGTCGCCAAACACGAGCCCGACAATAAATTAGCTATGCCCTATATGTGGGCgacgac                                                         >  2:1568081/1‑99 (MQ=255)

                                          tGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGTTATGCCCTATATGTGGGCGACGACCgg                                                      <  2:1388816/100‑1 (MQ=255)

                                            atggATTGGATGAATCTCGACCCGGAAATGCTGAAGCTGATATCCAAACAAGAGTGGGACAATAAATTAGTTAGGTCATATACGTGGGCGCCGACAgg                                                      >  2:1019197/3‑98 (MQ=255)

                                             ctgAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGAt                                                   <  2:1557125/98‑1 (MQ=255)

                                              ggAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGAtt                                                  >  1:1533765/1‑100 (MQ=255)

                                                 tctgaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                               <  1:139023/97‑1 (MQ=255)

                                                 gTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                               <  2:2167772/100‑1 (MQ=255)

                                                 gTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                               <  1:2349880/100‑1 (MQ=255)

                                                  tGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCt                                              <  2:2142361/100‑1 (MQ=255)

                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATACATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                             >  2:1228648/1‑100 (MQ=255)

                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                             >  1:823554/1‑100 (MQ=255)

                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                             >  1:1960306/1‑100 (MQ=255)

                                                   ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                             >  1:1327890/1‑100 (MQ=255)

                                                   ggAAGAATCTCGACACGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                             >  2:2290248/1‑100 (MQ=255)

                                                    gaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGcaga                                            >  2:139021/1‑97 (MQ=255)

                                                    gaagaaTCTCGACCCGGAACAGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtat                                            >  2:1571484/1‑100 (MQ=255)

                                                      agaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATaa                                          >  1:2080740/1‑100 (MQ=255)

                                                       gaaTCTCGACTCGGAACTGTTGAAGCTGGTTGCAAAACACGAACCCGACAACAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAAc                                         >  2:1999347/1‑100 (MQ=255)

                                                       gaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAAc                                         >  1:1112970/1‑100 (MQ=255)

                                                       gaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAAc                                         <  2:178181/100‑1 (MQ=255)

                                                           ctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTg                                     <  2:363983/100‑1 (MQ=255)

                                                               aCCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATaa                                 >  1:1228523/1‑100 (MQ=255)

                                                               aCCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATaa                                 >  1:1592599/1‑100 (MQ=255)

                                                                cccGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATaaa                                <  2:1533754/100‑1 (MQ=255)

                                                                  cGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGt                              <  1:1670766/100‑1 (MQ=255)

                                                                    gAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTa                            >  1:774953/1‑100 (MQ=255)

                                                                     aaCTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTaa                           >  1:1655680/1‑100 (MQ=255)

                                                                       ctgctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAg                         >  1:358991/1‑100 (MQ=255)

                                                                        tgctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGc                        <  1:2254348/100‑1 (MQ=255)

                                                                         gctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCg                       >  1:553315/1‑100 (MQ=255)

                                                                         gctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCg                       >  1:751647/1‑100 (MQ=255)

                                                                             aaGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGc                   <  1:986978/100‑1 (MQ=255)

                                                                              aGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCt                  >  1:1665542/1‑100 (MQ=255)

                                                                              aGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCt                  >  1:278574/1‑100 (MQ=255)

                                                                               ncTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTg                 >  1:2279580/2‑100 (MQ=255)

                                                                               gCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTg                 >  1:685468/1‑100 (MQ=255)

                                                                                  ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTggg               <  2:1183121/99‑1 (MQ=255)

                                                                                   gTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCg             <  2:803790/100‑1 (MQ=255)

                                                                                   gTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCg             <  2:2235004/100‑1 (MQ=255)

                                                                                    tCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa            >  2:432301/1‑100 (MQ=255)

                                                                                      gCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaa          >  1:1008871/1‑100 (MQ=255)

                                                                                       ccAAACACGAGCCCGACAATACATTTGCTATGCCCTATACGCGGGCGACGACCGAGATTGACTATAACATTGATAAAGTCAAAGCGTTGCTGGGCAaaaa         >  2:156887/1‑100 (MQ=255)

                                                                                       ccAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaaa         <  1:1138488/100‑1 (MQ=255)

                                                                                         aaaCACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAAcg       <  2:2004191/100‑1 (MQ=255)

                                                                                         aaaCACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAAcg       <  1:676429/100‑1 (MQ=255)

                                                                                             acGAGCCCGATAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGGTAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc   <  2:1678489/100‑1 (MQ=255)

                                                                                             acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgcn   <  1:571073/100‑2 (MQ=255)

                                                                                             acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgcn   <  1:1207104/100‑2 (MQ=255)

                                                                                             acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc   <  2:170208/100‑1 (MQ=255)

                                                                                             acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc   <  1:2550900/100‑1 (MQ=255)

                                                                                             acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc   <  1:2360264/100‑1 (MQ=255)

                                                                                             acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc   <  1:1749974/100‑1 (MQ=255)

                                                                                              cGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGccc  >  1:792785/1‑100 (MQ=255)

                                                                                              cGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGccc  >  1:1959165/1‑100 (MQ=255)

                                                                                              cGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGACTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGccc  >  1:1139951/1‑100 (MQ=255)

AGTTGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGCCC  >  NC_000913/894058‑894251

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