BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,005,601	+C	coding (237/282 nt)	b4568 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,005,600	1	.	C	99.0%	957.4 / 0.2	304	coding (236/282 nt)	b4568	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/3); new base C (154/147); total (154/150)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

CGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTGGC  >  NC_000913/4005508‑4005689

cgTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTgg                                                                                            >  1:4380135/1‑93 (MQ=255)

cgTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTgg                                                                                            >  2:1139079/1‑93 (MQ=255)

 gTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑c                                                                                          >  1:5453785/1‑93 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGCAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cc                                                                                         >  1:7134348/1‑92 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cc                                                                                         >  1:795407/1‑92 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cc                                                                                         >  2:3698129/1‑92 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cc                                                                                         >  2:5051125/1‑92 (MQ=255)

    aTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cctg                                                                                       <  2:1679114/93‑4 (MQ=255)

        tttGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGGGGCCTGTAtc                                                                                    >  1:1536118/1‑93 (MQ=255)

         ttGGGCTGGGATTGCAATGATTAGCCTGCATTTGATTATCGAAAACATGGCCTGGCTGCATACACTGATTATGGTGGGCGGTGGCCTGTAtct                                                                                   <  1:7048577/93‑1 (MQ=255)

          tGGGCTGGGATTGCGCTTCTTCGCCGGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtctc                                                                                  <  2:6147848/93‑1 (MQ=255)

          tGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTAtctc                                                                                  <  1:5845303/93‑1 (MQ=255)

            ggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctg                                                                                >  2:2278106/1‑93 (MQ=255)

               tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                                             >  2:1586626/1‑93 (MQ=255)

               tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                                             >  2:2567703/1‑93 (MQ=255)

               tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                                             <  2:6014837/93‑1 (MQ=255)

               tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                                             >  1:4809916/1‑93 (MQ=255)

                 ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                                           <  2:6337219/93‑1 (MQ=255)

                   aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                                         <  2:3741273/93‑1 (MQ=255)

                    ttGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATgg                                                                        <  2:1904308/93‑1 (MQ=255)

                     tGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATggg                                                                       >  1:3821447/1‑93 (MQ=255)

                      gcgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGt                                                                      >  2:7212898/1‑93 (MQ=255)

                      gcgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGt                                                                      <  1:2074186/93‑1 (MQ=255)

                       cgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGtt                                                                     >  2:3868953/1‑93 (MQ=255)

                        gctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTa                                                                    <  1:777125/93‑1 (MQ=255)

                        gctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTa                                                                    >  2:1634972/1‑93 (MQ=255)

                          tgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAcc                                                                  <  2:2199669/93‑1 (MQ=255)

                          tgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAcc                                                                  <  1:327624/93‑1 (MQ=255)

                          tgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAcc                                                                  >  1:6060139/1‑93 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCa                                                                 <  1:2962804/93‑1 (MQ=255)

                            ctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAg                                                                >  1:682758/1‑93 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                               <  1:2859321/93‑1 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                               <  2:3456921/93‑1 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                               <  1:6551203/93‑1 (MQ=255)

                              tgtccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGAt                                                              <  1:4238749/90‑1 (MQ=255)

                              tGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGAt                                                              <  1:3503200/93‑1 (MQ=255)

                               ggCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATg                                                             <  1:3319384/93‑1 (MQ=255)

                                  cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                          >  1:7423249/1‑93 (MQ=255)

                                  cTGCATTTGATTATCGAAAAAATGGCCTGGATGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATTCTa                                                          <  2:7639345/93‑1 (MQ=255)

                                   tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                         >  2:6795522/1‑93 (MQ=255)

                                   tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                         >  1:5268782/1‑93 (MQ=255)

                                   tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                         <  2:3033144/93‑1 (MQ=255)

                                   tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                         >  1:3542963/1‑93 (MQ=255)

                                     cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgt                                                       >  2:5285415/1‑93 (MQ=255)

                                     cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgt                                                       <  2:477977/93‑1 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                      >  1:5953281/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                      >  2:660965/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                      >  1:1216231/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                      <  1:5795213/93‑1 (MQ=255)

                                       tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtgg                                                     >  2:5613105/1‑93 (MQ=255)

                                        ttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggt                                                    >  2:4105183/1‑93 (MQ=255)

                                          gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                                                  <  2:3553883/93‑1 (MQ=255)

                                          gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                                                  <  2:3081940/93‑1 (MQ=255)

                                            ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                                >  2:1062025/1‑93 (MQ=255)

                                            ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                                >  1:1651429/1‑93 (MQ=255)

                                            ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                                <  2:2898073/93‑1 (MQ=255)

                                             tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                                               <  2:2913131/93‑1 (MQ=255)

                                              aTCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTg                                              >  2:5893874/1‑93 (MQ=255)

                                              aTCGAAAAAATGGACTGGCTGCATACGCTGCTTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTg                                              >  2:3114140/1‑93 (MQ=255)

                                                   aaaaaTGGTCTGGCTGTATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGGGCACGGtaaaa                                         >  2:7038218/1‑93 (MQ=255)

                                                    aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaaa                                        <  1:1956208/93‑1 (MQ=255)

                                                      aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                      <  1:612866/93‑1 (MQ=255)

                                                      aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                      <  2:143008/93‑1 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                     <  2:3362673/93‑1 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGGGGTGCACTGAAAAAagag                                     >  1:6573758/1‑93 (MQ=255)

                                                        tGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAgg                                    <  1:2758459/93‑1 (MQ=255)

                                                         ggCCGGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGc                                   <  1:1439130/93‑1 (MQ=255)

                                                          gCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg                                  >  2:538431/1‑93 (MQ=255)

                                                           ccTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCgg                                 <  2:7476326/93‑1 (MQ=255)

                                                            cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                                >  2:4438153/1‑93 (MQ=255)

                                                             tGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGtt                               >  1:6320454/1‑93 (MQ=255)

                                                              ggCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGttt                              >  2:4892633/1‑93 (MQ=255)

                                                              ggCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGttt                              >  1:1067468/1‑93 (MQ=255)

                                                              ggCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGgtt                              >  2:6965599/1‑93 (MQ=255)

                                                               gCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTc                             <  1:446867/93‑1 (MQ=255)

                                                                cTGCATACGCTGATTATGGTGGGCGGTGGCCTGTTTCTCTGCTGGATGGCTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCt                            <  1:6664308/93‑1 (MQ=255)

                                                                    aTACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAc                        >  2:7867231/1‑93 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                       >  2:4657403/1‑93 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                       >  1:963701/1‑93 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                       >  1:5609208/1‑93 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                       >  1:5121249/1‑93 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                       >  1:4898810/1‑93 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGGGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                       >  2:5566833/1‑93 (MQ=255)

                                                                      aCGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCt                      <  1:2652070/93‑1 (MQ=255)

                                                                      aCGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCt                      >  1:4635685/1‑93 (MQ=255)

                                                                         cTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcg                   >  1:1361232/1‑93 (MQ=255)

                                                                          tGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcgc                  >  2:6839773/1‑93 (MQ=255)

                                                                           gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGcc                 >  1:7477822/1‑93 (MQ=255)

                                                                           gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGcc                 >  1:6146942/1‑93 (MQ=255)

                                                                           gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGGGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGcc                 >  1:3128923/1‑93 (MQ=255)

                                                                             ttATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac               >  2:3178452/1‑93 (MQ=255)

                                                                             ttATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac               >  1:7624866/1‑93 (MQ=255)

                                                                              tATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcaca              >  1:2329343/1‑93 (MQ=255)

                                                                                tggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg            >  1:4056075/1‑93 (MQ=255)

                                                                                tggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg            >  2:1209407/1‑93 (MQ=255)

                                                                                tggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg            <  2:4326075/93‑1 (MQ=255)

                                                                                    gggCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGa        <  1:1110427/93‑1 (MQ=255)

                                                                                     ggCGGTGGCCTGCATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAg       <  1:2451671/93‑1 (MQ=255)

                                                                                      gCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGc      <  2:2409489/93‑1 (MQ=255)

                                                                                        ggTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTg    <  1:1905227/93‑1 (MQ=255)

CGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTGGC  >  NC_000913/4005508‑4005689

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