BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000731	81	N→G	intergenic (–/‑151)	– / → repA	–/replication protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000731	81	0	N	G	99.9%	4906.9 / NA	1265	intergenic (–/‑151)	–/repA	–/replication protein A
Reads supporting (aligned to +/- strand): ref base N (0/0); major base G (683/581); minor base T (0/1); total (683/582)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1267 total aligned reads displayed.

‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAA  >  CP000731/1‑170

tataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGa                                                                                          >  2:130100‑M1/10‑91 (MQ=255)

tataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGa                                                                                          <  1:2130049‑M1/82‑1 (MQ=255)

 ataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaa                                                                                         >  2:4802058‑M1/9‑91 (MQ=255)

  taagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTAGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGaaa                                                                                        <  1:5106930‑M1/84‑1 (MQ=255)

   aagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAc                                                                                       >  1:3730032‑M1/7‑91 (MQ=255)

   aagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAc                                                                                       >  1:5749677‑M1/7‑91 (MQ=255)

    agtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACt                                                                                      <  1:4469396‑M1/86‑1 (MQ=255)

     gtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTg                                                                                     >  1:1663737‑M1/5‑91 (MQ=255)

       ttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtt                                                                                   <  1:4956209‑M1/89‑1 (MQ=255)

            tAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattt                                                                              <  2:1969047/91‑1 (MQ=255)

             aaaGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattta                                                                             >  1:2307478/1‑91 (MQ=255)

             aaaGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattta                                                                             <  1:5512263/91‑1 (MQ=255)

              aaGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTttatttat                                                                            >  2:1658273/1‑91 (MQ=255)

                gTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtatat                                                                          <  1:4040100/91‑1 (MQ=255)

                   cTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAt                                                                       >  2:4982109/1‑91 (MQ=255)

                   cTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAt                                                                       <  1:2705558/91‑1 (MQ=255)

                   cTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAt                                                                       >  1:1893069/1‑91 (MQ=255)

                   cTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAt                                                                       >  1:1724668/1‑91 (MQ=255)

                   cTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAt                                                                       >  2:4590361/1‑91 (MQ=255)

                     ataataAACTTAATATAACAAAAACCCCAACTACTGGTAATACTCGGGGTTTGGTAATATAACGTTTAGAAACTGTTTATTTATATAAttt                                                                     <  1:5675158/91‑1 (MQ=255)

                      taataaACTTAATATAACAAAAACCCCAACTACTGGTAATCGTCGGGGTTTGATAATAAAACGTTTAGAAACTGTTTATTTATACAAtcct                                                                    >  1:4853944/1‑88 (MQ=255)

                      taataaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtttt                                                                    >  1:5450122/1‑91 (MQ=255)

                      taataaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtttt                                                                    >  2:2556020/1‑91 (MQ=255)

                          aaaCTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGAt                                                                >  2:2788073/1‑91 (MQ=255)

                          aaaCTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGAt                                                                <  1:5312786/91‑1 (MQ=255)

                          aaaCTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGAt                                                                >  1:5203681/1‑91 (MQ=255)

                             cTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACt                                                             <  1:3353749/91‑1 (MQ=255)

                               tAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTg                                                           >  1:4561657/1‑91 (MQ=255)

                                aaTATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGt                                                          <  2:2202735/91‑1 (MQ=255)

                                aaTATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGt                                                          <  1:1577911/91‑1 (MQ=255)

                                aaTATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGt                                                          >  2:188345/1‑91 (MQ=255)

                                 atatAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGta                                                         <  1:4298052/91‑1 (MQ=255)

                                  tataACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtat                                                        >  2:5472988/1‑91 (MQ=255)

                                  tataACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtat                                                        <  2:3334256/91‑1 (MQ=255)

                                    taACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtatta                                                      >  2:2858279/1‑91 (MQ=255)

                                    taACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtatta                                                      <  2:5122899/91‑1 (MQ=255)

                                     aaCAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtattat                                                     <  2:3268126/91‑1 (MQ=255)

                                      aCAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAc                                                                                       >  2:2085266/1‑56 (MQ=255)

                                      aCAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtata                                                    >  2:3151334/1‑91 (MQ=255)

                                       cAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatat                                                   >  1:4782817/1‑91 (MQ=255)

                                       cAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatat                                                   >  2:3026321/1‑91 (MQ=255)

                                        aaaaaCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatata                                                  <  1:953175/91‑1 (MQ=255)

                                          aaaCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAAt                                                >  1:5814737/1‑91 (MQ=255)

                                          aaaCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAAt                                                <  2:2464335/91‑1 (MQ=255)

                                            aCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGa                                              >  2:1120607/1‑91 (MQ=255)

                                             ccccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaa                                             >  1:2808758/1‑91 (MQ=255)

                                             ccccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaa                                             >  1:3149182/1‑91 (MQ=255)

                                              cccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaat                                            <  1:862863/91‑1 (MQ=255)

                                                cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                          <  2:767394/91‑1 (MQ=255)

                                                cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                          <  2:646851/91‑1 (MQ=255)

                                                cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                          <  1:4347201/91‑1 (MQ=255)

                                                  actactGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAataata                                        <  2:3700637/91‑1 (MQ=255)

                                                     actGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCt                                     <  2:4191764/91‑1 (MQ=255)

                                                        ggTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAg                                  <  1:1420787/91‑1 (MQ=255)

                                                        ggTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAg                                  >  2:2954290/1‑91 (MQ=255)

                                                        ggTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAg                                  <  1:1026730/91‑1 (MQ=255)

                                                         gTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGt                                 >  2:5264924/1‑91 (MQ=255)

                                                          tAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTg                                <  1:2927983/91‑1 (MQ=255)

                                                             tAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCaa                             >  1:3020690/1‑91 (MQ=255)

                                                               gTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGc                           >  2:1228103/1‑91 (MQ=255)

                                                               gTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGc                           >  1:1097137/1‑91 (MQ=255)

                                                                tCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                     >  1:4182011/1‑48 (MQ=40)

                                                                tCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCa                          <  2:1569429/91‑1 (MQ=255)

                                                                tCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCa                          <  2:3605753/91‑1 (MQ=255)

                                                                 cGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAc                         >  1:2023675/1‑91 (MQ=255)

                                                                 cGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAc                         <  2:265297/91‑1 (MQ=255)

                                                                  ggggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTGTATAATGAATAATATCTTAGTGCAAGCAcc                        <  1:4632826/91‑1 (MQ=255)

                                                                  ggggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        >  2:3830018/1‑91 (MQ=255)

                                                                  ggggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        >  2:864367/1‑91 (MQ=255)

                                                                  ggggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        >  1:1248204/1‑91 (MQ=255)

                                                                   gggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCt                       >  2:3967096/1‑91 (MQ=255)

                                                                    ggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTa                      >  1:3548700/1‑91 (MQ=255)

                                                                     gTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTaa                     >  1:3426284/1‑91 (MQ=255)

                                                                     gTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTaa                     >  2:5752510/1‑91 (MQ=255)

                                                                       ttGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAAt                   >  2:239058/1‑91 (MQ=255)

                                                                        tGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATc                  >  1:2563965/1‑91 (MQ=255)

                                                                        tGGTAATAAAACGTTTAGAAACTGTCTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTACTGCAAGCAGCTAAATc                  >  2:5544287/1‑91 (MQ=255)

                                                                         ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCa                 <  2:348872/91‑1 (MQ=255)

                                                                         ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCa                 >  1:643765/1‑91 (MQ=255)

                                                                         ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCa                 >  1:74218/1‑91 (MQ=255)

                                                                         ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCa                 <  1:3301723/91‑1 (MQ=255)

                                                                         ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCa                 <  2:5679772/91‑1 (MQ=255)

                                                                         ggTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCa                 >  1:2399853/1‑91 (MQ=255)

                                                                          gTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCAt                <  2:4331873/91‑1 (MQ=255)

                                                                           taataaAACGTATAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCAtt               >  2:4414805/1‑91 (MQ=255)

                                                                            aataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATta              <  1:3937845/91‑1 (MQ=255)

                                                                            aataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATta              >  1:3215204/1‑91 (MQ=255)

                                                                             ataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaa             >  1:3827306/1‑91 (MQ=255)

                                                                             ataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaa             >  2:700862/1‑91 (MQ=255)

                                                                              taaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaa            <  2:1773761/91‑1 (MQ=255)

                                                                               aaaaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaat           >  2:4691150/1‑91 (MQ=255)

                                                                                 aaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataa         >  2:4901005/1‑91 (MQ=255)

                                                                                  aCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataaa        <  1:2495112/91‑1 (MQ=255)

                                                                                  aCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataaa        <  2:4076158/91‑1 (MQ=255)

                                                                                  aCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataaa        >  1:1801803/1‑91 (MQ=255)

                                                                                  aCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaataaa        >  1:4697679/1‑91 (MQ=255)

                                                                                     tttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGt     <  1:2223106/91‑1 (MQ=255)

                                                                                     tttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGt     >  1:1186350/1‑91 (MQ=255)

                                                                                      ttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTa    >  2:5409052/1‑91 (MQ=255)

                                                                                       tAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTaa   >  2:1413776/1‑91 (MQ=255)

‑‑‑‑‑‑‑‑‑ATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTANAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGTAAA  >  CP000731/1‑170

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