BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,247,622	A→G	T164A (ACG→GCG)	malM →	MalM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,247,622	0	A	G	87.9%	656.7 / 72.5	264	T164A (ACG→GCG)	malM	MalM
Reads supporting (aligned to +/- strand): ref base A (8/24); new base G (123/109); total (131/133)
Fisher's exact test for biased strand distribution p-value = 4.12e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 266 total aligned reads displayed.

CCAGGCGTGATGAGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCC  >  NC_000913/4247532‑4247711

ccAGGCGTGATGAGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACcgcg                                                                                         >  2:5036923/1‑93 (MQ=255)

 cAGGCGTGATGAGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCgg                                                                                        >  1:7345691/1‑93 (MQ=255)

 cAGGCGTGATGAGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCgg                                                                                        >  2:225615/1‑93 (MQ=255)

 cAGGCGTGATGAGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACgg                                                                                        >  1:6334448/1‑93 (MQ=255)

    gCGTGATGAGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGaaa                                                                                     >  2:6529307/1‑93 (MQ=255)

      gTGATGAGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGaaaaa                                                                                   >  1:3269120/1‑93 (MQ=255)

      gTGATGAGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACAGCGGaaaaa                                                                                   >  1:3836642/1‑93 (MQ=255)

       tgatgaGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAg                                                                                  >  1:1089784/1‑93 (MQ=255)

       tgatgaGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAg                                                                                  <  2:1921966/93‑1 (MQ=255)

        gatgaGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGa                                                                                 >  1:1350975/1‑93 (MQ=255)

        gatgaGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGa                                                                                 >  1:5857758/1‑93 (MQ=255)

          tgaGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGAtc                                                                               <  1:6874565/93‑1 (MQ=255)

          tgaGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGAtc                                                                               <  1:4546952/93‑1 (MQ=255)

          tgaGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGAtc                                                                               >  1:3562582/1‑93 (MQ=255)

           gaGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGAtct                                                                              <  1:4059560/93‑1 (MQ=255)

            aGTGCAGATCGGCTGGAAGGCGTTATGGGCCTGACACCGGGGTTGGGGCAGCAAAAACTTTTTTTTCTGG‑CTTTACCGCGGAAAAAGATCTcc                                                                            <  1:2395347/93‑1 (MQ=255)

             gTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTcc                                                                            <  1:3408602/93‑1 (MQ=255)

             gTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTcc                                                                            >  1:3641937/1‑93 (MQ=255)

             gTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTcc                                                                            <  2:7050223/93‑1 (MQ=255)

              tGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCca                                                                           <  1:4019948/93‑1 (MQ=255)

               gCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCcag                                                                          <  2:1105199/93‑1 (MQ=255)

                 aGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCcagca                                                                        >  2:1684838/1‑93 (MQ=255)

                 aGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCcagca                                                                        <  2:3725321/93‑1 (MQ=255)

                  gATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCcagcag                                                                       >  1:6055170/1‑93 (MQ=255)

                  gATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCcagcag                                                                       >  1:4906287/1‑93 (MQ=255)

                    tCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAgac                                                                     <  1:3122884/93‑1 (MQ=255)

                    tCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAgac                                                                     <  1:7406704/93‑1 (MQ=255)

                     cGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAgacg                                                                    <  1:6456058/93‑1 (MQ=255)

                      ggCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAgacga                                                                   >  2:1331912/1‑93 (MQ=255)

                      ggCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAgacga                                                                   >  1:5783406/1‑93 (MQ=255)

                       gCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAgacgac                                                                  <  2:3628368/93‑1 (MQ=255)

                       gCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAgacgac                                                                  >  2:2519403/1‑93 (MQ=255)

                       gCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAgacgac                                                                  >  2:4018105/1‑93 (MQ=255)

                        cTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGAcc                                                                 <  1:5499433/93‑1 (MQ=255)

                         tGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGAccc                                                                <  1:5374750/93‑1 (MQ=255)

                          ggAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCa                                                               >  2:207112/1‑93 (MQ=255)

                           gAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCaa                                                              >  1:3946441/1‑93 (MQ=255)

                           gAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCaa                                                              <  1:1023871/93‑1 (MQ=255)

                            aaGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAAc                                                             <  2:6307805/93‑1 (MQ=255)

                            aaGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAAc                                                             <  1:1613636/93‑1 (MQ=255)

                             aGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACt                                                            >  1:5653387/1‑93 (MQ=255)

                             aGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACt                                                            <  2:6271145/93‑1 (MQ=255)

                             aGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACt                                                            <  1:3253413/93‑1 (MQ=255)

                              ggCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTg                                                           >  1:6720533/1‑93 (MQ=255)

                                cGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCt                                                         <  2:7506021/93‑1 (MQ=255)

                                  ttATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCg                                                       >  2:1850010/1‑93 (MQ=255)

                                   tATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGa                                                      >  2:4352785/1‑93 (MQ=255)

                                     tGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGAcc                                                    >  2:5348440/1‑93 (MQ=255)

                                         ccTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGc                                                >  1:3058711/1‑93 (MQ=255)

                                          cTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCt                                               >  2:1043363/1‑93 (MQ=255)

                                             acacCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTaaa                                            <  2:5799766/93‑1 (MQ=255)

                                             acacCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTaaa                                            <  1:5901215/93‑1 (MQ=255)

                                              cacCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAg                                           <  1:2030233/93‑1 (MQ=255)

                                                 cGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCt                                        <  2:3453077/93‑1 (MQ=255)

                                                  ggCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTa                                       <  1:629813/93‑1 (MQ=255)

                                                  ggCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTa                                       <  1:615591/93‑1 (MQ=255)

                                                  ggCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTa                                       <  2:3960500/93‑1 (MQ=255)

                                                   gCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTAt                                      >  2:5789164/1‑93 (MQ=255)

                                                   gCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTAt                                      <  1:3712221/93‑1 (MQ=255)

                                                   gCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTAt                                      <  1:4609156/93‑1 (MQ=255)

                                                     gTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGc                                    >  2:4878922/1‑93 (MQ=255)

                                                     gTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGc                                    <  1:6191957/93‑1 (MQ=255)

                                                     gTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGc                                    >  2:38823/1‑93 (MQ=255)

                                                       tGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCa                                  <  2:4776707/93‑1 (MQ=255)

                                                          ggCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAgg                               <  2:6991198/93‑1 (MQ=255)

                                                          ggCAGCAAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAgg                               <  1:3424171/93‑1 (MQ=255)

                                                           gcagcaAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAggg                              >  2:1419109/1‑93 (MQ=255)

                                                           gcagcaAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAggg                              >  1:2791256/1‑93 (MQ=255)

                                                           gcagcaAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAggg                              >  1:7124187/1‑93 (MQ=255)

                                                              gcaAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGt                           <  2:3205778/93‑1 (MQ=255)

                                                              gcaAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGt                           >  1:1541609/1‑93 (MQ=255)

                                                              gcaAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGt                           >  2:1222664/1‑93 (MQ=255)

                                                              gcaAAAACTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGt                           >  2:5164265/1‑93 (MQ=255)

                                                                aaaaaCTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCg                         >  2:2595865/1‑93 (MQ=255)

                                                                    aCTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTaa                     >  2:2965705/1‑93 (MQ=255)

                                                                    aCTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTaa                     <  2:5519982/93‑1 (MQ=255)

                                                                     cTTTATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAAc                    <  1:1245294/93‑1 (MQ=255)

                                                                      tttATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACt                   >  2:618847/1‑93 (MQ=255)

                                                                      tttATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACt                   <  2:6861807/93‑1 (MQ=255)

                                                                       ttATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTc                  <  2:3032525/93‑1 (MQ=255)

                                                                       ttATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTc                  <  2:6759441/93‑1 (MQ=255)

                                                                       ttATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTc                  <  1:1914881/93‑1 (MQ=255)

                                                                        tATGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCg                 >  1:803716/1‑93 (MQ=255)

                                                                          tGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGAt               >  2:5455049/1‑93 (MQ=255)

                                                                          tGTTCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGAt               >  1:4721742/1‑93 (MQ=255)

                                                                          tGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGAt               <  2:4687359/93‑1 (MQ=255)

                                                                           gTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATc              >  1:2168353/1‑93 (MQ=255)

                                                                            ttCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATcc             <  2:5712686/93‑1 (MQ=255)

                                                                            ttCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATcc             <  1:6904849/93‑1 (MQ=255)

                                                                             tCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATccc            <  1:4807684/93‑1 (MQ=255)

                                                                             tCTGGTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATccc            <  2:240689/93‑1 (MQ=255)

                                                                                ggTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGa         <  2:109003/93‑1 (MQ=255)

                                                                                 gTCTTTACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGat        >  2:2290/1‑93 (MQ=255)

                                                                                    tttACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCCAGGGCGTCGGTAACTCGATCCCGGATATc     >  2:7729/1‑93 (MQ=255)

                                                                                    tttACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATc     >  2:1078760/1‑93 (MQ=255)

                                                                                    tttACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATc     >  2:6828687/1‑93 (MQ=255)

                                                                                    tttACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATc     >  1:2620139/1‑93 (MQ=255)

                                                                                      tACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATccc   <  2:7282262/93‑1 (MQ=255)

                                                                                      tACCGCGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATccc   >  2:926440/1‑93 (MQ=255)

                                                                                      tACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATccc   <  1:1688723/93‑1 (MQ=255)

CCAGGCGTGATGAGTGCAGATCGGCTGGAAGGCGTTATGCGCCTGACACCGGCGTTGGGGCAGCAAAAACTTTATGTTCTGGTCTTTACCACGGAAAAAGATCTCCAGCAGACGACCCAACTGCTCGACCCGGCTAAAGCCTATGCCAAGGGCGTCGGTAACTCGATCCCGGATATCCCC  >  NC_000913/4247532‑4247711

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: `‑`