BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,976,399	C→A	8.1%	A41E (GCA→GAA)	wzxE →	O‑antigen translocase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,976,399	0	C	A	8.1%	436.0 / 21.3	155	A41E (GCA→GAA)	wzxE	O‑antigen translocase
Reads supporting (aligned to +/- strand): ref base C (77/65); new base A (9/4); total (86/69)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

Only 100 of 158 total aligned reads displayed.

GACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCT  >  NC_000913/3976252‑3976548

gACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCaa                                                                                                                                                     >  1:387054/1‑150 (MQ=255)

   tttGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAtt                                                                                                                                                  >  2:8977/1‑150 (MQ=255)

    ttgttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAttt                                                                                                                                                 >  2:1463896/1‑150 (MQ=255)

    ttgttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAttt                                                                                                                                                 <  2:501391/150‑1 (MQ=255)

        tgAACTAGTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTCGTCAAGATTGGTGCAGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGc                                                                                                                                             >  1:414719/1‑150 (MQ=255)

         gAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGcc                                                                                                                                            <  1:586/150‑1 (MQ=255)

          aaCTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCa                                                                                                                                           >  2:1788931/1‑150 (MQ=255)

           actactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAg                                                                                                                                          >  2:18941/1‑150 (MQ=255)

            gctctTTTCCTGAAATGTCGCTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCAgc                                                                                                                                         <  2:69229/147‑2 (MQ=255)

            ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCAgc                                                                                                                                         <  1:43455/150‑2 (MQ=255)

            ctacaTTTCCTCCTATGTCGTTGGCAAAAGCGGCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAATTTGCCGCCAgc                                                                                                                                         <  2:2015342/150‑2 (MQ=255)

              actTTTCCTGATAGCTCGTTGGCACACGCGTCCGTGTTGACGGCGGCCTGTACTATGGTCAAGTCTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGGGGTGTCATTTGGGCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTg                                                                                                                                       <  2:418650/150‑1 (MQ=255)

               ctTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCTGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAATTTTGCCGCCAgcaga                                                                                                                                      >  1:69229/1‑146 (MQ=255)

               ctTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGa                                                                                                                                      >  1:1237997/1‑150 (MQ=255)

                  tcccTGATATGGCGGTGGCACACGCGTCCATGTGGACGGCGGCCAGTACACTGGGCAAGAATGGTGCCGGGTTACTGGTCGGTAAGTGGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTa                                                                                                                                   <  2:153646/148‑1 (MQ=255)

                    ccTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAcc                                                                                                                                 >  2:1853187/1‑150 (MQ=255)

                     cTGATATGTCGTTGGCAAAAGCGGCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGATGATTACCg                                                                                                                                <  2:1280679/150‑1 (MQ=255)

                       gATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTg                                                                                                                              >  1:2093491/1‑150 (MQ=255)

                         tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTAGACGGGTCAAGATTGGTGCCCGGGTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGcg                                                                                                                            >  2:318848/1‑149 (MQ=255)

                         tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCGGGTGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGATGATTACCGTGCt                                                                                                                            >  1:503227/1‑150 (MQ=255)

                            gTCGTTGGCAAAAGCGGCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCgg                                                                                                                         <  1:944737/150‑1 (MQ=255)

                              cGTTGGCAAACGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCg                                                                                                                       >  2:2187100/1‑150 (MQ=255)

                              cGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCg                                                                                                                       >  1:651386/1‑150 (MQ=255)

                               gTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGt                                                                                                                      >  2:1001777/1‑150 (MQ=255)

                               gTTGGCAAAAGCGTCCTTGAGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGGTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGt                                                                                                                      <  2:703223/150‑1 (MQ=255)

                                 tGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGGCAAGATTGGTGCCGGGGTACTGGTCGGTAAGGTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGc                                                                                                                    <  1:1667881/150‑1 (MQ=255)

                                   gCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGGGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAAGTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                                                                                  >  2:644933/1‑150 (MQ=255)

                                        aGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGGTACTGGTCGGTAAGTTGCTGGCGGTGGCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgg                                                                                                             <  1:1451269/150‑1 (MQ=255)

                                         gCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                                                                                            <  1:1925771/150‑1 (MQ=255)

                                          cGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGGTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgggg                                                                                                           <  2:2176095/150‑1 (MQ=255)

                                           gTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                                                                                          >  2:1975355/1‑150 (MQ=255)

                                              cTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCAGGGGCTgg                                                                                                       >  1:1541/1‑150 (MQ=255)

                                               ttGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGc                                                                                                      >  2:1249470/1‑150 (MQ=255)

                                                  tgGACGGCGGCCAGTACACTGGTCAAGATTCGAGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGTCAGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATc                                                                                                   >  2:82351/1‑150 (MQ=255)

                                                   ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCt                                                                                                  >  1:1040907/1‑150 (MQ=255)

                                                       ggcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGGTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTaa                                                                                              >  2:1146366/1‑150 (MQ=255)

                                                        gcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAc                                                                                             >  2:1708560/1‑150 (MQ=255)

                                                          ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                                                                                           <  1:988592/150‑1 (MQ=255)

                                                            ccAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                                                                         <  1:1683213/150‑1 (MQ=255)

                                                            ccAGTACACTGGGCCTGAGTGGTGGCGGGGTACTGGTGGGTTAGTTGCTGGCGGGGTCTTTTGGGCCGGCGGGGCTTGGGCTGGCGGCAACATTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtg                                                                                         <  1:88102/150‑1 (MQ=255)

                                                              aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                                                                                       <  1:1024074/150‑1 (MQ=255)

                                                                  cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCa                                                                                   >  2:1331844/1‑150 (MQ=255)

                                                                   acTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaa                                                                                  <  2:1115572/150‑1 (MQ=255)

                                                                   acTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaa                                                                                  >  1:619362/1‑150 (MQ=255)

                                                                      ggTCAAGAGTGGTGCCGGGTTACTGGGCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                                                                               <  1:1817343/150‑1 (MQ=255)

                                                                         cAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt                                                                            >  1:1715655/1‑150 (MQ=255)

                                                                         cAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGt                                                                            >  1:2203781/1‑150 (MQ=255)

                                                                            gATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc                                                                         <  1:1518026/150‑1 (MQ=255)

                                                                            gATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc                                                                         <  2:1013678/150‑1 (MQ=255)

                                                                              ttGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGccc                                                                       <  1:1853187/150‑1 (MQ=255)

                                                                                ggTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAg                                                                     <  2:53330/150‑1 (MQ=255)

                                                                                ggTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAATACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCat                                                                     >  2:546639/1‑149 (MQ=255)

                                                                                   gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc                                                                  <  2:1529155/150‑1 (MQ=255)

                                                                                   gCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc                                                                  >  1:2164434/1‑150 (MQ=255)

                                                                                    ccGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAcc                                                                 <  1:1272936/150‑1 (MQ=255)

                                                                                    ccGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAcc                                                                 <  1:73168/150‑1 (MQ=255)

                                                                                        gTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGa                                                             <  2:750832/150‑1 (MQ=255)

                                                                                              ggTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc                                                       >  2:1241599/1‑150 (MQ=255)

                                                                                              ggTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc                                                       >  2:1299279/1‑150 (MQ=255)

                                                                                              ggTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc                                                       >  2:1820000/1‑150 (MQ=255)

                                                                                              ggTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAAGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc                                                       >  1:1100317/1‑150 (MQ=255)

                                                                                              ggTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGAAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCAAGTACCATGATAATcc                                                       >  1:374981/1‑150 (MQ=255)

                                                                                                tCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGc                                                     <  2:1169977/150‑1 (MQ=255)

                                                                                                tCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGc                                                     <  1:1755615/150‑1 (MQ=255)

                                                                                                  ggTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCaa                                                   >  1:696694/1‑150 (MQ=255)

                                                                                                   gTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTGGCCCAGGACCATGATTATCCGCAAc                                                  >  1:896217/1‑150 (MQ=255)

                                                                                                       gTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCACCAGCt                                              >  1:74795/1‑150 (MQ=255)

                                                                                                         tGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTgc                                            <  2:2201636/150‑1 (MQ=255)

                                                                                                         tGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTgc                                            >  1:1636066/1‑150 (MQ=255)

                                                                                                             ggCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGcgccg                                        >  2:1896389/1‑150 (MQ=255)

                                                                                                                ggTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGt                                     >  1:181857/1‑150 (MQ=255)

                                                                                                                  tgtCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTgg                                   >  1:544533/1‑150 (MQ=255)

                                                                                                                   gtCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGt                                  <  1:1445012/150‑1 (MQ=255)

                                                                                                                     cATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCg                                >  1:2042472/1‑150 (MQ=255)

                                                                                                                      aTTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGGTTATTACCGTGCTCGGCGTTCTTGCCGGGGCTGGCCTCTTTTACGGTGTAACCAAATACGTTGCCCAGGACCATGATAATCCGCCACAGCTGCGCCGCGTGGTggg                               >  2:1963464/1‑150 (MQ=255)

                                                                                                                       tttGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGc                              >  1:742916/1‑150 (MQ=255)

                                                                                                                       tttGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGc                              <  1:1873289/150‑1 (MQ=255)

                                                                                                                       tttGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCGCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGc                              <  1:1738528/150‑1 (MQ=255)

                                                                                                                       tgcggTCCGTCGGGGCTGGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGGTCGGCGGTCTTGCCGGGGCTGGGATCTTTAACGGTTTAACCAAATAGGTTGCCCACTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCggg                              <  2:1094523/147‑2 (MQ=255)

                                                                                                                         tGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCATGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCAc                            >  1:236712/1‑150 (MQ=255)

                                                                                                                          ggTCCGGCGGGGCTTGGGCTGGCGGCAGATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTATCCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACt                           >  1:1781479/1‑150 (MQ=255)

                                                                                                                          ggTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACt                           >  2:178841/1‑150 (MQ=255)

                                                                                                                          ggTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACt                           >  1:544984/1‑150 (MQ=255)

                                                                                                                           gTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACtt                          >  1:2126709/1‑150 (MQ=255)

                                                                                                                           gTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACtt                          >  1:1386635/1‑150 (MQ=255)

                                                                                                                             ccGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTtca                        <  2:387054/150‑1 (MQ=255)

                                                                                                                             ccGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTtca                        <  1:904228/150‑1 (MQ=255)

                                                                                                                                 cggGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAg                    <  2:745959/150‑1 (MQ=255)

                                                                                                                                 cggGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAg                    >  2:779732/1‑150 (MQ=255)

                                                                                                                                    ggCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGa                 >  2:906147/1‑150 (MQ=255)

                                                                                                                                      cTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATg               <  1:1107529/150‑1 (MQ=255)

                                                                                                                                      cTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATg               >  2:823934/1‑150 (MQ=255)

                                                                                                                                       tttggcgggcgGCAAATTTCCTCCTGTTGATTACCGTGCTCGGCGTGCTTGGCGGGGCTGGGATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCGGCATCAGCTGCGGCGCGTGGTCGGCACTTCATCAGCGATgg              <  2:897585/147‑1 (MQ=255)

                                                                                                                                         caacTGGCGGAAATTTGCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTa            >  2:1082149/4‑150 (MQ=255)

                                                                                                                                          aacTGGCGGAAATTTGCCGCCAGTTTATTACCGTGCTCGGCGTGCTTGCCGTGCCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAACATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTAc           >  2:1353681/3‑150 (MQ=255)

                                                                                                                                          aacTGGCGGAAATTTGCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTAc           >  2:645091/3‑150 (MQ=255)

                                                                                                                                           acTGGCGGAAATTTGCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATTATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCAGCAGCGATGGTACt          >  1:1633185/2‑150 (MQ=255)

                                                                                                                                           acTGGCGGAAATTTGCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAAGACGTTGCCCAGAACCATGATACTCCTCAACAGCTGCGCCGCGTGGTCGGCACTACATCAGCGATGGTAcg          >  2:2164798/2‑149 (MQ=255)

                                                                                                                                            cTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTGGCCGGGGCTAGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACtt         <  2:1693355/150‑1 (MQ=255)

                                                                                                                                                  gcAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTtc   >  2:1669132/1‑150 (MQ=255)

GACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCT  >  NC_000913/3976252‑3976548

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