BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,434,055	A→G	P404P (CCA→CCG)	sun →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,434,055	0	A	G	99.5%	623.9 / ‑4.8	198	P404P (CCA→CCG)	sun	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (106/90); total (106/92)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 201 total aligned reads displayed.

CTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGG  >  NC_000913/3433965‑3434144

 tcttCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCggag                                                                                        <  2:147601/91‑1 (MQ=255)

    tCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCaa                                                                                     >  1:3840512/1‑93 (MQ=255)

     cGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAAc                                                                                    <  1:1494039/93‑1 (MQ=255)

        tgtTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCgg                                                                                 >  1:1000483/1‑93 (MQ=255)

            aCCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTaa                                                                             >  2:3196682/1‑93 (MQ=255)

              cGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAAc                                                                           <  2:609958/93‑1 (MQ=255)

               ggAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACa                                                                          >  1:1786121/1‑93 (MQ=255)

                gAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACaa                                                                         >  1:3025530/1‑93 (MQ=255)

                 aaGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACaaa                                                                        <  1:3968966/93‑1 (MQ=255)

                   gagaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAAt                                                                      >  1:3411250/1‑93 (MQ=255)

                   gagaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAAt                                                                      >  2:4605861/1‑93 (MQ=255)

                    agaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATc                                                                     >  2:3043896/1‑93 (MQ=255)

                    agaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATc                                                                     >  1:2865934/1‑93 (MQ=255)

                     gaATCGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCt                                                                    >  2:2257318/1‑93 (MQ=255)

                     gaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCt                                                                    >  2:1650606/1‑93 (MQ=255)

                     gaATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCt                                                                    >  2:1761130/1‑93 (MQ=255)

                      aaTAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTa                                                                   <  1:3441757/93‑1 (MQ=255)

                       aTAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTAc                                                                  >  1:1758226/1‑93 (MQ=255)

                       aTAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTAc                                                                  >  2:390764/1‑93 (MQ=255)

                       aTAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTAc                                                                  >  1:3270986/1‑93 (MQ=255)

                        tAGCCTGCAGATTAAATCCTTTTTGCAGCGTACCGCTGATGCCGATCTTTGCGAAACAGGAACACCTGAGCAACCGGGTAAACAAAATCTAcc                                                                 <  2:3804614/93‑1 (MQ=255)

                        tAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTAcc                                                                 <  2:775782/93‑1 (MQ=255)

                         aGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCt                                                                <  1:4170980/93‑1 (MQ=255)

                            cTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGt                                                             <  2:4178348/93‑1 (MQ=255)

                             tGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATATACCTGGTg                                                            >  2:682325/1‑93 (MQ=255)

                                 gATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGa                                                        >  2:3140807/1‑93 (MQ=255)

                                  aTTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGaa                                                       >  1:1838153/1‑93 (MQ=255)

                                    tAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAagg                                                     >  2:561210/1‑92 (MQ=255)

                                    tAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAaga                                                     <  2:3019047/93‑1 (MQ=255)

                                     aaaGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAagag                                                    >  1:2693229/1‑93 (MQ=255)

                                      aaGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAgg                                                   <  1:3788667/93‑1 (MQ=255)

                                      aaGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAgg                                                   <  1:2549849/93‑1 (MQ=255)

                                       aGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAggg                                                  >  2:1252851/1‑93 (MQ=255)

                                       aGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAggg                                                  <  2:1422824/93‑1 (MQ=255)

                                        gCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGc                                                 >  1:4507220/1‑93 (MQ=255)

                                        gCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGc                                                 <  1:4285791/93‑1 (MQ=255)

                                          cTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGa                                               >  1:2414116/1‑93 (MQ=255)

                                          cTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGa                                               >  2:1919674/1‑93 (MQ=255)

                                           tttttGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGAc                                              >  2:1619749/1‑93 (MQ=255)

                                           tttttGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGAc                                              >  1:2505331/1‑93 (MQ=255)

                                              ttGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACACAATCTACCTGGTGCCGAAGAGGGCGACGGc                                           >  1:795810/1‑93 (MQ=255)

                                              ttGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGc                                           <  2:1669344/93‑1 (MQ=255)

                                               tGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGct                                          <  1:3348044/93‑1 (MQ=255)

                                               tGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGct                                          >  2:4104055/1‑93 (MQ=255)

                                                gCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGctt                                         >  2:1041038/1‑93 (MQ=255)

                                                gCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGctt                                         <  2:3236502/93‑1 (MQ=255)

                                                   aCGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGcttctt                                      >  1:3664008/1‑93 (MQ=255)

                                                       tccGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTAc                                  <  1:2966967/92‑1 (MQ=255)

                                                       aCCGCTGATGCCGAACTTTGCGAAACAGGAAGACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTAc                                  >  2:3756616/1‑93 (MQ=255)

                                                       aCCGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTAc                                  <  1:4311367/93‑1 (MQ=255)

                                                        ccGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACAGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACg                                 >  1:3303399/1‑93 (MQ=255)

                                                         cGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgc                                >  1:1068506/1‑93 (MQ=255)

                                                         cGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgc                                >  2:196007/1‑93 (MQ=255)

                                                         cGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgc                                <  1:2035739/93‑1 (MQ=255)

                                                         cGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgc                                <  2:1328073/93‑1 (MQ=255)

                                                         cGCTGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgc                                <  1:4457133/93‑1 (MQ=255)

                                                            tGATGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaa                             <  1:1708531/93‑1 (MQ=255)

                                                               tGCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagct                          >  2:3615424/1‑93 (MQ=255)

                                                                gCCGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagcta                         <  1:3166611/93‑1 (MQ=255)

                                                                 ccGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagctaa                        >  2:4037736/1‑93 (MQ=255)

                                                                 ccGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagctaa                        <  1:1480440/93‑1 (MQ=255)

                                                                 ccGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagctaa                        >  1:4073617/1‑93 (MQ=255)

                                                                 ccGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACgctaagctaa                        <  1:857743/93‑1 (MQ=255)

                                                                  cGAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAAt                       >  1:1633319/1‑93 (MQ=255)

                                                                   gAACTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATc                      <  2:2787133/93‑1 (MQ=255)

                                                                    aaCTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCa                     >  1:2547301/1‑93 (MQ=255)

                                                                    aaCTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCa                     >  2:2683347/1‑93 (MQ=255)

                                                                     aCTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaa                    >  1:899050/1‑93 (MQ=255)

                                                                     aCTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaa                    >  2:3254804/1‑93 (MQ=255)

                                                                     aCTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaa                    >  1:1672937/1‑93 (MQ=255)

                                                                     aCTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaa                    >  2:1902532/1‑93 (MQ=255)

                                                                      cTTTGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaa                   >  1:3081595/1‑93 (MQ=255)

                                                                       tttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaa                  >  2:253532/1‑93 (MQ=255)

                                                                       tttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaa                  >  2:552345/1‑93 (MQ=255)

                                                                       tttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaa                  >  2:2094576/1‑93 (MQ=255)

                                                                        ttGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCaaaaa                 <  2:394833/93‑1 (MQ=255)

                                                                         tGCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAg                <  1:2925621/93‑1 (MQ=255)

                                                                          gCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGt               <  1:1534963/93‑1 (MQ=255)

                                                                          gCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGt               <  2:2281558/93‑1 (MQ=255)

                                                                          gCGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGt               <  1:1215120/93‑1 (MQ=255)

                                                                           cGAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtg              <  1:3886893/93‑1 (MQ=255)

                                                                            gAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtga             <  2:3299028/93‑1 (MQ=255)

                                                                            gAAACAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtga             >  2:4358899/1‑93 (MQ=255)

                                                                             aaaCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgat            <  1:3494981/93‑1 (MQ=255)

                                                                              aaCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgatg           <  1:2229556/93‑1 (MQ=255)

                                                                              aaCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgatg           >  2:4645982/1‑93 (MQ=255)

                                                                              aaCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgatg           <  2:2598946/93‑1 (MQ=255)

                                                                               aCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgatga          >  2:3398641/1‑93 (MQ=255)

                                                                               aCAGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGtgatga          <  2:1007485/93‑1 (MQ=255)

                                                                                 aGGAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATgaga        <  1:2515991/93‑1 (MQ=255)

                                                                                  ggAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGAt       <  1:966802/93‑1 (MQ=255)

                                                                                  ggAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGAt       <  1:2131294/93‑1 (MQ=255)

                                                                                  ggAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGAt       <  1:764752/93‑1 (MQ=255)

                                                                                  ggAACACCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGAt       <  1:1811683/93‑1 (MQ=255)

                                                                                     acacCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAAc    >  1:577877/1‑93 (MQ=255)

                                                                                     acacCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAAc    >  1:1118897/1‑93 (MQ=255)

                                                                                      cacCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACg   <  1:1976013/93‑1 (MQ=255)

                                                                                       acCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACgg  <  2:4462347/93‑1 (MQ=255)

                                                                                       acCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACgg  >  2:3987561/1‑93 (MQ=255)

                                                                                       acCGGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACgg  <  1:2811365/93‑1 (MQ=255)

CTGTTCGGTGTTACCGGAAGAGAATAGCCTGCAGATTAAAGCCTTTTTGCAACGTACCGCTGATGCCGAACTTTGCGAAACAGGAACACCAGAGCAACCGGGTAAACAAAATCTACCTGGTGCCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGG  >  NC_000913/3433965‑3434144

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