BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	654,812	G→A	A127V (GCC→GTC)	dcuC ←	DcuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	654,812	0	G	A	99.3%	1001.6 / ‑4.7	283	A127V (GCC→GTC)	dcuC	DcuC
Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (121/159); total (123/160)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 289 total aligned reads displayed.

ATACCAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGGCGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCATTCGC  >  NC_000913/654722‑654900

aTACCAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACg                                                                                        <  1:5542260/93‑1 (MQ=255)

aTACCAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACg                                                                                        <  2:7612428/93‑1 (MQ=255)

   ccAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACa                                                                                     <  1:2361730/93‑1 (MQ=255)

    cAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGaaaa                                                                                    >  1:1647663/1‑93 (MQ=255)

    cAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACaa                                                                                    >  2:2702057/1‑93 (MQ=255)

    cAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACaa                                                                                    <  2:7159625/93‑1 (MQ=255)

     aaCGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACaaa                                                                                   >  1:3999132/1‑93 (MQ=255)

      aCGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACaaaa                                                                                  >  1:2894220/1‑93 (MQ=255)

       cGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAAt                                                                                 <  1:6664846/93‑1 (MQ=255)

        gTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATa                                                                                >  1:7738790/1‑93 (MQ=255)

        gTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATa                                                                                >  2:6977378/1‑93 (MQ=255)

         tttACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAg                                                                               >  1:4990068/1‑93 (MQ=255)

         tttACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAAAAg                                                                               >  2:6647838/1‑93 (MQ=255)

         tttACCATCACCGGAAATAGGGTGGCCATCAGCAAAACTCCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAg                                                                               <  2:6227389/93‑1 (MQ=255)

          ttACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAgg                                                                              >  1:3449843/1‑93 (MQ=255)

           tgccATCACCGGAAATAGGGTTGCCAGCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAggc                                                                             <  2:3486395/91‑1 (MQ=255)

           tACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAggc                                                                             <  1:5098184/93‑1 (MQ=255)

                tCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCaa                                                                        >  2:5725726/1‑93 (MQ=255)

                  aCCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatc                                                                      >  2:6786796/1‑93 (MQ=255)

                  aCCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAAGCGGCAatc                                                                      >  2:7842420/1‑93 (MQ=255)

                   cctgAAATAGGGTTGCCATCAGCAAACCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatca                                                                     <  1:3132568/90‑1 (MQ=255)

                   ccGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatca                                                                     >  2:6754925/1‑93 (MQ=255)

                     ggAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatcatc                                                                   <  1:4480690/93‑1 (MQ=255)

                     ggAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatcatc                                                                   <  1:2115724/93‑1 (MQ=255)

                     ggAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatcatc                                                                   <  2:310561/93‑1 (MQ=255)

                     ggAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAatcatc                                                                   <  2:1040339/93‑1 (MQ=255)

                      gAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAAtcatca                                                                  <  1:3739170/93‑1 (MQ=255)

                        aaTAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAAAAGGCGGCAATCATcaga                                                                >  2:4732965/1‑92 (MQ=255)

                          tAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATcagcag                                                              <  1:7372086/93‑1 (MQ=255)

                           aGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAgg                                                             <  2:1366155/93‑1 (MQ=255)

                            gggTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGt                                                            >  2:1189561/1‑93 (MQ=255)

                             ggTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTa                                                           <  1:1894598/93‑1 (MQ=255)

                              gTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTaa                                                          <  1:2280872/93‑1 (MQ=255)

                               ttGCCATCAGCAACACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAg                                                         <  1:807201/93‑1 (MQ=255)

                               ttGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAg                                                         >  2:4968228/1‑93 (MQ=255)

                               ttGCAATCAGCAAGACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAg                                                         <  2:2874003/93‑1 (MQ=255)

                                tGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAgg                                                        >  1:4578528/1‑93 (MQ=255)

                                tGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAgg                                                        <  1:7071953/93‑1 (MQ=255)

                                 gCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAggg                                                       >  2:3570020/1‑93 (MQ=255)

                                  ccATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAgggg                                                      <  1:1689266/93‑1 (MQ=255)

                                   cATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGa                                                     >  1:4902365/1‑93 (MQ=255)

                                   cATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGa                                                     <  1:6227579/93‑1 (MQ=255)

                                   cATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGa                                                     <  1:5805818/93‑1 (MQ=255)

                                   cATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGa                                                     <  1:1946057/93‑1 (MQ=255)

                                     tCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGt                                                   >  2:3968652/1‑93 (MQ=255)

                                     tCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGt                                                   >  2:5098457/1‑93 (MQ=255)

                                      cAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGtt                                                  <  2:4552720/93‑1 (MQ=255)

                                      cAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGtt                                                  <  1:2025679/93‑1 (MQ=255)

                                       aGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTa                                                 >  2:1427620/1‑93 (MQ=255)

                                        gCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTaa                                                <  2:927771/93‑1 (MQ=255)

                                         cAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAat                                               >  2:302935/1‑93 (MQ=255)

                                          aaaaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAata                                              >  2:5497968/1‑93 (MQ=255)

                                           aaaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAatat                                             <  1:7686194/93‑1 (MQ=255)

                                            aaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAAtata                                            <  1:127957/93‑1 (MQ=255)

                                            aaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAAtata                                            <  2:5572241/93‑1 (MQ=255)

                                            aaCACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAAtata                                            >  2:7646994/1‑93 (MQ=255)

                                             acacCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATAc                                           <  2:41553/93‑1 (MQ=255)

                                               acCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTg                                         >  2:6596100/1‑93 (MQ=255)

                                               acCCAGACCGGGTGCGGAAGAGACGGCAAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTg                                         <  2:7264078/93‑1 (MQ=255)

                                                cccAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCCGGTAAGGGGAGTTAATATACTGc                                        >  2:5380102/1‑93 (MQ=255)

                                                cccAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGc                                        <  2:4905737/93‑1 (MQ=255)

                                                 ccAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCa                                       >  2:134730/1‑93 (MQ=255)

                                                  cAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCaa                                      <  2:3401380/93‑1 (MQ=255)

                                                  cAGACCGGTTGAGGAAGAGATGGCCAGAGAAATCATACAGGCGACAAAATAGGGGGCAATCATCAGCAGGAAAGGGGAGATAATATACTGCaa                                      >  1:3649618/1‑93 (MQ=255)

                                                   aGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAAt                                     <  2:2217984/93‑1 (MQ=255)

                                                   aGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAAt                                     <  1:2479656/93‑1 (MQ=255)

                                                   aGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAAt                                     <  1:4783926/93‑1 (MQ=255)

                                                   aGACCGGTGGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAAt                                     <  2:1629474/93‑1 (MQ=255)

                                                     aCCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATgg                                   >  2:4700167/1‑93 (MQ=255)

                                                     aCCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATgg                                   >  1:6370234/1‑93 (MQ=255)

                                                       cGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGtt                                 >  1:626283/1‑93 (MQ=255)

                                                       cGGCTGCGGAAGAGACGCCCAGAGACATCAGGCAGTCTACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCCATGGtt                                 <  1:33849/93‑1 (MQ=255)

                                                        ggTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGttt                                <  2:7691609/93‑1 (MQ=255)

                                                         gTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGtttt                               >  1:1080822/1‑93 (MQ=255)

                                                         gTTGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGtttt                               >  2:6024465/1‑93 (MQ=255)

                                                          ttGCGGAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTg                              >  1:4541409/1‑93 (MQ=255)

                                                               gAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGcc                         >  1:6041864/1‑93 (MQ=255)

                                                               gAAGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGcc                         <  2:5507258/93‑1 (MQ=255)

                                                                aaGAGACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGGAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCa                        >  1:3227039/1‑93 (MQ=255)

                                                                 agagACGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAg                       >  1:5270312/1‑93 (MQ=255)

                                                                    gaCGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCtt                    >  1:6979682/1‑93 (MQ=255)

                                                                     aCGGCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTg                   <  2:2045423/93‑1 (MQ=255)

                                                                       ggCCAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGac                 <  1:718948/93‑1 (MQ=255)

                                                                         ccAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGacca               <  2:3238143/93‑1 (MQ=255)

                                                                          cAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGaccac              >  1:569752/1‑93 (MQ=255)

                                                                          cAGAGACATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGaccac              >  2:5547052/1‑93 (MQ=255)

                                                                            gagaCATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGAccacca            <  1:3521362/93‑1 (MQ=255)

                                                                            gagaCATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGAccacca            >  1:4347566/1‑93 (MQ=255)

                                                                            gagaCATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGAccacca            >  2:4245426/1‑93 (MQ=255)

                                                                            gagaCATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGAccacca            >  1:1542702/1‑93 (MQ=255)

                                                                                cATCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATc        >  1:267036/1‑93 (MQ=255)

                                                                                 aTCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCa       <  1:6859558/93‑1 (MQ=255)

                                                                                 aTCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCa       <  1:3263532/93‑1 (MQ=255)

                                                                                 aTCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCa       <  2:4471360/93‑1 (MQ=255)

                                                                                 aTCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCa       >  2:824435/1‑93 (MQ=255)

                                                                                  tCAGACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCAt      >  2:3816562/1‑93 (MQ=255)

                                                                                   cagacagaCGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCAtt     <  2:6538763/93‑1 (MQ=255)

                                                                                     gACAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCATTcg   <  2:3870115/93‑1 (MQ=255)

                                                                                      aCAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCATTcgc  <  2:6444280/93‑1 (MQ=255)

                                                                                      aCAGACGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCATTcgc  <  2:1768314/93‑1 (MQ=255)

ATACCAACGTTTACCATCACCGGAAATAGGGTTGCCATCAGCAAAACACCCAGACCGGTTGCGGAAGAGACGGCCAGAGACATCAGACAGGCGACAAAATAGGCGGCAATCATCAGCAGGTAAGGGGAGTTAATATACTGCAATGGTTTTGACGCCAGCTTGACCACCATATCATTCGC  >  NC_000913/654722‑654900

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