BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,976,402	A→T	7.7%	N42I (AAT→ATT)	wzxE →	lipid III flippase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,976,402	0	A	T	7.7%	476.1 / 22.0	156	N42I (AAT→ATT)	wzxE	lipid III flippase
Reads supporting (aligned to +/- strand): ref base A (71/73); new base T (8/4); total (79/77)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 159 total aligned reads displayed.

TGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCT  >  NC_000913/3976257‑3976550

tgttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑c                                                                                                                                                  >  2:1549321/1‑150 (MQ=255)

 gttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑cc                                                                                                                                                 <  1:181086/150‑1 (MQ=255)

  ttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCg                                                                                                                                                <  1:2426363/150‑1 (MQ=255)

    gAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGcc                                                                                                                                              <  2:978516/150‑1 (MQ=255)

    gAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGcc                                                                                                                                              >  2:2289677/1‑150 (MQ=255)

     aaCTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCAGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCa                                                                                                                                             >  1:1396694/1‑150 (MQ=255)

      actactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAg                                                                                                                                            >  2:838891/1‑150 (MQ=255)

       ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGt                                                                                                                                           >  1:1984181/1‑150 (MQ=255)

       ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTG‑CCGCCAgc                                                                                                                                           <  1:1189632/150‑2 (MQ=255)

        tactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGtt                                                                                                                                          >  1:536477/1‑150 (MQ=255)

         actTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTG‑CCGCCAgccc                                                                                                                                         <  2:668965/150‑4 (MQ=255)

               ccTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTAcc                                                                                                                                   >  1:1310058/1‑150 (MQ=255)

                 tGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGt                                                                                                                                 <  1:803864/150‑1 (MQ=255)

                  gATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTg                                                                                                                                >  1:262585/1‑150 (MQ=255)

                    tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACAATGGTCAAGATTGGTGCCGGGTTACTGGTCGTTAAGTTGCTGGCGGTGTCATATGGTCCGGCGGCGCTTCGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCt                                                                                                                              >  1:593862/1‑150 (MQ=255)

                       gTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCgg                                                                                                                           >  1:490443/1‑150 (MQ=255)

                         cGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCg                                                                                                                         <  1:1691232/150‑1 (MQ=255)

                           ttGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTg                                                                                                                       <  1:939176/150‑1 (MQ=255)

                              gCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                                                                                    >  1:1310920/1‑150 (MQ=255)

                               cAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTg                                                                                                                   >  1:2206665/1‑150 (MQ=255)

                                aaaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGc                                                                                                                  <  1:1200712/150‑1 (MQ=255)

                                 aaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGcc                                                                                                                 <  2:595953/150‑1 (MQ=255)

                                 aaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTG‑CCGCCAgc                                                                                                                                           <  1:1455366/124‑2 (MQ=255)

                                 aaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTG‑CCGCCAgc                                                                                                                                           >  2:1455366/1‑123 (MQ=255)

                                  aaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCg                                                                                                                <  2:496273/150‑1 (MQ=255)

                                   aGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGTTTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCGGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgg                                                                                                               <  1:653727/150‑1 (MQ=255)

                                    gCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                                                                                              <  2:2000575/150‑1 (MQ=255)

                                        ccTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTg                                                                                                          >  1:993223/1‑150 (MQ=255)

                                         cTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTgg                                                                                                         >  2:1749139/1‑150 (MQ=255)

                                         cTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTgg                                                                                                         <  2:29678/150‑1 (MQ=255)

                                          ttGTGGACGGCGGCTAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGc                                                                                                        <  2:1649465/150‑1 (MQ=255)

                                          ttGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTGGCCGGGGCTGGc                                                                                                        <  1:300873/150‑1 (MQ=255)

                                             gggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATc                                                                                                     <  2:1737093/149‑1 (MQ=255)

                                              ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCt                                                                                                    >  1:1671242/1‑150 (MQ=255)

                                              ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCt                                                                                                    >  2:2452292/1‑150 (MQ=255)

                                               gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                                                                                   >  2:1082497/1‑150 (MQ=255)

                                               gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                                                                                   <  1:638821/150‑1 (MQ=255)

                                                    cggcCAGTACACTGGTCAAGATTGGTGCCGGGGTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGGGCTTGCCGGGGCTGGCATCTTTAACg                                                                                              <  2:1585641/150‑1 (MQ=255)

                                                      gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                                                                                            <  2:2206665/150‑1 (MQ=255)

                                                       ccAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                                                                           <  2:98788/150‑1 (MQ=255)

                                                       ccAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                                                                           <  1:1082497/150‑1 (MQ=255)

                                                       ccAGTACACTAGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                                                                           >  2:258206/1‑150 (MQ=255)

                                                          gTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTaa                                                                                        >  2:1774130/1‑150 (MQ=255)

                                                           tACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAc                                                                                       <  2:2412774/150‑1 (MQ=255)

                                                             cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                                                                                     <  1:1783013/150‑1 (MQ=255)

                                                                tGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAAt                                                                                  >  1:70795/1‑150 (MQ=255)

                                                                tGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAAt                                                                                  <  1:759016/150‑1 (MQ=255)

                                                                 ggTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                                                                                 >  2:366075/1‑150 (MQ=255)

                                                                   tCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACg                                                                               >  1:232975/1‑150 (MQ=255)

                                                                    cAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt                                                                              >  2:2100293/1‑150 (MQ=255)

                                                                       gATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc                                                                           >  1:2429856/1‑150 (MQ=255)

                                                                       gATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc                                                                           >  2:968842/1‑150 (MQ=255)

                                                                            gTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGt                                                                      <  2:2038964/150‑1 (MQ=255)

                                                                            gTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGt                                                                      >  1:432826/1‑150 (MQ=255)

                                                                             tGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa                                                                     <  2:2180691/150‑1 (MQ=255)

                                                                              gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCTGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc                                                                    <  2:1776832/150‑1 (MQ=255)

                                                                              gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc                                                                    >  2:104370/1‑150 (MQ=255)

                                                                              gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc                                                                    <  1:31201/150‑1 (MQ=255)

                                                                                  ggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATg                                                                >  1:1574926/1‑150 (MQ=255)

                                                                                   ggtACTGGTCGGTAAGTTGCGGGCGGTGGCATTTGGTCCGGCGGGGCTTGGGCGGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGa                                                               <  2:1596953/148‑1 (MQ=255)

                                                                                    ttACTGGTCGGTAAGTTGCTGGCGGTGTCTTTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATAT‑CCGCCAGTTGATGTCCGGGCTCGGCGTGCTTGCCTGGGCGGGCAGCTTTAACGGTGTAACCAAATACGTTGCCCAGGACCATGAt                                                              <  2:1512855/150‑1 (MQ=255)

                                                                                    ttACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGAt                                                              <  1:258206/150‑1 (MQ=255)

                                                                                       cTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAAt                                                           >  1:38108/1‑150 (MQ=255)

                                                                                         ggTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc                                                         <  2:1621162/150‑1 (MQ=255)

                                                                                               tAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACa                                                   >  1:962161/1‑150 (MQ=255)

                                                                                                  gTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCt                                                >  2:2356232/1‑150 (MQ=255)

                                                                                                  gTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCt                                                <  1:2088087/150‑1 (MQ=255)

                                                                                                      cTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTgcgc                                            >  2:1314901/1‑150 (MQ=255)

                                                                                                          cGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCcgcg                                        <  1:70076/150‑1 (MQ=255)

                                                                                                          cGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCcgcg                                        <  1:61159/150‑1 (MQ=255)

                                                                                                           ggTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGt                                       <  2:2153066/150‑1 (MQ=255)

                                                                                                             tgtCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTgg                                     <  1:466190/150‑1 (MQ=255)

                                                                                                               tCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTc                                   <  1:2289677/150‑1 (MQ=255)

                                                                                                                cATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCg                                  <  2:1113920/150‑1 (MQ=255)

                                                                                                                 aTTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCgg                                 >  2:2071659/1‑150 (MQ=255)

                                                                                                                  tttGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGc                                <  1:1549321/150‑1 (MQ=255)

                                                                                                                    tGGTCCTGCGGGGCTTGGGCTGGCGGCAAATTT‑CAGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCAc                              >  2:1380874/1‑150 (MQ=255)

                                                                                                                    tGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCAc                              >  2:574940/1‑150 (MQ=255)

                                                                                                                    tGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCAc                              <  1:573743/150‑1 (MQ=255)

                                                                                                                    tGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCAc                              >  2:1282013/1‑150 (MQ=255)

                                                                                                                      gTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACtt                            >  1:1046452/1‑150 (MQ=255)

                                                                                                                        ccGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTtca                          >  1:207476/1‑150 (MQ=255)

                                                                                                                             ggggCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGc                     <  2:253194/150‑1 (MQ=255)

                                                                                                                                gCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGag                  >  2:72295/1‑149 (MQ=255)

                                                                                                                                gCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGAt                  >  2:1191681/1‑150 (MQ=255)

                                                                                                                                  ttGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATgg                <  1:1701131/150‑1 (MQ=255)

                                                                                                                                   tcaacTGGCGGAAATTTG‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGt               >  2:882242/5‑150 (MQ=255)

                                                                                                                                   tGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGt               <  1:1585393/150‑1 (MQ=255)

                                                                                                                                     aacTGGCGGAAATTTG‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTAc             >  1:767884/3‑150 (MQ=255)

                                                                                                                                      gCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt            >  2:1062296/1‑150 (MQ=255)

                                                                                                                                      acTGGCGGAAATTTG‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt            >  2:2073634/2‑150 (MQ=255)

                                                                                                                                      acTGGCGGAAATTTG‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt            >  1:1457183/2‑150 (MQ=255)

                                                                                                                                      acTGGCGGAAATTTG‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACAATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt            >  2:1932441/2‑150 (MQ=255)

                                                                                                                                       cTGGCGGAAATTTG‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAACAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACtt           >  1:2449979/1‑150 (MQ=255)

                                                                                                                                          gcggcAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGt        <  1:1980541/150‑1 (MQ=255)

                                                                                                                                           cggcAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGtt       <  2:1396694/150‑1 (MQ=255)

                                                                                                                                           cggcAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGtt       <  2:767036/150‑1 (MQ=255)

                                                                                                                                             gcAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCAATGGTACTTGGTTtc     <  1:1638351/150‑1 (MQ=255)

                                                                                                                                              gaaaTTTGCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTtc     >  1:176273/2‑150 (MQ=255)

                                                                                                                                                aaTTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTctct  <  2:145890/150‑1 (MQ=255)

TGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTT‑CCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCT  >  NC_000913/3976257‑3976550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`