BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,470,106	C→T	5.2%	intergenic (‑253/+39)	chiA ← / ← tufA	periplasmic endochitinase/translation elongation factor EF‑Tu 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,470,106	0	C	T	5.2%	497.7 / 8.6	149	intergenic (‑253/+39)	chiA/tufA	periplasmic endochitinase/translation elongation factor EF‑Tu 1
Reads supporting (aligned to +/- strand): ref base C (79/62); new base T (6/2); total (85/64)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 149 total aligned reads displayed.

CGCAGGAATTATTTACGTTTTACGAGAACCCTTCAATATGAATTAAATTACGGCATTAAAAATAAGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGCATTACCATCTCTACGCCTTCCGGCAGTTCGATGGTACCAGTCACGTCAGTAGTACGGAAGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGC  >  NC_000913/3469896‑3470377

     gAATTATTTACGTTTTACGAGAACCCTTCAATATGAATTAAATTACGGCATTAAAAATAAGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCcgc                                                                                                                                                                                                             >  2:424781/1‑274 (MQ=255)

                                                      aTTAAAAATAAGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCa                                                                                                                                                           <  2:497769/275‑1 (MQ=255)

                                                                   aaaaaGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAgc                                                                                                                                                                                                                                                                            >  5:1312926/1‑149 (MQ=255)

                                                                         ccTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCtt                                                                                                                                                                                                                                                                      <  7:1295408/149‑1 (MQ=255)

                                                                              cAAATGAAGCATTTTAAAAACAGAAACATTTATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttcaa                                                                                                                                                                                                                                                                 <  6:820633/149‑1 (MQ=255)

                                                                               aaaTGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttcaat                                                                                                                                                                                                                                                                <  5:796735/149‑1 (MQ=255)

                                                                                   gAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttcaattcaa                                                                                                                                                                                                                                                            >  5:237234/1‑149 (MQ=255)

                                                                                    aaGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCaaa                                                                                                                                                                                                                                                           >  7:1247119/1‑149 (MQ=255)

                                                                                      gCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCtt                                                                                                                                                                                                                                                                      >  7:493288/1‑136 (MQ=255)

                                                                                      gCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCTTCGGCGCCCtt                                                                                                                                                                                                                                                                      <  8:493288/136‑1 (MQ=255)

                                                                                          tttAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAAt                                                                                                                                                                                                                                                     <  8:1603651/149‑1 (MQ=255)

                                                                                               aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACg                                                                                                                                                                                                                                                >  5:529385/1‑149 (MQ=255)

                                                                                               aaaCAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACg                                                                                                                                                                                                                                                <  7:613732/149‑1 (MQ=255)

                                                                                                    gAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTATGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGc                                                                                                         <  1:78832/279‑1 (MQ=255)

                                                                                                     aaaCATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAgcg                                                                                                                                                                                                                                                                           >  6:16278/1‑116 (MQ=255)

                                                                                                       aCATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAGAACTAATTAACGTGTAATTa                                                                                                                                                                                                                                        >  5:9052/1‑149 (MQ=255)

                                                                                                         aTTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGc                                                                                                                                                                                                                                      >  5:1545730/1‑149 (MQ=255)

                                                                                                          ttCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGcc                                                                                                                                                                                                                                     >  8:1437528/1‑149 (MQ=255)

                                                                                                           tCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGccc                                                                                                                                                                                                                                    <  5:869611/149‑1 (MQ=255)

                                                                                                                   aaaaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGCATTACCATCCCTACg                                                                                          >  1:174738/1‑279 (MQ=255)

                                                                                                                    aaaTGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTa                                                                                                                                                                                                                           >  8:95175/1‑149 (MQ=255)

                                                                                                                         ttAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCaaaa                                                                                                                                                                                                                                                          <  8:1529777/113‑1 (MQ=255)

                                                                                                                         ttAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTACCGTGTAATTAGCCCAGAACTTTAGcaac                                                                                                                                                                                                                      >  6:467919/1‑149 (MQ=255)

                                                                                                                         ttAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGcaac                                                                                                                                                                                                                      >  6:1307364/1‑149 (MQ=255)

                                                                                                                              ttgaattgaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGACCAGAACTTTAGCAACAACGc                                                                                                                                                                                                                 >  7:696439/1‑149 (MQ=255)

                                                                                                                                 aattgaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACgt                                                                                                                                                                                                                                               >  5:1595973/1‑116 (MQ=255)

                                                                                                                                 aattgaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGCATTACCATCTCAACGCCTTCCGGcc                                                                                >  2:129701/1‑274 (MQ=255)

                                                                                                                                  attgaTATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCcgc                                                                                                                                                                                                             <  7:370914/149‑1 (MQ=255)

                                                                                                                                                tatGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTttcaatt                                                                                                                                                                                                                                                               >  2:122937/1‑85 (MQ=255)

                                                                                                                                                tatGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCaa                                                                                                                                                                                                        <  1:139820/140‑1 (MQ=255)

                                                                                                                                                    aataatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCaaa                                                                                                                                                                                                                                                           <  4:122935/85‑1 (MQ=255)

                                                                                                                                                    aataatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGgccgcc                                                                                                                                                                                           >  6:1096097/1‑149 (MQ=255)

                                                                                                                                                     ataatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCt                                                                                                                                                                                          >  8:798056/1‑149 (MQ=255)

                                                                                                                                                      taatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTa                                                                                                                                                                                                                           <  8:1665580/115‑1 (MQ=255)

                                                                                                                                                      taatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTa                                                                                                                                                                                                                           >  7:1665580/1‑115 (MQ=255)

                                                                                                                                                      taatTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCtt                                                                                                                                                                                         >  8:757197/1‑149 (MQ=255)

                                                                                                                                                        atTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCa                                                                                                                                                                                       >  5:803535/1‑149 (MQ=255)

                                                                                                                                                         tttATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAg                                                                                                                                                                                                                                       >  7:1194154/1‑100 (MQ=255)

                                                                                                                                                            aTTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGa                                                                                                                                                                                   <  7:798056/149‑1 (MQ=255)

                                                                                                                                                            aTTCGTACTGTCAGAACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCaaaa                                                                                                                                                                                                                                                          >  7:1287692/1‑78 (MQ=38)

                                                                                                                                                            aTTCGTACTGTCAGAACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCaaaa                                                                                                                                                                                                                                                          <  8:1287692/78‑1 (MQ=38)

                                                                                                                                                             ttCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAACACTAATTAACgt                                                                                                                                                                                                                                               >  8:1109959/1‑88 (MQ=255)

                                                                                                                                                             ttCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACgt                                                                                                                                                                                                                                               <  7:1109959/88‑1 (MQ=255)

                                                                                                                                                              tCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGt                                                                                                                                     <  3:43484/193‑1 (MQ=255)

                                                                                                                                                              tCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGt                                                                                                                                     >  1:43484/1‑193 (MQ=255)

                                                                                                                                                                gTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGc                                                                                                                                                                                                                                      >  5:629668/1‑94 (MQ=255)

                                                                                                                                                                gTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGc                                                                                                                                                                               <  8:1158775/149‑1 (MQ=255)

                                                                                                                                                                gTTCTGACAGTACGAATAAGATATGCCGACAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGc                                                                                                                                                                                                                                      <  6:629668/94‑1 (MQ=255)

                                                                                                                                                                    tGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGaaa                                                                                                                                                                           >  5:40478/1‑149 (MQ=255)

                                                                                                                                                                    tGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGaaa                                                                                                                                                                           <  7:95175/149‑1 (MQ=255)

                                                                                                                                                                    tGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCACTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGaaa                                                                                                                                                                           <  7:1087498/149‑1 (MQ=255)

                                                                                                                                                                      aCAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACg                                                                                                                                                                         >  5:115295/1‑149 (MQ=255)

                                                                                                                                                                          tACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAGAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGa                                                                                                                                                                     <  6:9052/149‑1 (MQ=255)

                                                                                                                                                                          tACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTaa                                                                                                                                                                                                                                                  >  8:1534890/1‑72 (MQ=255)

                                                                                                                                                                          tACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTaa                                                                                                                                                                                                                                                  <  7:1534890/72‑1 (MQ=255)

                                                                                                                                                                          tACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGa                                                                                                                                                                     >  8:1399892/1‑149 (MQ=255)

                                                                                                                                                                           aCGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGCATTACCATCTCTACGCCTTCCGGCAGTTCGATGGTACCAGTCACGTCAGTAGTACGGAAGTAGAACt                                      >  2:674432/1‑275 (MQ=255)

                                                                                                                                                                           aCGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATACAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGAc                                                                                                                                                                    >  6:1403432/1‑149 (MQ=255)

                                                                                                                                                                              aaTAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGc                                                                                                                                                                                                                 >  6:11158/1‑101 (MQ=255)

                                                                                                                                                                               aTAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgtc                                                                                                                                                                <  6:1472700/149‑1 (MQ=255)

                                                                                                                                                                               aTAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgtc                                                                                                                                                                >  7:1650975/1‑149 (MQ=255)

                                                                                                                                                                                   gATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGc                                                                                                                                                                               <  7:156647/130‑1 (MQ=255)

                                                                                                                                                                                   gATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGc                                                                                                                                                                               >  8:156647/1‑130 (MQ=255)

                                                                                                                                                                                   gATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCa                                                                                                                                                                       <  6:1723492/138‑1 (MQ=255)

                                                                                                                                                                                   gATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCa                                                                                                                                                                       >  5:1723492/1‑138 (MQ=255)

                                                                                                                                                                                   gATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTcc                                                                                                                                                            >  5:654563/1‑149 (MQ=255)

                                                                                                                                                                                    atatGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGCATTAc                                                                                                    <  4:148098/204‑1 (MQ=255)

                                                                                                                                                                                    atatGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGCATTAc                                                                                                    >  2:148100/1‑204 (MQ=255)

                                                                                                                                                                                    atatGCCGTAAACAACATCAAAAAGGCGGCCGAAGGTCCCTTTTCAATTCAAAACTAACTAGCGTGCAATTGGCCCAGACCTTTAGCAACAACTACCGCGCAAAAGGTACGGCCGCCTTCAAGGGTTAAGAAACGCATACCATGATCAATCACGATGGGGTAGATGAGGGTAACAACAAATTTGATGTTGTCGCCCGGCATTACCATCTCTAAGCCTGCCGGCAGTTCGATGGTACCAGTCAAGTCAGTAGTACAGAAGTAGAACTGCGGACGGt                             <  4:129699/275‑1 (MQ=9)

                                                                                                                                                                                       tGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATcg                                                                                                                                                        >  7:1149848/1‑149 (MQ=255)

                                                                                                                                                                                        gCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAAtt                                                                                                                                                                                                                                         >  4:688455/1‑67 (MQ=255)

                                                                                                                                                                                            tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGcc                                                                                                                                                                                                                                     <  2:688460/67‑1 (MQ=255)

                                                                                                                                                                                            tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgtcgt                                                                                                                                                              <  8:824216/138‑1 (MQ=255)

                                                                                                                                                                                            tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACcgtcgt                                                                                                                                                              >  7:824216/1‑138 (MQ=255)

                                                                                                                                                                                            tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATc                                                                                                                                                   >  5:1145897/1‑149 (MQ=255)

                                                                                                                                                                                            tCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGCATTACCATCTCTACGCCTTCCGGCAGTTCGATGGTACCAGTCACGt                                                           <  2:577165/237‑1 (MQ=255)

                                                                                                                                                                                              aacaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCgg                                                                                                                                                 >  7:1538913/1‑149 (MQ=255)

                                                                                                                                                                                              aacaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCgg                                                                                                                                                 >  6:730685/1‑149 (MQ=255)

                                                                                                                                                                                               acaaATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCggg                                                                                                                                                <  8:582240/149‑1 (MQ=255)

                                                                                                                                                                                                 aaaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATc                                                                                                                                                         >  5:654761/1‑138 (MQ=255)

                                                                                                                                                                                                 aaaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATc                                                                                                                                                         <  6:654761/138‑1 (MQ=255)

                                                                                                                                                                                                 aaaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTg                                                                                                                                              >  7:761574/1‑149 (MQ=255)

                                                                                                                                                                                                 aaaTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTg                                                                                                                                              >  8:699123/1‑149 (MQ=255)

                                                                                                                                                                                                   aTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGa                                                                                                                                            <  5:581301/149‑1 (MQ=255)

                                                                                                                                                                                                   aTGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGa                                                                                                                                            <  7:757197/149‑1 (MQ=255)

                                                                                                                                                                                                     gCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATc                                                                                                                                          >  8:1671218/1‑149 (MQ=255)

                                                                                                                                                                                                       gaaaaGGGCGCTTCGGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAg                                                                                                                                        >  5:1512721/2‑149 (MQ=255)

                                                                                                                                                                                                       aaaaaGGGCGCCGAAGCGCCCTTTTCAATTGAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTg                                                                                                                                              >  8:1266784/1‑143 (MQ=255)

                                                                                                                                                                                                       aaaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCa                                                                                                                                                                       >  5:153017/1‑118 (MQ=255)

                                                                                                                                                                                                       aaaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAg                                                                                                                                        >  7:933441/1‑149 (MQ=255)

                                                                                                                                                                                                       aaaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAg                                                                                                                                        >  6:1299148/1‑149 (MQ=255)

                                                                                                                                                                                                        aaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAgg                                                                                                                                       <  5:1455105/149‑1 (MQ=255)

                                                                                                                                                                                                        aaaaGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAgg                                                                                                                                       <  6:1545730/149‑1 (MQ=255)

                                                                                                                                                                                                           aGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGAt                                                                                                                                                    >  2:603854/1‑133 (MQ=255)

                                                                                                                                                                                                           aGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGt                                                                                                                                               <  7:1557262/138‑1 (MQ=255)

                                                                                                                                                                                                             ggCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaac                                                                                                                                  <  8:960771/149‑1 (MQ=255)

                                                                                                                                                                                                              gcgcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTaaca                                                                                                                                 >  6:1574320/1‑149 (MQ=255)

                                                                                                                                                                                                               cgcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCggg                                                                                                                                                <  4:603850/133‑1 (MQ=255)

                                                                                                                                                                                                               cgcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGCATTACCATCTCTACCCCTTCCGGCAATTCGATGGTACCAGTCACGTCCGTAATTCGGAAGTAGAACTGAGGACGGTAGCCTTTGAAGCACGGAGTATGACGgc  >  4:123622/1‑275 (MQ=255)

                                                                                                                                                                                                                gcCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGt                                                                                                                                     <  3:52255/143‑1 (MQ=255)

CGCAGGAATTATTTACGTTTTACGAGAACCCTTCAATATGAATTAAATTACGGCATTAAAAATAAGAAAAAAGCCTGACAAATGAAGCATTTTAAAAACAGAAACATTCATATTTAAAATGTTAAATTGAATTGATATTTTAAATATGAATAATTTATTCGTTCTGACAGTACGAATAAGATATGCCGTCAACAAATGCAAAAAGGGCGCCGAAGCGCCCTTTTCAATTCAAAACTAATTAACGTGTAATTAGCCCAGAACTTTAGCAACAACGCCCGCGCCAACGGTACGGCCGCCTTCACGGATTGCGAAACGCAGACCGTCGTCCATCGCGATCGGGTGGATCAGGGTAACAACCATTTTGATGTTGTCGCCCGGCATTACCATCTCTACGCCTTCCGGCAGTTCGATGGTACCAGTCACGTCAGTAGTACGGAAGTAGAACTGCGGACGGTAGCCTTTGAAGAACGGAGTATGACGGC  >  NC_000913/3469896‑3470377

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