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%	5843.2 / NA	1514	intergenic (–/‑151)	–/repA	–/replication protein A
Reads supporting (aligned to +/- strand): ref base N (0/0); major base G (794/719); minor base C (0/1); total (794/720)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1515 total aligned reads displayed.

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

tataagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGa                                                                                          >  1:2809026‑M1/10‑91 (MQ=255)

   aagtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAc                                                                                       >  2:2903951‑M1/7‑91 (MQ=255)

   aagtttATTTAAAGTGCTAGAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAAc                                                                                       >  2:3249170‑M1/7‑91 (MQ=255)

    agtttATTTAAAGTGCTATAATAAACTTAATATCACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACt                                                                                      <  1:6552436‑M1/86‑1 (MQ=255)

     gtttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTg                                                                                     <  2:3377716‑M1/87‑1 (MQ=255)

      tttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGt                                                                                    <  1:1228689‑M1/88‑1 (MQ=255)

      tttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGt                                                                                    <  1:5646811‑M1/88‑1 (MQ=255)

       ttATTTAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtt                                                                                   >  2:5360632‑M1/3‑91 (MQ=255)

            tAAAGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGtttattt                                                                              >  2:3713942/1‑91 (MQ=255)

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

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

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

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

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

               aGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtata                                                                           >  1:1003588/1‑91 (MQ=255)

               aGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtata                                                                           >  1:755733/1‑91 (MQ=255)

               aGTGCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTtata                                                                           >  2:810311/1‑91 (MQ=255)

                  gCTATAATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATaa                                                                        <  1:5178580/91‑1 (MQ=255)

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

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

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

                    tataATAAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtt                                                                      >  2:5110253/1‑91 (MQ=255)

                     ataataAACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAttt                                                                     >  2:5508400/1‑91 (MQ=255)

                      taataaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAAtttt                                                                    <  2:1759430/91‑1 (MQ=255)

                      taataaACTTAATATAAC‑AAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTc                                                                   <  2:4460884/91‑1 (MQ=255)

                        ataaACTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCg                                                                  <  1:180433/91‑1 (MQ=255)

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

                            aCTTAATATAACAAAACCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATAc                                                              <  1:291420/91‑1 (MQ=255)

                             cTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTACTAAAACGCTTAGAAACTGTTTATTTATATAATTTTCGATACt                                                             >  1:1897807/1‑91 (MQ=255)

                             cTTAATATAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACt                                                             >  1:1565016/1‑91 (MQ=255)

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

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

                                 atatAACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTCTCGATACTTGta                                                         >  2:4200781/1‑91 (MQ=255)

                                   ataACAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtatt                                                       >  2:2568066/1‑91 (MQ=255)

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

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

                                     aaCAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtattat                                                     >  1:5005771/1‑91 (MQ=255)

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

                                     aaCAAAAACCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtattat                                                     >  1:3973318/1‑91 (MQ=255)

                                     aaCAAAAACCCCAACTAATGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGtattat                                                     <  2:1394414/91‑1 (MQ=255)

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

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

                                        aaaaaCCCCAGCTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatata                                                  <  2:5246776/91‑1 (MQ=255)

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

                                        aaaaaCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATtatata                                                  >  1:632911/1‑91 (MQ=255)

                                         aaaaCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATaa                                                 >  1:930156/1‑91 (MQ=255)

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

                                            aCCCCAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTATATTATATAATGa                                              >  2:4677295/1‑91 (MQ=255)

                                              cccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaat                                            >  1:2677669/1‑91 (MQ=255)

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

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

                                              cccAACGACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaat                                            <  2:3823290/91‑1 (MQ=255)

                                               ccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaata                                           >  1:4283812/1‑91 (MQ=255)

                                               ccAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaata                                           <  2:4769471/91‑1 (MQ=255)

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

                                                cAACTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataa                                          >  2:6366698/1‑91 (MQ=255)

                                                 aaCTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataat                                         >  1:3837424/1‑91 (MQ=255)

                                                 aaCTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataat                                         >  2:362022/1‑91 (MQ=255)

                                                 aaCTACTGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGaataat                                         >  1:2376116/1‑91 (MQ=255)

                                                  actactGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAataata                                        >  2:117275/1‑91 (MQ=255)

                                                    tactGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATc                                      <  1:3642653/91‑1 (MQ=255)

                                                    tactGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATc                                      <  2:2032115/91‑1 (MQ=255)

                                                      ctcgtagtaGTCGGGGTTCGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCtt                                    <  1:5516970/88‑1 (MQ=255)

                                                      ctGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCtt                                    >  2:3907897/1‑91 (MQ=255)

                                                       tGGTAATAGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTa                                   <  1:3252259/91‑1 (MQ=255)

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

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

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

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

                                                              aGTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAg                            <  1:622174/91‑1 (MQ=255)

                                                               gTCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGc                           <  2:2971843/91‑1 (MQ=255)

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

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

                                                                tCGGGGTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCa                          >  2:5702111/1‑91 (MQ=255)

                                                                  ggggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCAcc                        <  1:6151484/91‑1 (MQ=255)

                                                                   gggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCt                       <  2:6002958/91‑1 (MQ=255)

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

                                                                   gggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCt                       <  2:2329813/91‑1 (MQ=255)

                                                                   gggTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTCTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCt                       <  2:5854221/91‑1 (MQ=255)

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

                                                                     gTTTGGTAATAAAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTaa                     <  2:189066/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                             ataaAACGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaa             <  1:45319/91‑1 (MQ=255)

                                                                                aaaCGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATtaaata          >  1:1980296/1‑91 (MQ=255)

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

                                                                                   cGTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTaaataaat       >  2:2475800/1‑91 (MQ=255)

                                                                                    gTTTAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCGTAGTGCAAGCACCTAAATCATTAAATAAATg      <  1:6440374/91‑1 (MQ=255)

                                                                                     tttAGAAACTGTTTATTTATATAATTTTCGATACTTGTATTATATAATGAATAATATCTTAGTGCAAGCACCTAAATCATTAAATAAATGt     >  1:1363020/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: `‑`