BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,382,000	C→T	6.0%	C18C (TGC→TGT)	ompG →	outer membrane porin G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,382,000	0	C	T	6.0%	542.8 / 16.7	167	C18C (TGC→TGT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand): ref base C (81/76); new base T (7/3); total (88/79)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 171 total aligned reads displayed.

TGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATT  >  NC_000913/1381853‑1382144

  aCGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATTAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGca                                                                                                                                               >  1:278953/1‑149 (MQ=255)

  aCGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGca                                                                                                                                               >  2:103347/1‑149 (MQ=255)

     aaaTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg                                                                                                                                            <  1:1059951/149‑1 (MQ=255)

      aaTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAgg                                                                                                                                           >  4:540217/1‑149 (MQ=255)

         gTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGcca                                                                                                                                        <  3:831528/149‑2 (MQ=255)

          tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGccat                                                                                                                                       <  1:1093059/149‑3 (MQ=255)

          tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGccat                                                                                                                                       <  2:325661/149‑3 (MQ=255)

             tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagga                                                                                                                                    >  3:1331858/1‑149 (MQ=255)

             tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGCGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagga                                                                                                                                    >  2:1528533/1‑149 (MQ=255)

               tCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                  <  3:799011/149‑1 (MQ=255)

               tCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                  <  1:484904/149‑1 (MQ=255)

                 tttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagg                                                                                                                                >  4:385640/1‑149 (MQ=255)

                 tttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagg                                                                                                                                >  1:82508/1‑149 (MQ=255)

                  ttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGACATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagga                                                                                                                               <  2:820718/149‑1 (MQ=255)

                      tGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGa                                                                                                                           >  2:37325/1‑149 (MQ=255)

                             aCGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTgc                                                                                                                                                  >  4:697954/1‑119 (MQ=255)

                             aCGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTgc                                                                                                                                                  <  3:697954/119‑1 (MQ=255)

                             aCGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCAc                                                                                                                    >  1:195629/1‑149 (MQ=255)

                                     aGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                  >  2:692767/1‑127 (MQ=255)

                                       cAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCg                                                                                                          <  2:1091888/149‑1 (MQ=255)

                                       cAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATCACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACCGGCACTTTAATATCg                                                                                                          <  4:1160251/149‑1 (MQ=255)

                                           tCGATAAATACAGGTGGAAGGCCTGCACCCAGGATAATATAAGGTACAATAATGAAAAAGTTATAACTCTGTACCGCACTGGTGATGTGTGCCGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCAGTTTAATATCGgcgc                                                                                                      <  2:217558/149‑1 (MQ=255)

                                             gCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGa                                                                                                    >  2:1154310/1‑149 (MQ=255)

                                                   tACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTAc                                                                                               <  1:130889/148‑1 (MQ=255)

                                                    aCAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                       >  1:859409/1‑107 (MQ=255)

                                                      tgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                           >  4:1435133/2‑149 (MQ=255)

                                                      aGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                           <  4:1061196/149‑1 (MQ=255)

                                                      aGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                           >  2:491558/1‑149 (MQ=255)

                                                       gggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAAt                                                                                          >  4:924480/1‑149 (MQ=255)

                                                        ggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATa                                                                                         <  1:552113/149‑1 (MQ=255)

                                                        ggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATa                                                                                         <  4:370712/149‑1 (MQ=255)

                                                             aGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGaaaa                                                                                    <  2:1418573/149‑1 (MQ=255)

                                                             aGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGaaaa                                                                                    >  1:1418573/1‑149 (MQ=255)

                                                                  ttCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                       >  3:44468/1‑93 (MQ=255)

                                                                  ttCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                       <  4:44468/93‑1 (MQ=255)

                                                                     ccccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAgg                                                                            <  4:343671/149‑1 (MQ=255)

                                                                     ccccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAgg                                                                            >  3:1411352/1‑149 (MQ=255)

                                                                      cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg                                                                                                                                            >  3:956287/1‑84 (MQ=255)

                                                                      cccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAggg                                                                           >  3:854339/1‑149 (MQ=255)

                                                                       ccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGt                                                                          >  1:729631/1‑149 (MQ=255)

                                                                         aGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTAGAGGGTTa                                                                        <  3:1103117/149‑1 (MQ=255)

                                                                          ggATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTAt                                                                       <  4:80595/149‑1 (MQ=255)

                                                                            ataataCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATgg                                                                     <  4:854339/149‑1 (MQ=255)

                                                                              aataCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCg                                                                   >  3:528785/1‑149 (MQ=255)

                                                                               ataCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                  >  2:1168577/1‑85 (MQ=255)

                                                                                taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGaa                                                                 <  1:895050/149‑1 (MQ=255)

                                                                                taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGaa                                                                 <  3:1435133/149‑1 (MQ=255)

                                                                                   aaGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcgc                                                                                                      >  3:565766/1‑109 (MQ=255)

                                                                                     ggAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGatat                                                            >  2:1369835/1‑149 (MQ=255)

                                                                                       aaCAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatgg                                                          <  4:1104664/149‑1 (MQ=255)

                                                                                         cAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggat                                                        <  2:880191/149‑1 (MQ=255)

                                                                                          aataatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatg                                                       >  1:475301/1‑149 (MQ=255)

                                                                                           ataatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatgg                                                      <  3:385640/149‑1 (MQ=255)

                                                                                            taatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATggg                                                     <  2:475301/149‑1 (MQ=255)

                                                                                            aaatGGAAAACTTAGTACCCTGTACCGCACTGGAGACGTGTGCGGGATTGGGCTGCGCACAGGAAGAGGGAAGGAAAGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAACACGTCGAGGGTTATGGCGAAGATATGGATggg                                                     >  2:405250/2‑149 (MQ=255)

                                                                                             aatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGc                                                    <  3:1495876/149‑1 (MQ=255)

                                                                                               tGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTg                                                  >  3:1201827/1‑149 (MQ=255)

                                                                                                   aaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCgg                                              >  2:1274013/1‑149 (MQ=255)

                                                                                                   aaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCgg                                              >  4:77017/1‑149 (MQ=255)

                                                                                                   aaaGTTATTACCCTGTACCGCACTGGTGATGAGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCgg                                              >  2:728297/1‑149 (MQ=255)

                                                                                                    aaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGa                                             >  4:183956/1‑149 (MQ=255)

                                                                                                     aGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAg                                            <  2:927933/149‑1 (MQ=255)

                                                                                                      gTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGc                                           <  2:195629/149‑1 (MQ=255)

                                                                                                           taCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAGGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATACAGAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCtt                                        >  4:145211/1‑147 (MQ=255)

                                                                                                           taCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCtt                                        <  3:145211/147‑1 (MQ=255)

                                                                                                           taCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCa                                      <  2:1525139/149‑1 (MQ=255)

                                                                                                           taCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCa                                      <  1:386066/149‑1 (MQ=255)

                                                                                                            aCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGCAGATATGGATGGGCTGGCGGAGCCCTCAg                                     >  2:1543082/1‑149 (MQ=255)

                                                                                                             cccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGt                                    >  4:1241102/1‑149 (MQ=255)

                                                                                                             cccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGt                                    >  2:849628/1‑149 (MQ=255)

                                                                                                              ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACATCGAGGGTTATGGCGAAGATatgg                                                          <  2:955324/126‑1 (MQ=255)

                                                                                                              ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACATCGAGGGTTATGGCGAAGATatgg                                                          >  1:955324/1‑126 (MQ=255)

                                                                                                               cTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAAt                                                                                          <  4:841992/93‑1 (MQ=255)

                                                                                                               cTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCt                                  >  4:752649/1‑149 (MQ=255)

                                                                                                                tGTACCGCCCTGGTGATGTGTGCGGGAATGGCCTGCGCACACGCCGAGGGAAGTAACGACTGGCACTTTAATATCGGCGCGAGGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATAGGCTGCCGGAGCCTTCAGGCTa                                 >  4:1083707/1‑149 (MQ=255)

                                                                                                                tGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatg                                                           >  4:1199469/1‑123 (MQ=255)

                                                                                                                tGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTa                                 <  1:706050/149‑1 (MQ=255)

                                                                                                                 gTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTAt                                <  3:924480/149‑1 (MQ=255)

                                                                                                                     cGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTa                            >  4:860764/1‑149 (MQ=255)

                                                                                                                      gCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTaa                           <  3:457271/149‑1 (MQ=255)

                                                                                                                      gCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTaa                           <  1:1519069/149‑1 (MQ=255)

                                                                                                                          tGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATgcc                       >  2:786452/1‑149 (MQ=255)

                                                                                                                               aTGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCaa                  >  3:923814/1‑149 (MQ=255)

                                                                                                                                tgtgtgCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAAc                 >  3:1291383/1‑149 (MQ=255)

                                                                                                                                  tgtgCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACgg               <  3:805485/149‑1 (MQ=255)

                                                                                                                                        ggAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTg         <  4:1499214/149‑1 (MQ=255)

                                                                                                                                         cctcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTgg        >  3:161451/5‑149 (MQ=255)

                                                                                                                                          ctcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGga       >  4:35379/4‑149 (MQ=255)

                                                                                                                                          ctcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGga       >  1:1383082/4‑149 (MQ=255)

                                                                                                                                          ctcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGga       >  3:1115509/4‑149 (MQ=255)

                                                                                                                                          aaTGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGga       <  1:1528533/149‑1 (MQ=255)

                                                                                                                                          aaTGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGga       <  1:786452/149‑1 (MQ=255)

                                                                                                                                           tcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGgag      >  1:1237446/3‑149 (MQ=255)

                                                                                                                                           tcggCCTGTGCGCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGgag      >  1:1277667/3‑149 (MQ=255)

                                                                                                                                           aTGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATgc                        >  1:234012/1‑131 (MQ=255)

                                                                                                                                           aTGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATgc                        <  2:234012/131‑1 (MQ=255)

                                                                                                                                             ggCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGaa    <  2:457086/149‑1 (MQ=255)

                                                                                                                                             ggCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGaa    <  4:1291383/149‑1 (MQ=255)

                                                                                                                                               ccTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAATATATGGATGGGCTGGCGGAGCCTTCAGTCCATTTTAATGCCGCCAACGGGCCGTGGAGAAtt  >  3:1226362/1‑149 (MQ=255)

                                                                                                                                               ccTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAAtt  >  1:905829/1‑149 (MQ=255)

TGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATT  >  NC_000913/1381853‑1382144

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