BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,073,511	A→G	L231P (CTG→CCG)	yggF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,073,511	0	A	G	99.3%	959.8 / ‑4.9	300	L231P (CTG→CCG)	yggF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (166/132); total (167/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

Only 100 of 314 total aligned reads displayed.

TCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCG  >  NC_000913/3073421‑3073600

tCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGc                                                                                         >  2:2275804/1‑93 (MQ=255)

gcAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGc                                                                                         <  2:709382/92‑1 (MQ=255)

 cAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCt                                                                                        >  2:3793379/1‑93 (MQ=255)

  aaCGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTc                                                                                       >  1:4766717/1‑93 (MQ=255)

   aCGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCg                                                                                      >  2:5472870/1‑93 (MQ=255)

     gCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAAtcggctcggc                                                                                    >  2:4752472/1‑93 (MQ=255)

       ccATTGCCTTGCACCGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTATCATGGCAAAAGTCAATCGGCTCGGCCt                                                                                  <  1:4543430/93‑1 (MQ=255)

       ccATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCt                                                                                  <  2:1099074/93‑1 (MQ=255)

        cATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTg                                                                                 >  2:1186312/1‑93 (MQ=255)

        cATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTg                                                                                 >  2:1460255/1‑93 (MQ=255)

        cATTGCCTTGCAACGCGTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGACAATCGGCTCGGCCTg                                                                                 <  1:2695710/93‑1 (MQ=255)

          ttGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCa                                                                               <  1:1616697/93‑1 (MQ=255)

            gCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCata                                                                             <  2:3305371/93‑1 (MQ=255)

             ccTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCatat                                                                            <  2:1240023/93‑1 (MQ=255)

              cTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCTTGGCAAAAGTCAATCGGCTCGGCCTGCATATc                                                                           <  1:3203599/93‑1 (MQ=255)

              cTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATc                                                                           <  1:1042192/93‑1 (MQ=255)

                tGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATccc                                                                         >  2:2986357/1‑93 (MQ=255)

                 tcAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGCCAAAAGTCAATCGGCTCGGCCTGCATATcccc                                                                        <  1:2256375/92‑1 (MQ=255)

                  cAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCg                                                                       <  2:5011794/93‑1 (MQ=255)

                    aCGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGcc                                                                     >  1:2975172/1‑93 (MQ=255)

                    aCGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGcc                                                                     <  1:3436653/93‑1 (MQ=255)

                    aCGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGcc                                                                     >  2:2516730/1‑93 (MQ=255)

                      gCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTa                                                                   >  2:6390342/1‑93 (MQ=255)

                       cTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTaa                                                                  <  1:606235/93‑1 (MQ=255)

                        ttGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAAt                                                                 <  2:6252257/93‑1 (MQ=255)

                            gcTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCtt                                                             <  1:2308517/93‑1 (MQ=255)

                            gcTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCtt                                                             >  1:5632750/1‑93 (MQ=255)

                             cTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCttt                                                            >  1:3586527/1‑93 (MQ=255)

                             cTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCttt                                                            >  1:747789/1‑93 (MQ=255)

                             cTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCttt                                                            >  1:1722817/1‑93 (MQ=255)

                              tCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTa                                                           >  2:1996475/1‑93 (MQ=255)

                              tCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTa                                                           >  2:730295/1‑93 (MQ=255)

                              tCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTa                                                           >  1:1421986/1‑93 (MQ=255)

                               ccTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTaa                                                          <  2:374076/93‑1 (MQ=255)

                                cTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAAc                                                         >  2:3405017/1‑93 (MQ=255)

                                cTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAAc                                                         >  1:4331852/1‑93 (MQ=255)

                                cTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAAc                                                         >  2:3325039/1‑93 (MQ=255)

                                cTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAAc                                                         <  1:501627/93‑1 (MQ=255)

                                 tGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACg                                                        >  2:6005125/1‑93 (MQ=255)

                                  gCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACgg                                                       >  2:5438273/1‑93 (MQ=255)

                                    tctGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGca                                                     <  1:5252968/93‑1 (MQ=255)

                                    tctGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGca                                                     >  2:3197657/1‑93 (MQ=255)

                                      tGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGcaca                                                   <  1:1096556/93‑1 (MQ=255)

                                       gCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAg                                                  >  1:5878785/1‑93 (MQ=255)

                                       gCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAg                                                  <  1:4169760/93‑1 (MQ=255)

                                        cAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTTATTCTTTAACGGCACAgg                                                 >  2:4044127/1‑93 (MQ=255)

                                         aaTCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTTACGGCACAGGc                                                >  2:6067040/1‑93 (MQ=255)

                                         aaTCTGCCGATTTTCCGTGTAATCCCCTTTACCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGc                                                <  1:5131565/93‑1 (MQ=255)

                                          aTCTTCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCa                                               <  2:2943728/93‑1 (MQ=255)

                                          aTCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATTCCCGCCTAATGCTTTTACGGCACAGGCa                                               >  1:3322388/1‑93 (MQ=255)

                                           tCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAg                                              >  1:2804532/1‑93 (MQ=255)

                                           tCTGCCGATTTTCCGTGTAATCACCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATCTCCCCGCCTAATGATTTAACGGCACAGGCAg                                              >  1:2217248/1‑93 (MQ=255)

                                             tGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAGTCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaa                                            >  2:4700815/1‑93 (MQ=255)

                                             tGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaa                                            >  2:3845479/1‑93 (MQ=255)

                                             tGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaa                                            >  1:2473708/1‑93 (MQ=255)

                                             tGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaa                                            >  1:2356088/1‑93 (MQ=255)

                                              gCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaaa                                           >  2:6701989/1‑93 (MQ=255)

                                                cGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATc                                         <  2:4446236/93‑1 (MQ=255)

                                                 gATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCa                                        >  1:4892073/1‑93 (MQ=255)

                                                  aTTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCAc                                       >  2:1593037/1‑93 (MQ=255)

                                                   ttttCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACg                                      >  2:189133/1‑93 (MQ=255)

                                                    tttCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGc                                     <  2:2653538/93‑1 (MQ=255)

                                                      tCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCt                                   <  1:4013667/93‑1 (MQ=255)

                                                      tCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCt                                   <  1:6317648/93‑1 (MQ=255)

                                                        cGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTc                                 <  2:5769841/93‑1 (MQ=255)

                                                        cGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTc                                 <  2:6297379/93‑1 (MQ=255)

                                                        cGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTc                                 <  2:3783588/93‑1 (MQ=255)

                                                         gtgtAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCt                                <  2:6488524/93‑1 (MQ=255)

                                                          tgtAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTg                               >  1:875285/1‑93 (MQ=255)

                                                          tgtAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTg                               <  1:2645640/93‑1 (MQ=255)

                                                            tAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATTCCCGCCTAATGCTTTAACGGCACAGGCAGCAATCACGCCTTCTGgc                             >  1:3694855/1‑93 (MQ=255)

                                                              aTCCCCTTTGTCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGgcgc                           <  2:1801414/93‑1 (MQ=255)

                                                              aTCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGgcgc                           >  1:3819696/1‑93 (MQ=255)

                                                                ccccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAc                         <  2:6439577/93‑1 (MQ=255)

                                                                ccccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCAGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCAGGCGCAc                         >  2:6614541/1‑93 (MQ=255)

                                                                 cccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAcc                        >  2:2815596/1‑93 (MQ=255)

                                                                 cccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAcc                        >  1:5051275/1‑93 (MQ=255)

                                                                 cccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAcc                        <  2:1939464/93‑1 (MQ=255)

                                                                 cccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAcc                        <  1:6592912/93‑1 (MQ=255)

                                                                    tttACCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgcc                     <  2:432751/93‑1 (MQ=255)

                                                                     ttAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgccg                    >  2:175296/1‑93 (MQ=255)

                                                                     ttAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgccg                    >  1:3902213/1‑93 (MQ=255)

                                                                      tAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGa                   >  1:1842683/1‑93 (MQ=255)

                                                                       aGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGAt                  <  1:5534454/93‑1 (MQ=255)

                                                                        gCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATg                 >  1:1204308/1‑93 (MQ=255)

                                                                         ccTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATgg                <  2:4915341/93‑1 (MQ=255)

                                                                          cTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgt               >  1:6620104/1‑93 (MQ=255)

                                                                          cTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgt               >  1:2054578/1‑93 (MQ=255)

                                                                           tGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtg              >  1:4044249/1‑93 (MQ=255)

                                                                            ggCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCGGAAATCACGCCTTCTGGCGCACCGCCGATGgtgt             <  1:950392/93‑1 (MQ=255)

                                                                            ggCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtgt             <  1:6397877/93‑1 (MQ=255)

                                                                            ggCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtga             >  2:517224/1‑92 (MQ=255)

                                                                             gCAAAAGTCAATCGGCTCGGCCTGCATATCCGCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTa            <  1:4798680/93‑1 (MQ=255)

                                                                               aaaaGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTaca          <  2:2338313/93‑1 (MQ=255)

                                                                                 aaGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatc        <  1:5906514/93‑1 (MQ=255)

                                                                                    tCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAATTCACGCCTTCTGGCGCACCGCCGATGGTGTacatcaca     <  2:1720141/93‑1 (MQ=255)

                                                                                     cAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcacat    <  2:1343358/93‑1 (MQ=255)

                                                                                     cAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcacat    >  2:498861/1‑93 (MQ=255)

                                                                                      aaTCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcacatc   >  1:6669352/1‑93 (MQ=255)

                                                                                       aTCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCg  <  2:2037557/93‑1 (MQ=255)

TCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCG  >  NC_000913/3073421‑3073600

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