BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	810,837	T→C	A31A (GCT→GCC)	bioC →	BioC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	810,837	0	T	C	98.0%	469.9 / ‑0.3	153	A31A (GCT→GCC)	bioC	BioC
Reads supporting (aligned to +/- strand): ref base T (3/0); new base C (64/83); total (69/84)
Fisher's exact test for biased strand distribution p-value = 8.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01

Only 100 of 159 total aligned reads displayed.

GGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCG  >  NC_000913/810747‑810926

ggCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGa                                                                                         <  2:2085260/93‑1 (MQ=255)

 gCAGCGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGAc                                                                                        >  1:3907747/1‑93 (MQ=255)

  cAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACg                                                                                       <  2:2654872/93‑1 (MQ=255)

   aaCGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGc                                                                                      >  1:1646872/1‑93 (MQ=255)

    aCGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGcc                                                                                     <  1:2024749/93‑1 (MQ=255)

      ggTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCtt                                                                                   <  2:2592655/93‑1 (MQ=255)

       gTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTa                                                                                  >  2:4291632/1‑93 (MQ=255)

       gTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTa                                                                                  <  1:1820585/93‑1 (MQ=255)

       gTTAATAAACAAGCCATAGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTa                                                                                  >  2:3757110/1‑93 (MQ=255)

        ttAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTAc                                                                                 <  2:385878/93‑1 (MQ=255)

        ttAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTAc                                                                                 >  2:1934720/1‑93 (MQ=255)

         taataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACt                                                                                <  1:1207358/93‑1 (MQ=255)

         taataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACt                                                                                <  2:317731/93‑1 (MQ=255)

          aataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTg                                                                               >  2:3579826/1‑93 (MQ=255)

           ataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTgg                                                                              <  2:1571102/93‑1 (MQ=255)

           agaaaaaaGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTgg                                                                              <  2:4667727/91‑1 (MQ=255)

            taaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGc                                                                             >  1:2765788/1‑93 (MQ=255)

             aaaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCa                                                                            <  2:4456980/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCATAGTGCCGACGCCTTACTGGCaa                                                                           >  1:2345553/1‑93 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTGGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  2:213415/93‑1 (MQ=255)

               aCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAAt                                                                          <  1:357018/93‑1 (MQ=255)

               aCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAAt                                                                          <  1:3099821/93‑1 (MQ=255)

               aCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAAt                                                                          <  1:2935912/93‑1 (MQ=255)

                    ccATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTc                                                                     <  2:3137824/93‑1 (MQ=255)

                     cATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTcc                                                                    <  1:4515452/93‑1 (MQ=255)

                     cATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTcc                                                                    <  1:3703522/93‑1 (MQ=255)

                      aTTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCca                                                                   <  1:2025552/93‑1 (MQ=255)

                      aTTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCca                                                                   >  1:4019021/1‑93 (MQ=255)

                       ttGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCcac                                                                  <  1:51903/93‑1 (MQ=255)

                       ttGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCcac                                                                  >  2:2661817/1‑93 (MQ=255)

                        tGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCcaca                                                                 <  2:698305/93‑1 (MQ=255)

                         gCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAg                                                                >  1:132748/1‑93 (MQ=255)

                          cAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGTCAATGCTTCcacacc                                                               <  2:3292255/93‑3 (MQ=255)

                            gCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGt                                                             >  2:3686507/1‑93 (MQ=255)

                            gCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGt                                                             >  1:4493493/1‑93 (MQ=255)

                              ggCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTaa                                                           <  2:2429268/93‑1 (MQ=255)

                               gCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTaaa                                                          <  1:2107630/93‑1 (MQ=255)

                               gCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTaaa                                                          >  1:2981267/1‑93 (MQ=255)

                                 aTTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATa                                                        <  1:437239/93‑1 (MQ=255)

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATac                                                       >  1:4184266/1‑93 (MQ=255)

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATac                                                       >  2:1775045/1‑93 (MQ=255)

                                    tGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATacac                                                     >  1:3748543/1‑93 (MQ=255)

                                     ggTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTAAAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTGCCACAGCGTAAATACAcc                                                    <  2:1154384/93‑1 (MQ=255)

                                      gTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACAccc                                                   <  2:424659/93‑1 (MQ=255)

                                      gTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACAccc                                                   <  2:569360/93‑1 (MQ=255)

                                        cGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCAc                                                 <  1:3437624/93‑1 (MQ=255)

                                         gggCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACg                                                <  1:3190256/93‑1 (MQ=255)

                                          ggCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGt                                               >  2:2047191/1‑93 (MQ=255)

                                          ggCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGt                                               >  1:2171628/1‑93 (MQ=255)

                                            cAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTAc                                             <  1:4161647/93‑1 (MQ=255)

                                             aGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACt                                            <  1:4509665/93‑1 (MQ=255)

                                              gCCGCACACTATGAGCAACATGCGGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTg                                           <  1:1757491/93‑1 (MQ=255)

                                              gCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTg                                           >  1:345694/1‑93 (MQ=255)

                                                cGCACACTAGGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGa                                         <  2:1835923/93‑1 (MQ=255)

                                                 gCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAc                                        >  2:4629173/1‑93 (MQ=255)

                                                 gCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAc                                        <  2:177473/93‑1 (MQ=255)

                                                 gCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAc                                        <  2:3709773/93‑1 (MQ=255)

                                                 gCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAc                                        >  1:2202689/1‑93 (MQ=255)

                                                   acacTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAcgc                                      >  1:1724530/1‑93 (MQ=255)

                                                    cacTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAcgcg                                     >  2:4153637/1‑93 (MQ=255)

                                                    cacTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCAGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAcgcg                                     >  1:396543/1‑93 (MQ=255)

                                                      cTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCggg                                   >  1:1041207/1‑93 (MQ=255)

                                                        aTGAGCAACATGCAGATCTACAGTGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGtt                                 <  1:218859/93‑1 (MQ=255)

                                                         tGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtg                                <  1:1125376/93‑1 (MQ=255)

                                                         tGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtg                                <  2:2531349/93‑1 (MQ=255)

                                                         tGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtg                                >  1:2752152/1‑93 (MQ=255)

                                                         tGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtg                                >  2:2023668/1‑93 (MQ=255)

                                                         tGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtg                                <  2:124342/93‑1 (MQ=255)

                                                         tGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtg                                >  2:3815411/1‑93 (MQ=255)

                                                           aGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtgtg                              >  2:4370759/1‑93 (MQ=255)

                                                            gCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTgg                             <  2:1552447/93‑1 (MQ=255)

                                                            gCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTgg                             >  1:800711/1‑93 (MQ=255)

                                                               aCATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGAcc                          <  2:515065/93‑1 (MQ=255)

                                                                 aTGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctg                        >  1:1559901/1‑93 (MQ=255)

                                                                  tGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctgg                       >  1:627720/1‑93 (MQ=255)

                                                                  tGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctgg                       >  2:3038710/1‑93 (MQ=255)

                                                                   gCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctggc                      <  2:1383859/93‑1 (MQ=255)

                                                                     aGATCTACATCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctggctg                    >  1:958230/1‑93 (MQ=255)

                                                                     aGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctggctg                    >  2:1243063/1‑93 (MQ=255)

                                                                     aGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctggctg                    <  1:1446610/93‑1 (MQ=255)

                                                                      gATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctggctgg                   >  1:2853233/1‑93 (MQ=255)

                                                                       aTCTACAGCTCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGa                  <  1:627409/93‑1 (MQ=255)

                                                                       aTCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGa                  >  2:3355361/1‑93 (MQ=255)

                                                                        tCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGAt                 >  2:729918/1‑93 (MQ=255)

                                                                          tACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGa               >  1:181726/1‑93 (MQ=255)

                                                                             aGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgcc            >  1:734822/1‑93 (MQ=255)

                                                                              gcgcCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccg           >  1:1477389/1‑93 (MQ=255)

                                                                               cgcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccgc          >  2:2355089/1‑93 (MQ=255)

                                                                                  cAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCAc       >  1:2521351/1‑93 (MQ=255)

                                                                                  cAGAGTGCAGACGCCTTACTTGCAATGCTTCCACAGCGTAAAGACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCAc       <  1:1372539/93‑1 (MQ=255)

                                                                                   agagTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACt      >  2:4325831/1‑93 (MQ=255)

                                                                                   agagTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACt      >  1:1561843/1‑93 (MQ=255)

                                                                                    gagaGCCGAGGCCTTGCTGGAAAGGCTTGCACCGCGTAAATACATCCACGTACTGGACGCGGGATGTGGACCTGGCCGGATGAGCCGTCACTg     >  2:812004/1‑93 (MQ=255)

                                                                                    gagTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTg     >  2:1832873/1‑93 (MQ=255)

                                                                                     agTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTgg    <  1:2519107/93‑1 (MQ=255)

                                                                                     agTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTgg    <  1:2559176/93‑1 (MQ=255)

                                                                                      gTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGTATGAGCCGCCACTGGc   <  2:1108370/93‑1 (MQ=255)

                                                                                      gTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGc   <  2:4165175/93‑1 (MQ=255)

                                                                                      gTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGc   <  2:602484/93‑1 (MQ=255)

                                                                                       tGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCg  <  2:3836514/93‑1 (MQ=255)

GGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCG  >  NC_000913/810747‑810926

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