BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,363	G→A	C189C (TGC→TGT)	menA ←	MenA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,363	0	G	A	99.6%	852.1 / NA	231	C189C (TGC→TGT)	menA	MenA
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (142/88); total (142/89)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

GCGCACCACCAGCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCGCATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGCCAAAGAAAAC  >  NC_000913/4117273‑4117452

   caccacCAGCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGc                                                                                      <  2:3030204/93‑1 (MQ=255)

    accaccAGCGTGTTTTTGCCAGTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCg                                                                                     >  2:4008062/1‑93 (MQ=255)

       accaGCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTc                                                                                  >  2:5087139/1‑93 (MQ=255)

       aacaGCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTc                                                                                  >  1:4525719/3‑93 (MQ=255)

         caGCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGc                                                                                >  2:1527366/1‑93 (MQ=255)

           gCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCg                                                                              >  1:614308/1‑93 (MQ=255)

           gCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCg                                                                              >  2:5866399/1‑93 (MQ=255)

           gCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCg                                                                              <  1:4052865/93‑1 (MQ=255)

           gCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCg                                                                              >  2:1189666/1‑93 (MQ=255)

           gCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCg                                                                              <  2:4293495/93‑1 (MQ=255)

             gtgtTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGa                                                                            >  2:1555844/1‑93 (MQ=255)

             gtgtTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGa                                                                            >  1:1367761/1‑93 (MQ=255)

              tgtTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGaa                                                                           >  1:5897732/1‑93 (MQ=255)

               gtTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAg                                                                          >  2:4802009/1‑93 (MQ=255)

                tttttGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAgg                                                                         >  1:5159100/1‑93 (MQ=255)

                 ttttGCCATTTTCGCGTTCTCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGa                                                                        <  1:4695511/93‑1 (MQ=255)

                 ttttGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTCCTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGa                                                                        <  2:2969951/93‑1 (MQ=255)

                 ttttGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGa                                                                        <  1:4625000/93‑1 (MQ=255)

                  tttGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGAt                                                                       <  2:5077001/93‑1 (MQ=255)

                    tGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCa                                                                     <  2:4912416/93‑1 (MQ=255)

                    tGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCa                                                                     >  1:5138785/1‑93 (MQ=255)

                    tGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCa                                                                     >  2:5443111/1‑93 (MQ=255)

                     gCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAg                                                                    <  2:2689859/93‑1 (MQ=255)

                     gCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAg                                                                    >  1:4913363/1‑93 (MQ=255)

                      ccATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGt                                                                   <  2:6297349/93‑1 (MQ=255)

                      ccATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGt                                                                   <  1:2115197/93‑1 (MQ=255)

                       cATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTg                                                                  >  2:3640640/1‑93 (MQ=255)

                        aTTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGc                                                                 >  2:3336085/1‑93 (MQ=255)

                        aTTTTCGCGGTCGCTACTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGc                                                                 <  2:4406659/93‑1 (MQ=255)

                         ttttCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGcc                                                                >  1:357811/1‑93 (MQ=255)

                           ttCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCgg                                                              <  2:2036361/93‑1 (MQ=255)

                           ttCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCgg                                                              <  1:6297026/93‑1 (MQ=255)

                                ggTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCa                                                         <  2:2393862/93‑1 (MQ=255)

                                ggTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGAACAGTGCCGGAATCa                                                         >  1:2604079/1‑93 (MQ=255)

                                 gTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCaa                                                        <  2:321984/93‑1 (MQ=255)

                                 gTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCaa                                                        >  2:3389414/1‑93 (MQ=255)

                                  tCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAAt                                                       <  1:5646833/93‑1 (MQ=255)

                                   cGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATg                                                      <  2:2443691/93‑1 (MQ=255)

                                    gCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGGCGGAATCTATga                                                     >  2:4094260/1‑92 (MQ=255)

                                    gCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGt                                                     >  1:3894304/1‑93 (MQ=255)

                                      tATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTAt                                                   >  1:5847624/1‑93 (MQ=255)

                                      tATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTAt                                                   <  1:787189/93‑1 (MQ=255)

                                       aTTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATg                                                  >  2:3537649/1‑93 (MQ=255)

                                       aTTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATg                                                  <  1:2850410/93‑1 (MQ=255)

                                       aTTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATg                                                  >  1:6155176/1‑93 (MQ=255)

                                        ttGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGa                                                 >  2:665253/1‑93 (MQ=255)

                                        ttGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGa                                                 <  1:3673733/93‑1 (MQ=255)

                                            tatCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCt                                             >  2:4495336/1‑93 (MQ=255)

                                              tCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGt                                           <  1:5795693/93‑1 (MQ=255)

                                              tCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGt                                           >  2:1275279/1‑93 (MQ=255)

                                               cACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTa                                          >  2:1870061/1‑93 (MQ=255)

                                               cACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTa                                          <  2:3032547/93‑1 (MQ=255)

                                                 cGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTaaa                                        <  2:2555202/93‑1 (MQ=255)

                                                 cGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTaaa                                        >  2:830260/1‑93 (MQ=255)

                                                   cAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATa                                      <  1:3479592/93‑1 (MQ=255)

                                                    aGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATAc                                     >  1:2867360/1‑93 (MQ=255)

                                                    aGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATAc                                     <  2:4829155/93‑1 (MQ=255)

                                                     ggTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATAcc                                    >  2:3789263/1‑93 (MQ=255)

                                                     ggTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATAcc                                    >  1:4202839/1‑93 (MQ=255)

                                                     ggTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATAcc                                    >  1:1481993/1‑93 (MQ=255)

                                                     ggTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATAcc                                    >  1:5638346/1‑93 (MQ=255)

                                                     ggTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATAcc                                    <  2:5691650/93‑1 (MQ=255)

                                                      gTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCa                                   <  2:5339417/93‑1 (MQ=255)

                                                      gTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCa                                   >  1:4721273/1‑93 (MQ=255)

                                                        tattaATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGc                                 <  1:4083032/93‑1 (MQ=255)

                                                         attaATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCt                                >  1:3163962/1‑93 (MQ=255)

                                                         attaATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCt                                >  1:6269266/1‑93 (MQ=255)

                                                          ttaATATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTc                               >  2:2267790/1‑93 (MQ=255)

                                                            aaTATTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTccc                             >  2:779293/1‑93 (MQ=255)

                                                             atatTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTcccc                            <  2:2513182/93‑1 (MQ=255)

                                                              tatTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCa                           <  1:1662545/93‑1 (MQ=255)

                                                              tatTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCa                           >  2:3582046/1‑93 (MQ=255)

                                                               atTCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCAt                          >  1:1601089/1‑93 (MQ=255)

                                                                ttCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATg                         >  1:2904640/1‑93 (MQ=255)

                                                                 tCAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGa                        <  1:3089028/93‑1 (MQ=255)

                                                                  cAGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGac                       <  1:1038219/93‑1 (MQ=255)

                                                                  cAGTACTGCCGTTGCCAGCAGGCCACATGCGCTCGCCGGAAGGATCAGTCCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGac                       <  2:4441497/93‑1 (MQ=255)

                                                                   aGTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATTTATGAGCCTGTAAATACCAGCTCCCCATGaca                      >  1:2431199/1‑93 (MQ=255)

                                                                    gTACTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGacac                     <  1:6199829/93‑1 (MQ=255)

                                                                      aCTGCCGTTGCCAGCAGGCGACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTc                   >  2:6166144/1‑93 (MQ=255)

                                                                       cTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCa                  >  2:3488098/1‑93 (MQ=255)

                                                                       cTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCa                  >  1:435380/1‑93 (MQ=255)

                                                                       cTGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCa                  >  2:4109882/1‑93 (MQ=255)

                                                                        tGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCaa                 <  1:128241/93‑1 (MQ=255)

                                                                        tGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCaa                 >  1:521983/1‑93 (MQ=255)

                                                                        tGCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCaa                 >  2:277914/1‑93 (MQ=255)

                                                                         gCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAAc                >  2:1752256/1‑93 (MQ=255)

                                                                         gCCGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAAc                >  2:5358745/1‑93 (MQ=255)

                                                                           cGTTGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCa              <  1:2622747/93‑1 (MQ=255)

                                                                           cGTGGCAAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCa              <  2:3266804/93‑1 (MQ=255)

                                                                             ttGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGc            >  1:1943011/1‑93 (MQ=255)

                                                                             ttGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGc            >  1:6085098/1‑93 (MQ=255)

                                                                             ttGCCAGCAGGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGc            >  1:4414347/1‑93 (MQ=255)

                                                                                 cagcagGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGCCaaa        >  2:599307/1‑93 (MQ=255)

                                                                                  agcagGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGCCAAAg       >  1:2309283/1‑93 (MQ=255)

                                                                                  agcagGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGCCAAAg       <  1:3279389/93‑1 (MQ=255)

                                                                                   gcagGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGCCAAAGa      >  1:5395955/1‑93 (MQ=255)

                                                                                    cagGCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGCCAAAGaa     <  1:4863293/93‑1 (MQ=255)

                                                                                      ggCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGCCAAAGaaac   >  1:2092155/1‑92 (MQ=255)

                                                                                      ggCCACATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGCCAAAGaaaa   >  2:6185654/1‑93 (MQ=255)

GCGCACCACCAGCGTGTTTTTGCCATTTTCGCGGTCGCTATTGATATCACGCAGGTTATTAATATTCAGTACTGCCGTTGCCAGCAGGCCGCATGCGGTCGCCGGAAGGATCAGTGCCGGAATCAATGTATGAGCCTGTAAATACCAGCTCCCCATGACACTCAACCAGCCAAAGAAAAC  >  NC_000913/4117273‑4117452

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