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 →	O‑antigen translocase

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%	449.4 / 20.7	157	N42I (AAT→ATT)	wzxE	O‑antigen translocase
Reads supporting (aligned to +/- strand): ref base A (80/64); new base T (8/4); total (88/69)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

Only 100 of 162 total aligned reads displayed.

TTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTA  >  NC_000913/3976255‑3976551

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

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

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

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

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

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

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

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

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

             ttttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAt                                                                                                                                        >  2:1738528/1‑150 (MQ=255)

               ttCCTGATATGTCGTTGGCAAACGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTa                                                                                                                                      <  1:1731839/150‑1 (MQ=255)

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

                 ccTGATATCTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAcc                                                                                                                                    >  2:793037/1‑150 (MQ=255)

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

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

                      tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGTGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCt                                                                                                                               >  2:1793182/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)

                           cGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCggtg                                                                                                                          <  2:536035/150‑3 (MQ=255)

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

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

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

                                    aaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGCTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGCTTACCGTGCTCGGCGTGCTTGCCg                                                                                                                 >  2:1215766/1‑150 (MQ=255)

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

                                       cGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgggg                                                                                                              <  2:1460045/150‑1 (MQ=255)

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

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

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

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

                                           attGTGGACGGCGGCCAGTACACTGGGCAAGATGGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTgg                                                                                                          <  2:1325252/149‑1 (MQ=255)

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

                                               tgGACGGCGGACAGTACACTGGTCAAGATAGGTGCCGGGTTACTGGTAGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGATTGCCGGGGCAGGCATc                                                                                                      >  1:1669989/1‑150 (MQ=255)

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

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

                                                      cggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACg                                                                                               >  1:713459/1‑150 (MQ=255)

                                                      cggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGGCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACg                                                                                               <  1:1319204/150‑1 (MQ=255)

                                                       ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                                                                                              <  2:1034649/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                                                                                     >  1:619362/1‑150 (MQ=255)

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

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

                                                                       aaGATTGGTGCCGGGTTACTGGGCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt                                                                              <  2:1703813/150‑1 (MQ=255)

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

                                                                          atgggTGCCGGGTATCTGGTCGGTAAGGTGCGGGCGTGGAGAATTGGTCCGGCGGGGCTTGGGCTGTCTGCAAATTTCCTCCAGTTGATTAGCGGGCTCGGCGTGCATCCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGcc                                                                           <  2:435296/147‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                     gTTACTGGGCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGa                                                                <  1:1811178/150‑1 (MQ=255)

                                                                                        aCTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATaa                                                             >  1:745959/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)

                                                                                                  aaGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAg                                                   <  1:963762/150‑1 (MQ=255)

                                                                                                     ttGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTg                                                >  2:233550/1‑150 (MQ=255)

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

                                                                                                      tGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTgc                                               <  1:8977/150‑1 (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)

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

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

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

                                                                                                                    tttGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGc                                 >  1:742916/1‑150 (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)

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

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

                                                                                                                                gggCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCg                     <  1:41829/150‑1 (MQ=255)

                                                                                                                                 ggCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACGATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGa                    >  1:1407967/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)

                                                                                                                                        gCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt             <  1:793037/150‑1 (MQ=255)

                                                                                                                                        acTGGCGGAAATTTGCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt             >  1:577135/2‑150 (MQ=255)

                                                                                                                                        acTGGCGGAAATTTGCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACt             >  2:1553316/2‑150 (MQ=255)

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

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

                                                                                                                                                cAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTtct     >  1:1459744/1‑150 (MQ=255)

                                                                                                                                                 aaaTTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTtctc    <  2:1273824/150‑1 (MQ=255)

                                                                                                                                                  aaTTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTctct   <  2:1882046/150‑1 (MQ=255)

                                                                                                                                                   aTTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCAAAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGCTTTCTCTa  >  2:672588/1‑150 (MQ=255)

TTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAACAGCTGCGCCGCGTGGTCGGCACTTCATCAGCGATGGTACTTGGTTTCTCTA  >  NC_000913/3976255‑3976551

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