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.6%	451.7 / ‑2.8	145	A31A (GCT→GCC)	bioC	BioC
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (76/67); total (78/67)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 146 total aligned reads displayed.

GGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCG  >  NC_000913/810747‑810926

 gCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGAc                                                                                        <  1:3756979/93‑1 (MQ=255)

  cAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACg                                                                                       >  1:4212936/1‑93 (MQ=255)

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

     cgtttAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCt                                                                                    <  2:4200092/90‑1 (MQ=255)

     cgtttAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCt                                                                                    <  1:4569116/90‑1 (MQ=255)

     cGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCt                                                                                    >  1:2180526/1‑93 (MQ=255)

      gtttAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCtt                                                                                   <  1:2289118/91‑1 (MQ=255)

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

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

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

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

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

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

         taataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACt                                                                                >  1:3278230/1‑93 (MQ=255)

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

         taataaACAAGCCATTGCAGCGACATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACt                                                                                <  2:2890655/93‑1 (MQ=255)

          aataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTg                                                                               <  2:2528564/93‑1 (MQ=255)

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

          aagaaaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTg                                                                               <  1:3658485/90‑1 (MQ=255)

           ataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTgg                                                                              >  2:1403009/1‑93 (MQ=255)

             aaaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTGCAGCGCCAGAGTGCCGACGCCTTACTGGCa                                                                            <  1:2634999/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  1:4094187/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  2:1707930/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  1:2241685/93‑1 (MQ=255)

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

               aCAACCCATAGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAAt                                                                          <  2:1194554/93‑1 (MQ=255)

                cAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATg                                                                         >  1:4196234/1‑93 (MQ=255)

                 aaGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGc                                                                        <  2:1988630/93‑1 (MQ=255)

                  aGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCt                                                                       >  1:1951287/1‑93 (MQ=255)

                  aGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCt                                                                       <  2:2336931/93‑1 (MQ=255)

                  aGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAAATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCt                                                                       <  2:2831345/93‑1 (MQ=255)

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

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

                     cATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTAGCAATGCTTcc                                                                    >  2:3828802/1‑93 (MQ=255)

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

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

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

                         gCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAg                                                                <  2:4421902/93‑1 (MQ=255)

                          cAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGc                                                               <  1:3485884/93‑1 (MQ=255)

                          cAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGc                                                               <  2:2073986/93‑1 (MQ=255)

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

                             cGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTa                                                            >  2:4311073/1‑93 (MQ=255)

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

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

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

                                cATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAAt                                                         <  2:3892563/93‑1 (MQ=255)

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

                                 aTTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATa                                                        >  2:2400075/1‑93 (MQ=255)

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

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATac                                                       <  2:1861467/93‑1 (MQ=255)

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

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

                                  tttGGGCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATac                                                       >  2:3855936/1‑93 (MQ=255)

                                   tgggTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATaca                                                      <  2:3349655/91‑1 (MQ=255)

                                     ggTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACAcc                                                    >  2:1054502/1‑93 (MQ=255)

                                       tCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCa                                                  <  2:1342784/93‑1 (MQ=255)

                                       tCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCa                                                  <  1:4474432/93‑1 (MQ=255)

                                       tCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCa                                                  <  1:3831163/93‑1 (MQ=255)

                                         gggCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACg                                                >  2:4022314/1‑93 (MQ=255)

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

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

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

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

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

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

                                                cGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGa                                         <  1:4034803/93‑1 (MQ=255)

                                                  cacaCTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAcg                                       >  1:1003881/1‑93 (MQ=255)

                                                  cacaCTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCAACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAcg                                       >  2:1613726/1‑93 (MQ=255)

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

                                                      cTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCggg                                   <  1:4635920/93‑1 (MQ=255)

                                                       tATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGt                                  >  1:2954201/1‑93 (MQ=255)

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

                                                          gAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtgt                               >  2:2477308/1‑93 (MQ=255)

                                                          gAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTGACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtgt                               <  2:501989/93‑1 (MQ=255)

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

                                                               aCATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGGcc                          >  2:3853250/1‑93 (MQ=255)

                                                                 aTGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTCCTGGACGCGGGTTGTGGACctg                        <  1:2612067/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                        tCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGAt                 <  2:67661/93‑1 (MQ=255)

                                                                        tCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGAt                 <  1:2066952/93‑1 (MQ=255)

                                                                         cTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATg                <  2:4147182/93‑1 (MQ=255)

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

                                                                              gcgcCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAggcg           >  2:3981847/1‑93 (MQ=255)

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

                                                                               cgcCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccgc          <  1:3250460/93‑1 (MQ=255)

                                                                                gcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccgcc         >  2:2119295/1‑93 (MQ=255)

                                                                                gcCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccgcc         >  1:4403587/1‑93 (MQ=255)

                                                                                gcCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccgcc         >  2:4637369/1‑93 (MQ=255)

                                                                                gcCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccgcc         >  2:3069708/1‑93 (MQ=255)

                                                                                 ccAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCa        >  2:7301/1‑93 (MQ=255)

                                                                                   agagTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACt      >  1:4103327/1‑93 (MQ=255)

                                                                                   agagTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACt      >  2:647641/1‑93 (MQ=255)

                                                                                   agagTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACt      >  1:4204818/1‑93 (MQ=255)

                                                                                    gagTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTg     <  1:1387166/93‑1 (MQ=255)

                                                                                       tGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCg  >  1:1528361/1‑93 (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` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`