BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	186,043	T→C	N870D (AAC→GAC)	glnD ←	protein PII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	186,043	0	T	C	99.6%	735.4 / ‑5.2	225	N870D (AAC→GAC)	glnD	protein PII
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (120/104); total (121/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 237 total aligned reads displayed.

AACTCTTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTTATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCGAAA  >  NC_000913/185953‑186132

aaCTCTTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGtcat                                                                                         >  2:896812/1‑93 (MQ=255)

aaCTCTTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGtcat                                                                                         <  1:3874462/93‑1 (MQ=255)

  ctctTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGtcatta                                                                                       >  1:1450500/1‑93 (MQ=255)

   tctTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTaa                                                                                      >  2:6532923/1‑93 (MQ=255)

   tctTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTaa                                                                                      <  2:1049677/93‑1 (MQ=255)

   tctTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTaa                                                                                      >  1:1978134/1‑93 (MQ=255)

    ctTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAg                                                                                     >  1:3164519/1‑93 (MQ=255)

     ttttCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAgc                                                                                    >  1:66635/1‑93 (MQ=255)

     ttttCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAgc                                                                                    >  2:827122/1‑93 (MQ=255)

     ttttCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAgc                                                                                    >  2:6285781/1‑93 (MQ=255)

     ttttCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAgc                                                                                    <  1:266944/93‑1 (MQ=255)

      tttCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAgcg                                                                                   <  1:1019317/93‑1 (MQ=255)

      tttCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAgcg                                                                                   >  2:4953233/1‑93 (MQ=255)

        tCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCa                                                                                 <  2:393487/93‑1 (MQ=255)

        tCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCa                                                                                 <  2:3969997/93‑1 (MQ=255)

         cATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCAc                                                                                <  1:4658924/93‑1 (MQ=255)

         cATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCAc                                                                                >  1:4185998/1‑93 (MQ=255)

          atatCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACg                                                                               >  2:4059626/1‑93 (MQ=255)

           tatCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGc                                                                              >  1:2861084/1‑93 (MQ=255)

           tatCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGc                                                                              <  1:3008675/93‑1 (MQ=255)

           tatCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGc                                                                              >  2:591207/1‑93 (MQ=255)

            atCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGcc                                                                             <  2:1888236/93‑1 (MQ=255)

               aGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTTATTAAGCGCACGCCGGt                                                                          >  1:5942403/1‑93 (MQ=255)

               aGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGt                                                                          >  2:2848242/1‑93 (MQ=255)

               aGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGt                                                                          >  2:5120261/1‑93 (MQ=255)

               aGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGt                                                                          >  1:1645255/1‑93 (MQ=255)

                 tAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCa                                                                        <  1:6700641/93‑1 (MQ=255)

                   aaCACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGc                                                                      <  2:3191710/93‑1 (MQ=255)

                    acacATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCCCGCCGGTCAGCg                                                                     >  1:1805629/1‑93 (MQ=255)

                    acacATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCg                                                                     >  1:2269805/1‑93 (MQ=255)

                      acaTCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGt                                                                   >  2:1542691/1‑93 (MQ=255)

                        aTCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTgg                                                                 >  2:2656759/1‑93 (MQ=255)

                         tCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGc                                                                <  2:2530669/93‑1 (MQ=255)

                          cACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCa                                                               >  1:3391252/1‑93 (MQ=255)

                           aCCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCaa                                                              >  1:2697660/1‑93 (MQ=255)

                           aCCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCaa                                                              >  2:3106573/1‑93 (MQ=255)

                           aCCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCaa                                                              >  1:2493169/1‑93 (MQ=255)

                            cccTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAAt                                                             <  2:2143774/93‑1 (MQ=255)

                              cTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTa                                                           <  2:6079312/93‑1 (MQ=255)

                              cTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTa                                                           <  1:747757/93‑1 (MQ=255)

                               tttATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTAt                                                          <  1:1424945/93‑1 (MQ=255)

                                ttATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATg                                                         >  2:1420597/1‑93 (MQ=255)

                                 tATCGTTTGGATTGAGGGCGTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGa                                                        <  1:5391224/93‑1 (MQ=255)

                                 tATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGa                                                        <  1:2938151/93‑1 (MQ=255)

                                  aTCGTTTGGGTTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGaa                                                       >  2:1339030/1‑93 (MQ=255)

                                   tCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAAt                                                      >  1:2748449/1‑93 (MQ=255)

                                    cGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATa                                                     >  1:5801817/1‑93 (MQ=255)

                                    cGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATa                                                     >  1:6070393/1‑93 (MQ=255)

                                    cGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATa                                                     >  1:379888/1‑93 (MQ=255)

                                    cGTTTGGATTGAGGGCATCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATa                                                     <  1:603853/93‑1 (MQ=255)

                                      tttGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATaaa                                                   >  1:5047057/1‑93 (MQ=255)

                                      tttGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATaaa                                                   <  2:4327735/93‑1 (MQ=255)

                                       ttGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAAt                                                  <  1:788058/93‑1 (MQ=255)

                                        tGGATTGATGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAAtc                                                 <  1:5583924/93‑1 (MQ=255)

                                        tGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAAtc                                                 <  1:4272336/93‑1 (MQ=255)

                                         ggATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAAtct                                                <  1:1559122/93‑1 (MQ=255)

                                         ggATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAAtct                                                <  1:8687/93‑1 (MQ=255)

                                         ggATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAAtct                                                >  1:2605383/1‑93 (MQ=255)

                                         ggATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAAtct                                                <  1:6601031/93‑1 (MQ=255)

                                           aTTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAAtcttc                                              <  2:4845180/93‑1 (MQ=255)

                                             tGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTa                                            >  1:4536840/1‑93 (MQ=255)

                                              gAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTAc                                           >  2:4285205/1‑93 (MQ=255)

                                               aGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTAct                                          <  2:5588230/93‑1 (MQ=255)

                                                gggCCTCTGTCAACCGCTGATGCGCTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTctcctc                                         <  2:1169029/93‑1 (MQ=255)

                                                gggCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTActc                                         >  2:684729/1‑93 (MQ=255)

                                                gggCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTActc                                         >  2:6391430/1‑93 (MQ=255)

                                                gggCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTActc                                         >  2:4758894/1‑93 (MQ=255)

                                                 ggCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTActcg                                        >  1:4327465/1‑93 (MQ=255)

                                                    ctctGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTActcgctc                                     >  2:2411707/1‑93 (MQ=255)

                                                     tctGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTActcgctcg                                    >  1:4973197/1‑93 (MQ=255)

                                                        gTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCa                                 <  2:1241431/93‑1 (MQ=255)

                                                        gTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCa                                 <  2:5246714/93‑1 (MQ=255)

                                                         tCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCaa                                <  2:6635350/93‑1 (MQ=255)

                                                           aaCCGCTGATGCACTTCCTGCTGCAACTCGTCATTACGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATg                              <  1:2205705/93‑1 (MQ=255)

                                                             ccGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGt                            <  2:5429390/93‑1 (MQ=255)

                                                              cGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCTCTCGCCAATGGtt                           >  1:4993699/1‑93 (MQ=255)

                                                              cGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGtt                           <  2:5655781/93‑1 (MQ=255)

                                                              cGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGtt                           <  1:6204224/93‑1 (MQ=255)

                                                               gCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTg                          >  2:3723230/1‑93 (MQ=255)

                                                                   aTGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAAt                      >  1:1320429/1‑93 (MQ=255)

                                                                     gCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTc                    <  1:3544433/93‑1 (MQ=255)

                                                                      cACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCg                   <  1:1151773/93‑1 (MQ=255)

                                                                       aCTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCgg                  >  1:6391249/1‑93 (MQ=255)

                                                                       aCTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCgg                  >  1:2049234/1‑93 (MQ=255)

                                                                        cTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCggg                 <  2:303277/93‑1 (MQ=255)

                                                                         ttCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGc                >  1:6488025/1‑93 (MQ=255)

                                                                         ttCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGc                >  1:1002040/1‑93 (MQ=255)

                                                                          tCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCa               >  1:5458735/1‑93 (MQ=255)

                                                                          tCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCa               >  2:4179201/1‑93 (MQ=255)

                                                                           ccTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCAc              >  1:4409383/1‑93 (MQ=255)

                                                                             tgctgcAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCa            <  2:3304282/93‑1 (MQ=255)

                                                                              gctgcAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCAt           <  2:4644152/93‑1 (MQ=255)

                                                                              gctgcAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCAt           <  1:4745784/93‑1 (MQ=255)

                                                                               ctgcAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATg          >  2:3530105/1‑93 (MQ=255)

                                                                                 gcAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGaa        >  2:4440342/1‑93 (MQ=255)

                                                                                 gcAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGaa        >  1:5663950/1‑93 (MQ=255)

                                                                                   aaCTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGc      <  1:5917795/93‑1 (MQ=255)

                                                                                    aCTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCg     <  2:5849370/93‑1 (MQ=255)

                                                                                    aCTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCg     >  1:225525/1‑93 (MQ=255)

                                                                                       cGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGCTGTAATTCGGGCACCATGAAGCGaaa  >  2:1716482/1‑93 (MQ=255)

AACTCTTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTTATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCGAAA  >  NC_000913/185953‑186132

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