BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,689	C→T	D135D (GAC→GAT)	ytfI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,689	0	C	T	99.5%	723.3 / ‑5.4	197	D135D (GAC→GAT)	ytfI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (105/91); total (105/92)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

Only 100 of 201 total aligned reads displayed.

TTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑AAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTG  >  NC_000913/4435601‑4435778

ttGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAgatgat                                                                                        <  1:382109/93‑1 (MQ=255)

  gAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTg                                                                                      <  1:3967912/93‑1 (MQ=255)

   aaTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGt                                                                                     >  1:3900141/1‑93 (MQ=255)

    aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTa                                                                                    <  2:812363/93‑1 (MQ=255)

    aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTa                                                                                    <  2:376413/93‑1 (MQ=255)

     tACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTAt                                                                                   >  2:2353866/1‑93 (MQ=255)

     aaCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTAt                                                                                   <  2:1733602/92‑1 (MQ=255)

      aCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATg                                                                                  >  1:3594038/1‑93 (MQ=255)

      aCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTGCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATg                                                                                  <  1:2014035/93‑1 (MQ=255)

          agaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTg                                                                              <  1:389049/93‑1 (MQ=255)

          agaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTg                                                                              <  1:1285825/93‑1 (MQ=255)

           gaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTgg                                                                             >  2:4318708/1‑93 (MQ=255)

            aaaTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTggg                                                                            >  2:1110019/1‑93 (MQ=255)

              aTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGat                                                                          <  1:1731885/93‑1 (MQ=255)

               tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGtta                                                                         <  2:2935695/93‑1 (MQ=255)

                aCCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGatat                                                                        >  1:2361379/1‑93 (MQ=255)

                 cccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATAtt                                                                       >  1:3094186/1‑93 (MQ=255)

                 cccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATAtt                                                                       >  1:1505065/1‑93 (MQ=255)

                 cccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATAtt                                                                       <  2:1493561/93‑1 (MQ=255)

                  ccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTc                                                                      >  1:251035/1‑93 (MQ=255)

                  ccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTc                                                                      >  2:4289001/1‑93 (MQ=255)

                   cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATAttat                                                                     >  1:3697354/1‑91 (MQ=255)

                   cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCt                                                                     <  2:656858/93‑1 (MQ=255)

                     tAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTAc                                                                   >  1:2681605/1‑93 (MQ=255)

                     tAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTAc                                                                   >  2:247388/1‑93 (MQ=255)

                      aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACa                                                                  <  2:4216899/93‑1 (MQ=255)

                       aCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAg                                                                 >  1:1500340/1‑93 (MQ=255)

                           tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaa                                                             >  2:1802319/1‑93 (MQ=255)

                           tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaa                                                             >  1:1027078/1‑93 (MQ=255)

                           tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaa                                                             >  1:3494979/1‑93 (MQ=255)

                           tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaa                                                             >  2:3202509/1‑93 (MQ=255)

                             aaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaaaa                                                           <  2:4615635/93‑1 (MQ=255)

                              aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                          >  1:439065/1‑93 (MQ=255)

                              aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                          <  1:2644880/93‑1 (MQ=255)

                               tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                                         >  2:2629792/1‑93 (MQ=255)

                               tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                                         <  1:2126373/93‑1 (MQ=255)

                                  atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                                      >  2:2678517/1‑93 (MQ=255)

                                  atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                                      <  2:1413008/93‑1 (MQ=255)

                                  atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                                      <  2:578112/93‑1 (MQ=255)

                                   tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                                     <  1:1889585/93‑1 (MQ=255)

                                   tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                                     <  1:804017/93‑1 (MQ=255)

                                   tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                                     >  1:1883720/1‑93 (MQ=255)

                                   tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                                     >  2:3295694/1‑93 (MQ=255)

                                     taACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                                                   >  1:2033687/1‑93 (MQ=255)

                                       aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                                 <  2:284264/93‑1 (MQ=255)

                                       aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                                 >  1:294536/1‑93 (MQ=255)

                                        cTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                                                >  1:1563096/1‑93 (MQ=255)

                                          cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAc                                              >  2:7917/1‑93 (MQ=255)

                                           attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAcc                                             <  2:4125525/93‑1 (MQ=255)

                                            ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                            >  1:521869/1‑93 (MQ=255)

                                            ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                            >  1:1313742/1‑93 (MQ=255)

                                            ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                            >  1:1358475/1‑93 (MQ=255)

                                             tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGa                                           >  2:1911394/1‑93 (MQ=255)

                                              catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAATCTAACATACCGaa                                          >  2:3055162/1‑93 (MQ=255)

                                               atGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                                         >  2:2788810/1‑93 (MQ=255)

                                               atGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                                         >  1:4062030/1‑93 (MQ=255)

                                                tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                                        <  2:4413295/93‑1 (MQ=255)

                                                  gATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg                                      >  1:4642839/1‑93 (MQ=255)

                                                   aTTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa                                     >  1:3203886/1‑93 (MQ=255)

                                                   aTTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa                                     >  2:1301530/1‑93 (MQ=255)

                                                    tttCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt                                    <  2:2247964/93‑1 (MQ=255)

                                                     ttCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt                                   <  2:1815564/93‑1 (MQ=255)

                                                     ttCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt                                   <  2:2064156/93‑1 (MQ=255)

                                                     ttCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt                                   <  1:112656/93‑1 (MQ=255)

                                                      tCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa                                  <  1:4257709/93‑1 (MQ=255)

                                                       cAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaa                                 <  2:510773/93‑1 (MQ=255)

                                                       cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaa                                 <  2:708729/93‑1 (MQ=255)

                                                        aaaaaCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa                                <  1:286816/93‑1 (MQ=255)

                                                         aaaaCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAg                               <  2:296939/93‑1 (MQ=255)

                                                         aaaaCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAg                               >  2:1131394/1‑93 (MQ=255)

                                                           aaCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                             >  2:3965290/1‑93 (MQ=255)

                                                            aCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGtt                            >  1:4468401/1‑93 (MQ=255)

                                                             ccGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTta                           <  1:3850917/93‑1 (MQ=255)

                                                             ccGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTta                           <  1:2958287/93‑1 (MQ=255)

                                                              cGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtat                          >  2:3486918/1‑93 (MQ=255)

                                                               gTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtata                         <  1:2578754/93‑1 (MQ=255)

                                                                tGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTatat                        >  1:3380917/1‑93 (MQ=255)

                                                                tGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTatat                        >  1:3230337/1‑93 (MQ=255)

                                                                  aaTTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTc                      <  1:2660420/93‑1 (MQ=255)

                                                                    tttCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAc                    >  2:2106225/1‑93 (MQ=255)

                                                                    tttCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAc                    >  2:4497207/1‑93 (MQ=255)

                                                                    tttCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAc                    <  1:1627014/93‑1 (MQ=255)

                                                                     ttCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAcc                   <  1:3708480/93‑1 (MQ=255)

                                                                        cGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAgg                <  2:4060682/93‑1 (MQ=255)

                                                                         gTTAGAACAAAACAGATGATTTTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGc               <  1:1932344/93‑1 (MQ=255)

                                                                         gTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGc               >  1:2773544/1‑93 (MQ=255)

                                                                         gTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGc               >  2:4373424/1‑93 (MQ=255)

                                                                          ttAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCg              >  1:3021248/1‑93 (MQ=255)

                                                                          ttAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCg              <  2:2288043/93‑1 (MQ=255)

                                                                            aGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCgct            >  2:2817197/1‑91 (MQ=255)

                                                                               acaaaacaGATGATTGTATGCCTGGGTTATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTgg         <  1:572835/90‑1 (MQ=255)

                                                                               acaaaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTgg         >  1:513890/4‑93 (MQ=255)

                                                                               acaaaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTgg         <  1:788375/90‑1 (MQ=255)

                                                                                caaaacaGATGATTGTATGCCTTGGATATTCTACAGAAAAAACAGTCTGACATCCCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGt        <  2:454034/91‑1 (MQ=255)

                                                                                 aaaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTg       >  1:393067/2‑93 (MQ=255)

                                                                                 aaaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTg       >  2:630450/2‑93 (MQ=255)

                                                                                  aaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGa      <  1:2437978/93‑1 (MQ=255)

                                                                                   aacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGAc     >  1:1283776/1‑93 (MQ=255)

                                                                                     caGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACgt   <  2:4034508/93‑1 (MQ=255)

                                                                                      aGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACgtg  >  1:1158603/1‑93 (MQ=255)

TTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑AAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTG  >  NC_000913/4435601‑4435778

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