BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,693,594	(C)_5→4	intergenic (‑4/+32)	yhjQ ← / ← yhjR	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,693,590	0	C	.	82.5%	926.1 / 201.2	315	coding (1/729 nt)	yhjQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (28/27); new base . (138/122); total (166/149)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

Only 100 of 333 total aligned reads displayed.

GCAGGCATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAA  >  NC_000913/3693498‑3693682

   ggcaTCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCAcccc                                                                                           <  2:1731280/93‑1 (MQ=255)

    gcaTCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCAcccct                                                                                          <  1:3626323/93‑2 (MQ=255)

    gcaTCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCAccccc                                                                                          <  1:6868691/93‑1 (MQ=255)

      aTCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCAcccctgc                                                                                        >  1:5595902/1‑90 (MQ=255)

        cGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTa                                                                                     <  2:5231401/93‑1 (MQ=255)

            caccaGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATccc                                                                                 >  1:5512557/1‑93 (MQ=255)

            caccaGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATccc                                                                                 >  2:4983985/1‑93 (MQ=255)

               caGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGt                                                                              <  2:388214/93‑1 (MQ=255)

                aGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTa                                                                             <  1:4421984/93‑1 (MQ=255)

                 ggACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTAc                                                                            >  2:4088566/1‑93 (MQ=255)

                 ggACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTAc                                                                            >  1:2140311/1‑93 (MQ=255)

                 ggACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTAc                                                                            >  1:383691/1‑93 (MQ=255)

                     aTTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGcc                                                                        >  1:3316969/1‑93 (MQ=255)

                      ttttCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCa                                                                       >  1:4938685/1‑93 (MQ=255)

                        ttCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCAt                                                                      >  2:6747918/1‑93 (MQ=255)

                          ctcCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTc                                                                    <  1:4537715/93‑1 (MQ=255)

                          ctcCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCa                                                                   <  2:7236036/93‑1 (MQ=255)

                           tcCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCAc                                                                  >  2:2607996/1‑93 (MQ=255)

                            cccAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCAcc                                                                 >  2:5840239/1‑93 (MQ=255)

                            cccAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCAcc                                                                 >  1:5041856/1‑93 (MQ=255)

                             ccAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCa                                                                >  2:2761461/1‑93 (MQ=255)

                              cAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCaa                                                               <  1:2177025/93‑1 (MQ=255)

                               aaCATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAAt                                                              <  1:3257218/93‑1 (MQ=255)

                                aCATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATc                                                             <  2:6322606/93‑1 (MQ=255)

                                  aTTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATccc                                                           >  2:7527806/1‑93 (MQ=255)

                                     tGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTa                                                        >  1:912154/1‑93 (MQ=255)

                                      gTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTAc                                                       <  2:6145947/93‑1 (MQ=255)

                                        aaTGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACtt                                                     <  1:7420771/93‑1 (MQ=255)

                                         aTGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACttt                                                    <  1:4845163/93‑1 (MQ=255)

                                         aTGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACAGCCATTCACCAATCCCCTACttt                                                    >  1:1544954/1‑93 (MQ=255)

                                          tGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACtttt                                                   <  2:6203593/93‑1 (MQ=255)

                                           gACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTttg                                                  >  2:1413642/1‑93 (MQ=255)

                                           gACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCACAGTACGGCCATTCACCAATCCCCTACTTttg                                                  >  2:3578790/1‑93 (MQ=255)

                                            aCCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTttgt                                                 >  2:5322233/1‑93 (MQ=255)

                                               aGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTgc                                              <  1:978155/93‑1 (MQ=255)

                                               aGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTgc                                              >  2:4647107/1‑93 (MQ=255)

                                               aGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTgc                                              <  2:1930961/93‑1 (MQ=255)

                                                ggCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTgc                                              <  2:3211131/93‑1 (MQ=255)

                                                ggCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTgc                                              <  1:6895048/93‑1 (MQ=255)

                                                 gCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTgcgc                                            <  2:3137849/93‑1 (MQ=255)

                                                 gCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTgcgc                                            <  2:6925096/93‑1 (MQ=255)

                                                  cTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCa                                           >  2:5083485/1‑93 (MQ=255)

                                                    aaTGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCaa                                          <  2:3925279/93‑1 (MQ=255)

                                                    aaTGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCaaa                                         >  2:4314687/1‑93 (MQ=255)

                                                    aaTGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCaaa                                         <  1:1899075/93‑1 (MQ=255)

                                                    aaTGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCaaa                                         <  2:2977344/93‑1 (MQ=255)

                                                     aTGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCaaa                                         >  2:2104874/1‑93 (MQ=255)

                                                     aTGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAAc                                        <  2:1567920/93‑1 (MQ=255)

                                                      tGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCaann                                        >  1:1401357/1‑91 (MQ=255)

                                                      tGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAAct                                       >  1:1201666/1‑93 (MQ=255)

                                                       gcggcgGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAAct                                       <  2:2544681/93‑1 (MQ=255)

                                                        cggcggTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAActct                                     <  1:2283388/93‑1 (MQ=255)

                                                        cggcggTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAActct                                     <  1:6523654/93‑1 (MQ=255)

                                                          gcggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTg                                    <  2:743827/93‑1 (MQ=255)

                                                          gcggTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGc                                   >  1:4762892/1‑93 (MQ=255)

                                                          gcggTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGc                                   <  2:2260797/93‑1 (MQ=255)

                                                           cggTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGcc                                  <  1:6296704/93‑1 (MQ=255)

                                                             gTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAg                                >  2:7336986/1‑93 (MQ=255)

                                                             gTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAg                                >  2:4779622/1‑93 (MQ=255)

                                                             gTGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAg                                >  1:1146790/1‑93 (MQ=255)

                                                              tGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAg                                >  2:7059094/1‑93 (MQ=255)

                                                              tGATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGc                               <  2:4514700/93‑1 (MQ=255)

                                                               gATGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCa                              >  1:7103634/1‑93 (MQ=255)

                                                                aTGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCaa                             <  1:1460408/93‑1 (MQ=255)

                                                                aTGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCaa                             >  1:4083283/1‑93 (MQ=255)

                                                                aTGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCaa                             <  2:6051134/93‑1 (MQ=255)

                                                                 tGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCaa                             <  2:3784338/93‑1 (MQ=255)

                                                                 tGGTTGTTGTCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAAc                            <  1:2945301/93‑1 (MQ=255)

                                                                   gttgttGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACg                           >  2:35220/1‑93 (MQ=255)

                                                                    ttgttgTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACgg                          >  1:3475661/1‑93 (MQ=255)

                                                                      gttgtCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCa                       <  1:586105/93‑1 (MQ=255)

                                                                      gttgtCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCa                       <  2:6576372/93‑1 (MQ=255)

                                                                       ttgtCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAg                      <  2:5736497/93‑1 (MQ=255)

                                                                        tgtCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGc                     <  1:3150717/93‑1 (MQ=255)

                                                                          tCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGt                   <  1:2609673/93‑1 (MQ=255)

                                                                          tCCCCACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGt                   >  1:7224215/1‑93 (MQ=255)

                                                                           ccccACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGt                   <  1:3362213/93‑1 (MQ=255)

                                                                            cccACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGttt                 <  2:3340503/93‑1 (MQ=255)

                                                                            cccACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGttt                 >  2:2835804/1‑93 (MQ=255)

                                                                             ccACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGtttt                >  2:4186120/1‑93 (MQ=255)

                                                                             ccACGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGtttt                >  1:4212890/1‑93 (MQ=255)

                                                                               aCGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTa               <  2:5929026/93‑1 (MQ=255)

                                                                               aCGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTaa              <  1:4824493/93‑1 (MQ=255)

                                                                               aCGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTaa              >  2:2483206/1‑93 (MQ=255)

                                                                                cGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTaa              <  1:963648/93‑1 (MQ=255)

                                                                                cGCCTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAAt             <  2:5568582/93‑1 (MQ=255)

                                                                                  ccTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGc           <  1:2816983/93‑1 (MQ=255)

                                                                                  ccTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGc           <  1:615191/93‑1 (MQ=255)

                                                                                  ccTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGc           >  1:6695192/1‑93 (MQ=255)

                                                                                   cTCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGcc          >  2:81423/1‑93 (MQ=255)

                                                                                    tCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcg         >  2:4798608/1‑93 (MQ=255)

                                                                                    tCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcg         <  1:5238620/93‑1 (MQ=255)

                                                                                    tCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcg         >  1:53259/1‑93 (MQ=255)

                                                                                    tCCCCGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcg         <  1:6573993/93‑1 (MQ=255)

                                                                                     ccccGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcgc        <  1:469044/93‑1 (MQ=255)

                                                                                     ccccGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcgc        <  1:1830440/93‑1 (MQ=255)

                                                                                     ccccGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcgc        <  2:1061692/93‑1 (MQ=255)

                                                                                      cccGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcgcg       >  2:6005811/1‑93 (MQ=255)

                                                                                      cccGCA‑CCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCcgcg       >  2:933837/1‑93 (MQ=255)

                                                                                          cACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGcc    >  1:7440877/1‑93 (MQ=255)

GCAGGCATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAA  >  NC_000913/3693498‑3693682

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