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	95.9%	436.4 / 6.2	145	N870D (AAC→GAC)	glnD	protein PII
Reads supporting (aligned to +/- strand): ref base T (4/2); new base C (76/63); total (80/65)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 158 total aligned reads displayed.

TAAACTCTTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTTATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCGAAA  >  NC_000913/185951‑186132

tAAACTCTTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGtt                                                                                           >  2:3628297/1‑93 (MQ=255)

tAAACTCTTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGtt                                                                                           >  2:3628296/1‑93 (MQ=255)

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

   aCTCTTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGtcatt                                                                                        >  1:322562/1‑93 (MQ=255)

    ctctTTTCATATGAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGtcatta                                                                                       >  1:982379/1‑93 (MQ=255)

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

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

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

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

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

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

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

         ttCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAgcgc                                                                                  <  1:358922/93‑1 (MQ=255)

         ttCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAgcgc                                                                                  >  2:4579124/1‑93 (MQ=255)

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

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

               tCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCg                                                                            <  2:1191512/93‑1 (MQ=255)

               tCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCg                                                                            >  2:3502359/1‑93 (MQ=255)

               tCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCg                                                                            <  2:741608/93‑1 (MQ=255)

               tCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCg                                                                            >  1:3266547/1‑93 (MQ=255)

                 aGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGt                                                                          <  1:2662498/93‑1 (MQ=255)

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

                  gTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTc                                                                         >  1:3079297/1‑93 (MQ=255)

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

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

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

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

                       cacaTCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCgg                                                                    >  1:237643/1‑93 (MQ=255)

                       cacaTCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCgg                                                                    >  2:1416829/1‑93 (MQ=255)

                        acaTCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGt                                                                   <  2:3110055/93‑1 (MQ=255)

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

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

                         caTCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTTATTAAGCGCACGCCGGTCAGCGGTg                                                                  >  1:3896009/1‑93 (MQ=255)

                         caTCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTg                                                                  <  1:2206368/93‑1 (MQ=255)

                         caTCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTg                                                                  >  2:3774610/1‑93 (MQ=255)

                           tCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGc                                                                >  2:485795/1‑93 (MQ=255)

                           tCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGc                                                                >  1:3947688/1‑93 (MQ=255)

                           tCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGc                                                                >  2:155243/1‑93 (MQ=255)

                            cACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCa                                                               <  2:4228061/93‑1 (MQ=255)

                            cACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTCGCa                                                               >  1:3480829/1‑93 (MQ=255)

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

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

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

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

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

                               ccTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAAtt                                                            <  2:1724726/93‑1 (MQ=255)

                               ccTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAAtt                                                            <  2:2340368/93‑1 (MQ=255)

                               ccTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAAtt                                                            <  1:1462066/93‑1 (MQ=255)

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

                                cTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTa                                                           >  1:1812644/1‑93 (MQ=255)

                                 tttATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTAt                                                          >  2:1374271/1‑93 (MQ=255)

                                  ttATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTTATTAAGCGCACGCCGGTCAGCGGTGGCAATTATg                                                         >  1:3091817/1‑93 (MQ=255)

                                   tATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGa                                                        >  1:1783075/1‑93 (MQ=255)

                                   tATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGa                                                        >  1:458461/1‑93 (MQ=255)

                                    aTCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAAATCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGaa                                                       >  1:250010/1‑93 (MQ=255)

                                     tCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATTAAt                                                      <  2:2772043/93‑1 (MQ=255)

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

                                       gTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATaa                                                    >  1:647162/1‑93 (MQ=255)

                                       gTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATaa                                                    <  2:529721/93‑1 (MQ=255)

                                        tttGGATTGAGGGCCTCTGTCACACGCTTATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATaaa                                                   <  1:1606664/93‑1 (MQ=255)

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

                                         ttGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAAATCGTTATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAAt                                                  >  1:3917443/1‑93 (MQ=255)

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

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

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

                                               tGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTTTa                                            >  1:1161405/1‑93 (MQ=255)

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

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

                                                     ccTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTActcgct                                      >  2:2122778/1‑93 (MQ=255)

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

                                                         tGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGcc                                  >  2:159620/1‑93 (MQ=255)

                                                         tGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGcc                                  >  2:1191799/1‑93 (MQ=255)

                                                          gTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCa                                 >  2:212724/1‑93 (MQ=255)

                                                          gTCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCa                                 >  2:487250/1‑93 (MQ=255)

                                                           tCAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCat                                >  2:4715116/1‑92 (MQ=255)

                                                            cAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTTCTCGCTCGCCAAt                               <  2:1285846/93‑1 (MQ=255)

                                                            cAACCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAAt                               >  1:283613/1‑93 (MQ=255)

                                                             aaCCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATg                              >  2:4702402/1‑93 (MQ=255)

                                                              aCCGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATgg                             <  1:546265/93‑1 (MQ=255)

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

                                                               ccGCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGt                            >  2:1132386/1‑93 (MQ=255)

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

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

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

                                                                 gCTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTg                          <  2:2227298/93‑1 (MQ=255)

                                                                  cTGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGt                         <  1:2766144/93‑1 (MQ=255)

                                                                   tGATGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTa                        <  2:584890/93‑1 (MQ=255)

                                                                      tGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAAtt                     >  1:52317/1‑93 (MQ=255)

                                                                      tGCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAAtt                     <  2:3090479/93‑1 (MQ=255)

                                                                       gCACTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTc                    >  2:111673/1‑93 (MQ=255)

                                                                          cTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCggg                 >  1:3464605/1‑93 (MQ=255)

                                                                          cTTCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCggg                 >  2:749429/1‑93 (MQ=255)

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

                                                                           ttCCTGCTGCAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGCTGTAATTCTGGc                >  2:388229/1‑93 (MQ=255)

                                                                                  tgcAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGCATACATCTTCTACTCGCTCGCCAATGGTTGTCCTTCGGGCACAATGa         >  2:1144569/1‑93 (MQ=255)

                                                                                  tgcAACTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGa         <  2:1210415/93‑1 (MQ=255)

                                                                                       cTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCGa    >  1:812099/1‑93 (MQ=255)

                                                                                       cTCGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCGa    >  2:3264237/1‑93 (MQ=255)

                                                                                         cGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCGaaa  <  2:1871500/93‑1 (MQ=255)

                                                                                         cGTCATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCGaaa  <  1:2011401/93‑1 (MQ=255)

TAAACTCTTTTCATATCAGTAAACACATCACCCTTTATCGTTTGGATTGAGGGCCTCTGTCAACCGCTGATGCACTTCCTGCTGCAACTCGTTATTAAGCGCACGCCGGTCAGCGGTGGCAATTATGAATAAATCTTCTACTCGCTCGCCAATGGTTGTAATTCGGGCACCATGAAGCGAAA  >  NC_000913/185951‑186132

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