BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,397,618	T→G	I60L (ATT→CTT)	rng ←	ribonuclease G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,397,618	0	T	G	99.3%	463.3 / ‑4.0	138	I60L (ATT→CTT)	rng	ribonuclease G
Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (72/64); total (73/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 138 total aligned reads displayed.

CGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAATATG  >  NC_000913/3397522‑3397714

tgcgcACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATc                                                                                               <  2:1603640/99‑1 (MQ=255)

cgcgCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATc                                                                                               <  2:1652861/100‑1 (MQ=255)

 gcgcACCGCGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCt                                                                                              <  2:844698/100‑1 (MQ=255)

  cgcACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTa                                                                                             >  1:2965733/1‑100 (MQ=255)

  cgcACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTa                                                                                             >  1:3504886/1‑100 (MQ=255)

   gcACCGTGAATTGCTTTTGTTCTTCACCCGCCACACCTTCTGTATGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTAc                                                                                            >  1:3185881/1‑100 (MQ=255)

   gcACCGTGAATTGCTGTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTAc                                                                                            <  2:3305080/100‑1 (MQ=255)

     aCCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACaa                                                                                          >  1:1744181/1‑100 (MQ=255)

     aCCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACaa                                                                                          >  2:3162010/1‑100 (MQ=255)

      ccGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACaaa                                                                                         >  1:2257384/1‑100 (MQ=255)

      ccGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACaaa                                                                                         >  1:3670008/1‑100 (MQ=255)

      ccGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACaaa                                                                                         >  1:1235760/1‑100 (MQ=255)

       cGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACaaaa                                                                                        >  2:2734480/1‑100 (MQ=255)

           aaTTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAgccg                                                                                    >  1:2660119/1‑100 (MQ=255)

             ttGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAgccgcc                                                                                  <  2:360461/100‑1 (MQ=255)

             cggCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAgccgcc                                                                                  <  2:1401863/98‑1 (MQ=255)

               gCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTg                                                                                >  2:399382/1‑100 (MQ=255)

               gCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTg                                                                                <  2:316713/100‑1 (MQ=255)

                cgtttGTTCTTCACCCCCCGCACAGTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGc                                                                               <  1:707019/98‑1 (MQ=255)

                cTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGc                                                                               >  1:1466378/1‑100 (MQ=255)

                 ttttGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCa                                                                              <  1:3392225/100‑1 (MQ=255)

                    taatcttcACACGCCACACATTCGGTTTGCGGCCTGATGTCGGATGCATGAAGAAACGCGGCTTTATCCATCCCAAGGTATACAAAAGCCGCCAGCAtgc                                                                           >  2:1770984/4‑98 (MQ=255)

                    tGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATAc                                                                           >  2:757829/1‑100 (MQ=255)

                     gTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATAcc                                                                          <  1:2627087/100‑1 (MQ=255)

                     gTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATAcc                                                                          >  1:288131/1‑100 (MQ=255)

                     gTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATAcc                                                                          >  1:820344/1‑100 (MQ=255)

                      ttcttcACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATAccc                                                                         >  1:2818387/1‑100 (MQ=255)

                      ttcttcACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCCTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATAccc                                                                         >  2:810206/1‑100 (MQ=255)

                       tcttcACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCg                                                                        >  1:154795/1‑100 (MQ=255)

                       tcttcACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCAGGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCg                                                                        >  2:230113/1‑100 (MQ=255)

                          tcACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGaa                                                                     >  2:2051181/1‑100 (MQ=255)

                           cACCCGCTACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCGAGATCTACAAAAGCCGCCTGCATACCCGGAAg                                                                    <  2:221735/100‑1 (MQ=255)

                           cACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAg                                                                    >  2:2486429/1‑100 (MQ=255)

                            aCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGt                                                                   >  1:3305749/1‑100 (MQ=255)

                            aCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGt                                                                   >  1:3419745/1‑100 (MQ=255)

                            aCCCGCCACACATTCGGTGTGCGGCATGAAGTCGGATGCATGAAGAAACGCGGCTTTAACCAGCCCAAGATCTACAAAAGCCGCCTGCATATCCGGAAGt                                                                   >  1:3906657/1‑100 (MQ=255)

                             cccGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTa                                                                  >  1:2587734/1‑100 (MQ=255)

                              ccGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTac                                                                 >  1:274144/1‑100 (MQ=255)

                              ccGCCACACATGCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTac                                                                 <  2:1970573/100‑1 (MQ=255)

                                  cacacaTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGa                                                             >  2:2800046/1‑100 (MQ=255)

                                   acacaTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGAc                                                            >  1:2227444/1‑100 (MQ=255)

                                      caTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACCACTTa                                                         <  2:3521040/100‑1 (MQ=255)

                                      caTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCAAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTa                                                         >  2:3179058/1‑100 (MQ=255)

                                       aTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGACGCACACGACTTac                                                        >  1:2119515/1‑100 (MQ=255)

                                       aTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTac                                                        >  2:3845365/1‑100 (MQ=255)

                                       aTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTac                                                        <  2:1858481/100‑1 (MQ=255)

                                         tCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTacac                                                      <  1:3022148/100‑1 (MQ=255)

                                             tgtgCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGAcc                                                  <  1:309933/100‑1 (MQ=255)

                                              gtgCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGAccc                                                 <  1:38605/100‑1 (MQ=255)

                                               tgCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCt                                                <  1:3510008/100‑1 (MQ=255)

                                                  ggCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGAACCTTgc                                             >  2:2572270/1‑99 (MQ=255)

                                                   gCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTa                                            <  2:3268841/100‑1 (MQ=255)

                                                   gCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTa                                            <  2:3395781/100‑1 (MQ=255)

                                                     atgatgTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGa                                          >  1:1603597/1‑100 (MQ=255)

                                                      tgatgTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGat                                         <  2:2659266/100‑1 (MQ=255)

                                                        atgTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGatat                                       <  1:2768682/100‑1 (MQ=255)

                                                         tgTCGGATGCATGAAGAATCGCGGCTTTATCCACCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATAtt                                      <  2:2142531/100‑1 (MQ=255)

                                                           gcGGATGCATGAAGAAACGCGGCTTTACCCAGCCCAAAATCTACAAAAGCCGCCTGCATCCCCGGCAGGACACGTCTTACACGACCCTTGTAGATATTGc                                    <  2:924881/99‑1 (MQ=255)

                                                             ggATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCct                                  >  2:1803477/1‑100 (MQ=255)

                                                             ggATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCct                                  <  1:1954056/100‑1 (MQ=255)

                                                             ggATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCCGCATACCCGGAAGTACACGACGTACACGACCCTTGTAGATATTGCct                                  <  1:3648338/100‑1 (MQ=255)

                                                               aTGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCctac                                <  1:953767/100‑1 (MQ=255)

                                                               aTGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCctac                                <  1:1230216/100‑1 (MQ=255)

                                                                tGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCctact                               >  1:1051949/1‑100 (MQ=255)

                                                                 gCATGAAGAAACGCGGCTTTATCCAGCCCAAGCTCTACAAAAGCCGCCTGCATACACGGAAGTACACGACTTAGACGAACCTTGTAGATATTGCctacta                              >  2:2081505/1‑100 (MQ=255)

                                                                 gCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGCTATTGCctacta                              >  2:756/1‑100 (MQ=255)

                                                                 gCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCctacta                              >  1:2461320/1‑100 (MQ=255)

                                                                  cATGGAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTAt                             >  1:3958428/1‑100 (MQ=255)

                                                                  cATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTAt                             >  2:1810269/1‑100 (MQ=255)

                                                                  cATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTAt                             >  2:179355/1‑100 (MQ=255)

                                                                    tGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTc                           <  1:1784183/100‑1 (MQ=255)

                                                                     gaagaaACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTcn                          <  1:1317526/100‑2 (MQ=255)

                                                                     gaagaaACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTcc                          >  2:3075608/1‑100 (MQ=255)

                                                                      aagaaACGCGGCTTTATCCAGCCCAAGATCTACAAATTCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCcg                         <  2:3875074/100‑1 (MQ=255)

                                                                       agaaACGCGGTTTTATCCAGCCCAAGACCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCcgc                        <  2:2846978/100‑1 (MQ=255)

                                                                          aaCGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTc                     >  1:480200/1‑100 (MQ=255)

                                                                               ggCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTCGATATTGCCTACTATTCCGCGTCGCGcc                <  1:3212261/100‑1 (MQ=255)

                                                                               ggCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGcc                <  1:3312772/100‑1 (MQ=255)

                                                                               ggCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGcc                <  1:3211317/100‑1 (MQ=255)

                                                                                gCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCt               <  1:1735512/100‑1 (MQ=255)

                                                                                  tttATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCa             <  1:2604903/100‑1 (MQ=255)

                                                                                   ttATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCAc            <  2:3180348/100‑1 (MQ=255)

                                                                                   ttATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCAc            >  2:2374817/1‑100 (MQ=255)

                                                                                    tATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACg           >  2:538202/1‑100 (MQ=255)

                                                                                      tCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCCCGtt         >  1:1381791/1‑100 (MQ=255)

                                                                                      tCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGtt         >  1:2961/1‑100 (MQ=255)

                                                                                       ccAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTc        >  2:1105352/1‑100 (MQ=255)

                                                                                       ccAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTc        <  2:2130141/100‑1 (MQ=255)

                                                                                        cAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCa       >  1:3859460/1‑100 (MQ=255)

                                                                                          gCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCaan     <  1:3239710/100‑2 (MQ=255)

                                                                                           cccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAata    <  2:1218388/100‑1 (MQ=255)

                                                                                           cccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAata    >  2:867797/1‑100 (MQ=255)

                                                                                           cccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAata    <  2:2459937/100‑1 (MQ=255)

                                                                                           cccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAata    <  1:2145032/100‑1 (MQ=255)

                                                                                            ccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAatat   >  2:460649/1‑100 (MQ=255)

                                                                                            ccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAatat   >  2:382913/1‑100 (MQ=255)

                                                                                            ccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAatat   >  2:903414/1‑100 (MQ=255)

                                                                                             cAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAATATg  >  1:551021/1‑100 (MQ=255)

                                                                                             cAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAATATg  >  1:550356/1‑100 (MQ=255)

                                                                                             cAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAATATg  >  2:2778006/1‑100 (MQ=255)

CGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAATATG  >  NC_000913/3397522‑3397714

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