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.5%	749.2 / NA	210	I60L (ATT→CTT)	rng	ribonuclease G
Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (108/100); total (109/101)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

TCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAATATG  >  NC_000913/3397521‑3397714

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

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

 cgcgCACCGTGAATTGCTTTTATTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATc                                                                                               >  2:2450610/1‑100 (MQ=255)

    gcACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTAc                                                                                            <  1:2811427/100‑1 (MQ=255)

     cACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACa                                                                                           >  1:2593885/1‑100 (MQ=255)

     cACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACa                                                                                           >  2:1530538/1‑100 (MQ=255)

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

       ccGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACaaa                                                                                         <  1:4595793/100‑1 (MQ=255)

         gTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAg                                                                                       >  2:2122483/1‑100 (MQ=255)

         gTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAg                                                                                       <  2:1346994/100‑1 (MQ=255)

            aaTTGCTTTTGTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAgccg                                                                                    <  1:5065669/100‑1 (MQ=255)

             aTTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAgccgc                                                                                   >  1:4318566/1‑100 (MQ=255)

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

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

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

               tGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCt                                                                                 >  1:1899877/1‑100 (MQ=255)

                gCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAGGCCGCCTg                                                                                <  2:3436433/100‑1 (MQ=255)

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

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

                gCTTTTGTTCTTCACCCGCCACACATTCGGTGAGCGGCCTGATGTCCGATGCATGAAGAAACGCGGCTTTATCCAGCCCAGGATCTACAAAAGCCGCCTg                                                                                >  1:4954905/1‑100 (MQ=255)

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

                   tttGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCAt                                                                             <  1:5464558/100‑1 (MQ=255)

                   tttGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCAt                                                                             <  1:2693227/100‑1 (MQ=255)

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

                      gTTCTTCACCCTCCACAGTTTCGGTGTGCGACATGATGTCAGATGCATCAAGACACGCCGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATAcc                                                                          <  2:788065/100‑1 (MQ=255)

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

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

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

                         cttcACCCGCCACACGTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGTTTTGTGCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCgg                                                                       <  2:491742/100‑1 (MQ=255)

                          ttcACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGa                                                                      <  2:4238596/100‑1 (MQ=255)

                          ttcACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGa                                                                      >  1:2750254/1‑100 (MQ=255)

                            cACCCGCCACACATTCGGTGTGCGGTATGATTTCGGATACATGTAGAATCAAGACTTTACCAAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAg                                                                    <  2:1301679/100‑1 (MQ=255)

                            cACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAg                                                                    <  2:1893392/100‑1 (MQ=255)

                            cACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAg                                                                    <  1:1070621/100‑1 (MQ=255)

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

                             aCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGt                                                                   <  2:3983378/100‑1 (MQ=255)

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

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

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

                               ccGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTac                                                                 <  2:3452146/100‑1 (MQ=255)

                               ccGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTac                                                                 <  1:3135847/100‑1 (MQ=255)

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

                                     cacaTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACt                                                           >  2:395144/1‑100 (MQ=255)

                                        aTTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACTCGGAAGTACACGACTTac                                                        <  1:403078/100‑1 (MQ=255)

                                         ttCGGTATGCGGCATGATGTCGGATGCATGAAGAACCGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTaca                                                       <  1:1836196/100‑1 (MQ=255)

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

                                             gtgtGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCCAGATCTACAAAAGCCGCCTGCATAGCCGGAAGTACACGACTTACACGAc                                                   >  2:3249921/1‑100 (MQ=255)

                                             gtgtGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGAc                                                   <  1:3930909/100‑1 (MQ=255)

                                               gtgCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGAATTACAAGAccc                                                 >  1:2258970/1‑100 (MQ=255)

                                                 gCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCtt                                               >  1:3357236/1‑100 (MQ=255)

                                                 gCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCAAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCtt                                               <  2:1713514/100‑1 (MQ=255)

                                                   ggCATGATGTCGGATGCATGAAGCAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGt                                             <  2:3901233/100‑1 (MQ=255)

                                                      atgatgTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATTCCCGGAAGTACACGACTTACACGGCCCTTGTAGa                                          <  2:5076572/100‑1 (MQ=255)

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

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

                                                      atgatgTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGa                                          <  1:607619/100‑1 (MQ=255)

                                                      atgatgTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTGGTAGa                                          >  2:2501773/1‑100 (MQ=255)

                                                       tgatgTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGat                                         >  2:3712648/1‑100 (MQ=255)

                                                         atgTCGGATGCATGAAGAAACGCGGCTTTATCCAGTCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTatatat                                       <  2:4647959/100‑1 (MQ=255)

                                                          tgTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATAtt                                      <  1:2949159/100‑1 (MQ=255)

                                                             cGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGcc                                   >  1:1423886/1‑100 (MQ=255)

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

                                                               gATGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCcta                                 >  2:1437997/1‑100 (MQ=255)

                                                                 tGCATGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCctact                               <  1:2548157/100‑1 (MQ=255)

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

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

                                                                     tGAAGAAACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTc                           >  2:3063649/1‑100 (MQ=255)

                                                                      gaagaaACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTTCCTACTATTcc                          >  2:113581/1‑100 (MQ=255)

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

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

                                                                       aagaaACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCcg                         >  2:2096105/1‑100 (MQ=255)

                                                                        agaaACGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCcgc                        >  2:3491842/1‑100 (MQ=255)

                                                                          aaaCGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGt                      >  2:3790252/1‑100 (MQ=255)

                                                                          aaaCGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGt                      >  1:3846945/1‑100 (MQ=255)

                                                                           aaCGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTc                     <  2:4823407/100‑1 (MQ=255)

                                                                           aaCGCGGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTc                     <  1:2793045/100‑1 (MQ=255)

                                                                              gcgGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTcgcg                  >  1:1397438/1‑100 (MQ=255)

                                                                              gcgGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTcgcg                  <  1:1119637/100‑1 (MQ=255)

                                                                               cgGCTTTATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCgcgc                 >  2:1269820/1‑100 (MQ=255)

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

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

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

                                                                                    gatctcAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCAc            <  2:3135736/95‑1 (MQ=255)

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

                                                                                     tATCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACg           <  1:3780600/100‑1 (MQ=255)

                                                                                      aTCCAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGt          >  2:1966736/1‑100 (MQ=255)

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

                                                                                         cAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCAcagatc       >  1:3165025/1‑95 (MQ=255)

                                                                                         cAGCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCa       <  1:2284591/100‑1 (MQ=255)

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

                                                                                           gCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAat     <  1:2530687/100‑1 (MQ=255)

                                                                                           gCCCAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAat     <  2:2960972/100‑1 (MQ=255)

                                                                                            cccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCaaag    >  1:1461256/1‑98 (MQ=255)

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

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

                                                                                             ccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCCCGTTCCatat   >  2:3180897/1‑100 (MQ=255)

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

                                                                                             ccAAGATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAatat   <  2:4064338/100‑1 (MQ=255)

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

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

TCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTACTATTCCGCGTCGCGCCTCACGTTCAATATG  >  NC_000913/3397521‑3397714

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