BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,562,351	G→A	G118S (GGC→AGC)	b2449 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,562,351	0	G	A	99.6%	988.9 / ‑5.9	270	G118S (GGC→AGC)	b2449	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (126/143); total (127/143)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 271 total aligned reads displayed.

ATCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGG  >  NC_000913/2562261‑2562440

aTCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGc                                                                                         >  2:85650/1‑93 (MQ=255)

  ccTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCtt                                                                                       <  1:7803591/93‑1 (MQ=255)

  ccTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCtt                                                                                       >  1:3903373/1‑93 (MQ=255)

  ccTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCtt                                                                                       <  1:3168247/93‑1 (MQ=255)

   cTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCttt                                                                                      >  1:1621774/1‑93 (MQ=255)

    tGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTg                                                                                     <  1:4624822/93‑1 (MQ=255)

    tGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTg                                                                                     <  2:2354582/93‑1 (MQ=255)

      gTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGtt                                                                                   <  1:6690773/93‑1 (MQ=255)

        aaaaTACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGtttt                                                                                 >  1:4443505/1‑93 (MQ=255)

           aTACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggt                                                                              >  2:6582142/1‑93 (MQ=255)

            tACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggta                                                                             >  1:3868708/1‑93 (MQ=255)

             aCCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggtag                                                                            <  2:4739259/93‑1 (MQ=255)

              ccATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggtagg                                                                           >  2:145458/1‑93 (MQ=255)

                aTTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggtaggta                                                                         <  1:2306384/93‑1 (MQ=255)

                aTTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggtaggta                                                                         <  1:5637900/93‑1 (MQ=255)

                aTTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggtaggta                                                                         <  2:1493278/93‑1 (MQ=255)

                aTTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggtaggta                                                                         >  1:4160324/1‑93 (MQ=255)

                 tttCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGgtaggtag                                                                        <  2:2992305/93‑1 (MQ=255)

                 tttCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGgtaggtag                                                                        >  1:1737893/1‑93 (MQ=255)

                  ttCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTaga                                                                       <  1:7169996/93‑1 (MQ=255)

                   tCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTagag                                                                      <  2:2587689/93‑1 (MQ=255)

                   tCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTagag                                                                      <  2:3347003/93‑1 (MQ=255)

                   tCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTagag                                                                      <  2:1142541/93‑1 (MQ=255)

                    cAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGt                                                                     >  1:263500/1‑93 (MQ=255)

                      gCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTAt                                                                   <  1:301822/93‑1 (MQ=255)

                      gCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTAt                                                                   >  1:5170970/1‑93 (MQ=255)

                       cAATGAAACTGAGTTCAAACTTAAAGCGGGAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTAtt                                                                  <  1:2003989/93‑1 (MQ=255)

                        aaTGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTAttt                                                                 >  2:3070434/1‑93 (MQ=255)

                        aaTGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTAttt                                                                 >  2:1051755/1‑93 (MQ=255)

                        aaTGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTAttt                                                                 <  2:4096030/93‑1 (MQ=255)

                          tGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGa                                                               <  2:6615781/93‑1 (MQ=255)

                          tGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGa                                                               >  2:5195950/1‑93 (MQ=255)

                             aaCTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGaaaa                                                            >  2:4990035/1‑93 (MQ=255)

                             aaCTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGaaaa                                                            >  1:1171954/1‑93 (MQ=255)

                              aCTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGaaaaa                                                           <  2:313278/93‑1 (MQ=255)

                                 gAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGt                                                        <  2:7382088/93‑1 (MQ=255)

                                 gAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGt                                                        >  1:5434786/1‑93 (MQ=255)

                                    ttCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAat                                                     <  2:7647774/93‑1 (MQ=255)

                                    ttCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAat                                                     >  2:6431140/1‑93 (MQ=255)

                                    ttCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAat                                                     >  1:827989/1‑93 (MQ=255)

                                       aaaCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattga                                                  <  2:3449498/93‑1 (MQ=255)

                                        aaCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgat                                                 <  2:6085330/93‑1 (MQ=255)

                                        aaCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgat                                                 <  1:7498712/93‑1 (MQ=255)

                                        aaCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgat                                                 >  2:2643118/1‑93 (MQ=255)

                                        aaCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTATGGTAGGTAGAGTATTTGAAAAACGTAattgat                                                 >  2:7030016/1‑93 (MQ=255)

                                          cTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgattg                                               <  1:1118630/93‑1 (MQ=255)

                                          cTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgattg                                               <  2:218562/93‑1 (MQ=255)

                                           ttAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGc                                              >  1:5812933/1‑93 (MQ=255)

                                           ttAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGc                                              >  1:2024447/1‑93 (MQ=255)

                                            tAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCg                                             >  1:91930/1‑93 (MQ=255)

                                             aaaGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGa                                            >  1:2419343/1‑93 (MQ=255)

                                             aaaGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGa                                            <  2:4874195/93‑1 (MQ=255)

                                             aaaGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGa                                            >  1:6875656/1‑93 (MQ=255)

                                              aaGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGaa                                           <  2:3805480/93‑1 (MQ=255)

                                              aaGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGaa                                           <  1:2522365/93‑1 (MQ=255)

                                                 cGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAg                                        >  2:3616068/1‑93 (MQ=255)

                                                  ggCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGc                                       <  1:3650188/93‑1 (MQ=255)

                                                   gCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGca                                      >  2:6364459/1‑93 (MQ=255)

                                                   gCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGca                                      >  1:3306094/1‑93 (MQ=255)

                                                    cAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGcaa                                     >  2:776272/1‑93 (MQ=255)

                                                     aGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGcaac                                    >  1:4804441/1‑93 (MQ=255)

                                                     aGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGcaac                                    >  1:6237391/1‑93 (MQ=255)

                                                      gCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGcaaca                                   >  1:7830402/1‑93 (MQ=255)

                                                        ggTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGcaacaaa                                 >  2:6680159/1‑92 (MQ=255)

                                                         gTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACt                                >  1:1913315/1‑93 (MQ=255)

                                                          tCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGAAACAACTa                               >  2:2486470/1‑93 (MQ=255)

                                                            agagAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTaaa                             >  1:4334821/1‑93 (MQ=255)

                                                             gagaATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaac                            <  2:3964871/93‑1 (MQ=255)

                                                             gagaATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaac                            >  2:6211547/1‑93 (MQ=255)

                                                             gagaATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaac                            >  1:4470230/1‑93 (MQ=255)

                                                             gagaATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaac                            <  2:1324723/93‑1 (MQ=255)

                                                               gaATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaacaa                          >  1:5049157/1‑93 (MQ=255)

                                                                aaTGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaacaac                         <  2:7090729/93‑1 (MQ=255)

                                                                 aTGCGCCACCGTAGCAGAAAACAACAGCTTTGTTGTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAAcc                        <  2:3378978/93‑1 (MQ=255)

                                                                 aTGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAAcc                        <  2:4908298/93‑1 (MQ=255)

                                                                  tGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCa                       >  2:5920878/1‑93 (MQ=255)

                                                                    cgcCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAAt                     <  1:6147704/93‑1 (MQ=255)

                                                                     gcCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATa                    <  1:6436820/93‑1 (MQ=255)

                                                                      ccAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATac                   <  1:5389531/93‑1 (MQ=255)

                                                                      ccAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATac                   <  2:1696630/93‑1 (MQ=255)

                                                                       cAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATaca                  >  2:3699633/1‑93 (MQ=255)

                                                                       cAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATaca                  <  1:743095/93‑1 (MQ=255)

                                                                        aaCGTAGCAGAAAACGACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATacac                 >  2:7821253/1‑93 (MQ=255)

                                                                        aaCGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATacac                 <  2:2062147/93‑1 (MQ=255)

                                                                           gTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCg              >  2:6060628/1‑93 (MQ=255)

                                                                           gTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCg              <  2:5847747/93‑1 (MQ=255)

                                                                              gCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGcc           >  2:5740026/1‑93 (MQ=255)

                                                                              gCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGcc           >  1:1385764/1‑93 (MQ=255)

                                                                               cAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGccc          >  2:7213702/1‑93 (MQ=255)

                                                                                aGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGcccc         >  1:1427364/1‑93 (MQ=255)

                                                                                 gAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCt        >  2:570084/1‑93 (MQ=255)

                                                                                 gAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCt        >  1:2963434/1‑93 (MQ=255)

                                                                                 gAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCt        >  2:1218841/1‑93 (MQ=255)

                                                                                 gAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCt        <  2:2398707/93‑1 (MQ=255)

                                                                                 gAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCt        <  1:6992385/93‑1 (MQ=255)

                                                                                  aaaaCAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTc       >  2:2984016/1‑93 (MQ=255)

                                                                                  aaaaCAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTc       >  2:4091305/1‑93 (MQ=255)

                                                                                   aaaCAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCg      >  1:5909910/1‑93 (MQ=255)

                                                                                    aacaacAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGc     >  2:7845286/1‑93 (MQ=255)

                                                                                      caacAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCg   >  1:4085572/1‑93 (MQ=255)

ATCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGG  >  NC_000913/2562261‑2562440

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: `‑`