BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,610,653	A→G	P18P (CCA→CCG)	yhhU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,610,653	0	A	G	86.9%	758.1 / 93.7	314	P18P (CCA→CCG)	yhhU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (24/17); new base G (177/96); total (201/113)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 323 total aligned reads displayed.

GATACGCACTTTCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCACCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACC  >  NC_000913/3610561‑3610745

  tacgCACTTTCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTgccgcc                                                                                            >  1:2489951/1‑93 (MQ=255)

     gCACTTTCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCggg                                                                                         <  1:5542497/93‑1 (MQ=255)

     gCACTTTCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCggg                                                                                         >  2:3429339/1‑93 (MQ=255)

       aCTTTCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGtt                                                                                       >  1:3041867/1‑93 (MQ=255)

         tttCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTa                                                                                     >  1:7292625/1‑93 (MQ=255)

         tttCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTa                                                                                     >  1:579927/1‑93 (MQ=255)

         tttCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTa                                                                                     <  1:2798858/93‑1 (MQ=255)

          ttCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTAc                                                                                    <  2:838979/93‑1 (MQ=255)

          ttCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTAc                                                                                    >  1:6783521/1‑93 (MQ=255)

          ttCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTAc                                                                                    >  2:1020457/1‑93 (MQ=255)

            cATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTAcgc                                                                                  <  2:5113072/93‑1 (MQ=255)

            cATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTAcgc                                                                                  <  2:4421908/93‑1 (MQ=255)

             aTTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTAcgcg                                                                                 >  1:1908758/1‑93 (MQ=255)

             aTTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTAcgcg                                                                                 >  2:5030316/1‑93 (MQ=255)

              ttttCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGa                                                                                >  2:1314987/1‑93 (MQ=255)

              ttttCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGa                                                                                <  2:5868423/93‑1 (MQ=255)

              ttttCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGa                                                                                <  1:3461412/93‑1 (MQ=255)

               tttCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGag                                                                               >  1:4006291/1‑93 (MQ=255)

                 tCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGagca                                                                             >  1:1424333/1‑93 (MQ=255)

                  ccATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGagcaa                                                                            >  2:7518258/1‑93 (MQ=255)

                  ccATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGagcaa                                                                            <  2:1541773/93‑1 (MQ=255)

                   cATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGagcaag                                                                           >  1:6303774/1‑93 (MQ=255)

                   cATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGGTTACGCGagcaag                                                                           >  1:174294/1‑93 (MQ=255)

                    aTTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGagcaagc                                                                          <  1:2937545/93‑1 (MQ=255)

                     ttAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGagcaagca                                                                         >  1:4705549/1‑93 (MQ=255)

                      tAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCAc                                                                        <  2:6297159/93‑1 (MQ=255)

                      tAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCAc                                                                        <  1:1240510/93‑1 (MQ=255)

                       aaaCGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCAcc                                                                       >  1:1519286/1‑93 (MQ=255)

                         aCGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGc                                                                     <  1:6664915/93‑1 (MQ=255)

                         aCGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGc                                                                     <  2:6457639/93‑1 (MQ=255)

                          cGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCa                                                                    <  2:4083268/93‑1 (MQ=255)

                          cGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCa                                                                    >  1:7325001/1‑93 (MQ=255)

                           gTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCaa                                                                   <  1:7357509/93‑1 (MQ=255)

                           gTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCACCGGGTTTACGCGAGCAAGCACCGCag                                                                   >  2:2201672/1‑92 (MQ=255)

                             tGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAgg                                                                 <  2:5979520/93‑1 (MQ=255)

                              ggCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGt                                                                <  2:400222/93‑1 (MQ=255)

                               gCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTc                                                               >  2:5263294/1‑93 (MQ=255)

                               gCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTc                                                               >  2:3117223/1‑93 (MQ=255)

                                cccTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTcc                                                              >  2:6139714/1‑93 (MQ=255)

                                 ccTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCa                                                             <  2:7330151/93‑1 (MQ=255)

                                 ccTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCa                                                             <  1:270741/93‑1 (MQ=255)

                                  cTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCac                                                            >  1:3777822/1‑93 (MQ=255)

                                   tGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCacg                                                           >  1:1097693/1‑93 (MQ=255)

                                   tGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCacg                                                           >  1:4997558/1‑93 (MQ=255)

                                     atatGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCacgac                                                         >  1:7498812/1‑93 (MQ=255)

                                      tatGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCacgacg                                                        >  1:7463532/1‑93 (MQ=255)

                                       atGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGAcgc                                                       >  1:4882690/1‑93 (MQ=255)

                                       atGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGAcgc                                                       <  2:2836189/93‑1 (MQ=255)

                                        tGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGAcgcg                                                      >  2:494887/1‑93 (MQ=255)

                                        tGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGAcgcg                                                      >  1:4425865/1‑93 (MQ=255)

                                          tATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGaa                                                    <  1:6938344/93‑1 (MQ=255)

                                          tATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGaa                                                    <  1:5973698/93‑1 (MQ=255)

                                           aTCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAAc                                                   >  1:6425203/1‑93 (MQ=255)

                                            tCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACg                                                  >  1:1644836/1‑93 (MQ=255)

                                            tCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACg                                                  >  1:5178819/1‑93 (MQ=255)

                                             cGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgc                                                 <  2:6007298/93‑1 (MQ=255)

                                             cGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgc                                                 >  2:268700/1‑93 (MQ=255)

                                              ggATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgct                                                <  2:4144721/93‑1 (MQ=255)

                                              ggATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgct                                                <  2:1709124/93‑1 (MQ=255)

                                               gATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgctg                                               >  2:2472194/1‑93 (MQ=255)

                                               gATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgctg                                               >  1:3361544/1‑93 (MQ=255)

                                               gATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgctg                                               >  2:6597230/1‑93 (MQ=255)

                                                 tAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgctggc                                             >  1:4370437/1‑93 (MQ=255)

                                                   gTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgctggctg                                           >  2:7478782/1‑93 (MQ=255)

                                                   gTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgctggctg                                           <  1:6506867/93‑1 (MQ=255)

                                                   gTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACgctggctg                                           >  2:5561579/1‑93 (MQ=255)

                                                    ttCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGctggctgg                                          >  2:1344105/1‑93 (MQ=255)

                                                      cTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCg                                        <  1:4203696/93‑1 (MQ=255)

                                                       tGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCgg                                       <  1:2128468/93‑1 (MQ=255)

                                                        ggggAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCggg                                      <  1:2205862/93‑1 (MQ=255)

                                                        ggggAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCggg                                      >  1:5325019/1‑93 (MQ=255)

                                                         gggAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCgggg                                     <  1:6109718/93‑1 (MQ=255)

                                                         gggAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCgggg                                     >  2:3749274/1‑93 (MQ=255)

                                                         gggAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCgggg                                     >  1:2406929/1‑93 (MQ=255)

                                                          ggAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGc                                    >  1:7220269/1‑93 (MQ=255)

                                                           gAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCg                                   >  2:1148/1‑93 (MQ=255)

                                                            aaaGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGt                                  <  1:3433136/93‑1 (MQ=255)

                                                             aaGTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTg                                 >  1:2075386/1‑93 (MQ=255)

                                                               gTTTCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCa                               <  2:3939936/93‑1 (MQ=255)

                                                                 ttCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCAtt                             >  2:2774354/1‑93 (MQ=255)

                                                                  tCGACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTg                            <  1:4554443/93‑1 (MQ=255)

                                                                    gACCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCt                          <  2:1909714/93‑1 (MQ=255)

                                                                     aCCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCtt                         <  1:915977/93‑1 (MQ=255)

                                                                     aCCTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCtt                         <  2:3849667/93‑1 (MQ=255)

                                                                      ccTTAAGCGCAGCTCCACTGCCACCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCttt                        >  2:4278647/1‑93 (MQ=255)

                                                                      ccTTAAGCGCAGCTCCACTGCCACCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCttt                        >  2:3352537/1‑93 (MQ=255)

                                                                       cTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTc                       >  2:5451103/1‑93 (MQ=255)

                                                                       cTTAAGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTc                       >  2:6230891/1‑93 (MQ=255)

                                                                        ttAAGCGCAGCTCCACTGCCACCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCg                      >  2:6702694/1‑93 (MQ=255)

                                                                          aaGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGca                    <  2:1267241/93‑1 (MQ=255)

                                                                           aGCGCAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGcac                   >  2:465078/1‑93 (MQ=255)

                                                                              gcAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACg                >  2:3095993/1‑93 (MQ=255)

                                                                              gcAGCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACg                <  2:2307066/93‑1 (MQ=255)

                                                                                aGCTCCACTGCCGCAGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCCTTGCTTTCGCACACGCt              >  1:7143866/1‑93 (MQ=255)

                                                                                 gCTCCACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCtt             >  1:3593857/1‑93 (MQ=255)

                                                                                   tCCACTGCCACCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTc           >  2:7007831/1‑93 (MQ=255)

                                                                                     cACTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTccc         <  1:2909100/93‑1 (MQ=255)

                                                                                       cTGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCg       >  2:303066/1‑93 (MQ=255)

                                                                                        tGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGc      >  1:6800878/1‑93 (MQ=255)

                                                                                        tGCCGCCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGc      >  1:1022978/1‑93 (MQ=255)

GATACGCACTTTCATTTTCCATTAAACGTTGGCCCTGATATGTATCGGATAGTTCTGGGGAAAGTTTCGACCTTAAGCGCAGCTCCACTGCCACCGGGTTTACGCGAGCAAGCACCGCAAGGTCCACGACGCGAACGCTGGCTGGCGGGGCGTGCATTGCTTTCGCACACGCTTTCCCCGCTACC  >  NC_000913/3610561‑3610745

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