BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,976,402	A→T	6.6%	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	6.6%	379.3 / 10.5	120	N42I (AAT→ATT)	wzxE	lipid III flippase
Reads supporting (aligned to +/- strand): ref base A (59/53); new base T (4/4); total (63/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 121 total aligned reads displayed.

CTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTC  >  NC_000913/3976254‑3976549

cTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAt                                                                                                                                                     <  1:1238966/150‑1 (MQ=255)

cTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAt                                                                                                                                                     <  2:579203/150‑1 (MQ=255)

    gttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTcc                                                                                                                                                 <  2:1282822/150‑1 (MQ=255)

      tgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGc                                                                                                                                               >  2:702006/1‑150 (MQ=255)

        aaCTACTTTTCCTGATATGGCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCa                                                                                                                                             <  2:1426358/150‑1 (MQ=255)

         actactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑g                                                                                                                                           >  1:347713/1‑150 (MQ=255)

         actactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑g                                                                                                                                           >  2:235435/1‑150 (MQ=255)

          ctactTTTCCTGATATGTCGTTGGCAAAAGCGGCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGGGCCGGGTGACTGGTCGGTAAGTTGCTGGCGGGGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCA‑gc                                                                                                                                          <  1:1000793/150‑2 (MQ=255)

           tactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCA‑gcc                                                                                                                                         <  1:539126/150‑3 (MQ=255)

           tactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCA‑gcc                                                                                                                                         <  1:686329/150‑3 (MQ=255)

              ttttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGAt                                                                                                                                      <  1:839926/150‑1 (MQ=255)

              ttttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGAt                                                                                                                                      >  1:637692/1‑150 (MQ=255)

                ttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTa                                                                                                                                    <  1:428169/150‑1 (MQ=255)

                   cTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCg                                                                                                                                 >  1:800447/1‑150 (MQ=255)

                        atGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTc                                                                                                                            >  2:362884/1‑150 (MQ=255)

                         tGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCg                                                                                                                           <  2:515029/150‑1 (MQ=255)

                          gTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCgg                                                                                                                          >  2:191915/1‑150 (MQ=255)

                          gTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCgg                                                                                                                          >  2:922933/1‑150 (MQ=255)

                               tGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGc                                                                                                                     >  1:322297/1‑150 (MQ=255)

                                 gCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCtt                                                                                                                   >  2:772612/1‑150 (MQ=255)

                                  cAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTg                                                                                                                  >  1:402146/1‑150 (MQ=255)

                                  cAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTg                                                                                                                  >  1:556978/1‑150 (MQ=255)

                                          tCCTTGTGGACGGCGGCCAGTACACTGGTCTAGTTTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCTTTTGGTCCGGCGGGGCTTTGTCTTGCGGCAAATTTCCGCGA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                                                                                          <  2:179790/150‑1 (MQ=255)

                                          tCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                                                                                          >  2:458241/1‑150 (MQ=255)

                                            cTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTgg                                                                                                        <  1:867760/150‑1 (MQ=255)

                                             ttGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGATTGCCGGGGCTGGc                                                                                                       >  1:494555/1‑150 (MQ=255)

                                              tgtgGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCa                                                                                                      <  1:643827/150‑1 (MQ=255)

                                               gtgGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCAt                                                                                                     >  1:701435/1‑150 (MQ=255)

                                                 ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCt                                                                                                   >  1:585476/1‑150 (MQ=255)

                                                  gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                                                                                  >  1:1034895/1‑150 (MQ=255)

                                                   aCGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCttt                                                                                                 >  2:274465/1‑150 (MQ=255)

                                                    cggcggCCAGTACACTGGTCAAGATTGGTTCCTTGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTTGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGTCGTGCTTGCCGGGGCTGGCATCTTTa                                                                                                <  2:780480/150‑1 (MQ=255)

                                                     ggcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTaa                                                                                               >  1:1393129/1‑150 (MQ=255)

                                                      gcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAc                                                                                              >  1:762198/1‑150 (MQ=255)

                                                       cggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACg                                                                                             <  1:274465/150‑1 (MQ=255)

                                                        ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTTGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                                                                                            <  2:172105/150‑1 (MQ=255)

                                                        ggcCAGTACACTGGTCAAGATTGGGGCCGGGGTACTGGTCGGTAAGTTGCTGGCGGTGGCATTTGGTCCGGCGGGGCTTGGGCTGGCGCCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                                                                                            <  1:941926/150‑1 (MQ=255)

                                                          ccAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                                                                          <  2:1390965/150‑1 (MQ=255)

                                                          ccAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                                                                          <  2:910280/150‑1 (MQ=255)

                                                             gTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTaa                                                                                               >  2:742557/1‑142 (MQ=255)

                                                             gTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTaa                                                                                               <  1:742557/142‑1 (MQ=255)

                                                              tACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTGCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAc                                                                                      <  2:1329512/150‑1 (MQ=255)

                                                              tACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAc                                                                                      <  2:1034895/150‑1 (MQ=255)

                                                                cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTg                                                                                                                      >  2:43857/1‑116 (MQ=255)

                                                                cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTg                                                                                                                      <  1:43857/116‑1 (MQ=255)

                                                                cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                                                                                    >  2:1075409/1‑150 (MQ=255)

                                                                cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                                                                                    >  2:86923/1‑150 (MQ=255)

                                                                cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                                                                                    >  2:1019055/1‑150 (MQ=255)

                                                                   tGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAAt                                                                                 <  1:772612/150‑1 (MQ=255)

                                                                    ggTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                                                                                <  1:1242227/150‑1 (MQ=255)

                                                                    ggTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                                                                                >  2:494777/1‑150 (MQ=255)

                                                                    ggTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                                                                                <  1:255214/150‑1 (MQ=255)

                                                                         aGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTg                                                                           >  2:792420/1‑150 (MQ=255)

                                                                            ttGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTTCTCGGCGTGCTTGCCGGTGCTGGAATCTTTAACGGTGTAACCAAATAGGTTgact                                                                        <  2:701435/150‑4 (MQ=255)

                                                                              ggTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCCTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAg                                                                      >  1:442387/1‑150 (MQ=255)

                                                                                tGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa                                                                    <  1:1058434/150‑1 (MQ=255)

                                                                                tGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa                                                                    >  1:1335558/1‑150 (MQ=255)

                                                                                tGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa                                                                    >  2:1226224/1‑150 (MQ=255)

                                                                                 gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCggg                                                                                                             <  1:983271/108‑1 (MQ=255)

                                                                                 gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCggg                                                                                                             >  2:983271/1‑108 (MQ=255)

                                                                                 gCCGGGTAACTGGTCGGTAAGTTGCTGGCGGTCTCATTTGGTCCGGCGGGGATTGGGCTGGAGGCAAAATTCCGCCA‑GTTCATAACCGTGATCGGAGTGCTAGCCGCGCCTGGCATCTTTAAAGGAGTAACCAAATACGTTGCCCAGTAc                                                                   >  2:144975/1‑150 (MQ=255)

                                                                                  ccGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAcc                                                                  >  1:757541/1‑150 (MQ=255)

                                                                                   ggggTTACTGGTCGGTAAGTTGCTGGCTGTGTCATTTGGTCCGGCGGGGCTTGGGCTTGCGGCAAATTTCCGCCA‑GTTGATTACCGTTCTCGGCGTTCTTGCCGTGGCTGGCATCTTTATCTGTGTAACGAAATACGTTGCCCAGTAGCa                                                                 <  1:144975/149‑1 (MQ=255)

                                                                                   cGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCa                                                                 >  2:287006/1‑150 (MQ=255)

                                                                                     ggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATg                                                               <  1:260399/150‑1 (MQ=255)

                                                                                      gTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGa                                                              <  2:299409/150‑1 (MQ=255)

                                                                                         aCTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATaa                                                           >  1:1088369/1‑150 (MQ=255)

                                                                                         aCTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATaa                                                           <  1:235435/150‑1 (MQ=255)

                                                                                           tGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATc                                                         >  2:461422/1‑150 (MQ=255)

                                                                                            ggTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc                                                        >  1:285506/1‑150 (MQ=255)

                                                                                              tCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGc                                                      <  2:1393129/150‑1 (MQ=255)

                                                                                                   aaGTTGGTGGCGGGGTCTTTTGGTCGGGGGGGGCTTGGGCTGGCGGC‑AATTTTCGCCATTTTTATTACGGTGCTGGGGGTGCTTGCGGGGGGTGG‑ATTATTAACGGTGTAACCAAATACTTTGCCCTGTACCATGATAATCCGCAACAGc                                                <  2:1223831/150‑1 (MQ=255)

                                                                                                   aaGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAg                                                 >  2:1360479/1‑150 (MQ=255)

                                                                                                        gCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTgcg                                            >  2:749686/1‑150 (MQ=255)

                                                                                                        gCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTgcg                                            <  1:764582/150‑1 (MQ=255)

                                                                                                         cTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc                                                        <  1:3094/137‑1 (MQ=255)

                                                                                                         cTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc                                                        >  2:3094/1‑137 (MQ=255)

                                                                                                           ggCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGcgccg                                         <  2:1301167/150‑1 (MQ=255)

                                                                                                             cGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGc                                                      >  1:918239/1‑135 (MQ=255)

                                                                                                             cGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGc                                                      <  2:918239/135‑1 (MQ=255)

                                                                                                                    aTTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCgg                                >  1:1194566/1‑150 (MQ=255)

                                                                                                                     tttGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGc                               <  2:285506/150‑1 (MQ=255)

                                                                                                                        ggTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACt                            >  2:1276762/1‑150 (MQ=255)

                                                                                                                         gTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACtt                           >  2:13273/1‑150 (MQ=255)

                                                                                                                         gTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACtt                           >  1:263291/1‑150 (MQ=255)

                                                                                                                           ccGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTtca                         >  1:809920/1‑150 (MQ=255)

                                                                                                                                ggggCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGc                    <  2:762198/150‑1 (MQ=255)

                                                                                                                                     ttGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATgg               <  1:1290275/150‑1 (MQ=255)

                                                                                                                                      tGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGt              <  1:191915/150‑1 (MQ=255)

                                                                                                                                        aacTGGCGGAAATTTGCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTAc            >  2:225679/3‑150 (MQ=255)

                                                                                                                                         gCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt           >  2:710149/1‑150 (MQ=255)

                                                                                                                                         gCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt           >  2:54174/1‑150 (MQ=255)

                                                                                                                                         acTGGCGGAAATTTGCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt           >  2:1179048/2‑150 (MQ=255)

                                                                                                                                         acTGGCGGAAATTTGCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt           >  1:16054/2‑150 (MQ=255)

                                                                                                                                           tGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTg         <  2:1194566/150‑1 (MQ=255)

                                                                                                                                            ggcggcAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTgg        <  1:702006/150‑1 (MQ=255)

                                                                                                                                              cggcAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGtt      <  1:458241/150‑1 (MQ=255)

                                                                                                                                                gcAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTtc    <  1:1075409/150‑1 (MQ=255)

                                                                                                                                                gcAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGATTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTtc    <  2:494555/150‑1 (MQ=255)

                                                                                                                                                  aaaTTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTtctc  <  1:664703/150‑1 (MQ=255)

CTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCA‑GTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTC  >  NC_000913/3976254‑3976549

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