BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,852,039	G→T	7.6%	Q51H (CAG→CAT) ‡	pncA →	nicotinamidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,852,039	0	G	T	7.6%	489.2 / 16.1	147	Q51H (CAG→CAT) ‡	pncA	nicotinamidase
Reads supporting (aligned to +/- strand): ref base G (81/54); new base T (4/8); total (85/62)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 151 total aligned reads displayed.

CCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGC  >  NC_000913/1851890‑1852187

cccccTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAg                                                                                                                                                      <  1:1518599/150‑1 (MQ=255)

cccccTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAg                                                                                                                                                      <  1:1261560/150‑1 (MQ=255)

  cccTCGCGCCCTGTTACTGGGCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGa                                                                                                                                                    <  2:1207729/150‑1 (MQ=255)

       gcgcCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGc                                                                                                                                               <  1:968146/150‑1 (MQ=255)

       gcgcCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGc                                                                                                                                               <  2:2481828/150‑1 (MQ=255)

         gcCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcga                                                                                                                                             <  2:970936/150‑3 (MQ=255)

          cccTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcgat                                                                                                                                            <  1:1422743/150‑4 (MQ=255)

          cccTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCAcc                                                                                                                                            >  2:1822081/1‑150 (MQ=255)

           ccTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCAccc                                                                                                                                           <  2:1813465/150‑1 (MQ=255)

             tGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcgataac                                                                                                                                         <  2:703030/150‑7 (MQ=255)

              gTTCCTGGGCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCGCGCCGGGCCGGAAGGGGACAGTACGGTGGATGGCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGc                                                                                                                                        <  1:2423933/150‑1 (MQ=255)

              gTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGTTTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcgataacc                                                                                                                                        <  2:729702/150‑8 (MQ=255)

               ttACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCg                                                                                                                                       <  2:1988025/150‑1 (MQ=255)

                tACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGa                                                                                                                                      >  1:556153/1‑150 (MQ=255)

                  cTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAAt                                                                                                                                    >  1:46762/1‑150 (MQ=255)

                  cTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAAt                                                                                                                                    >  1:1784514/1‑150 (MQ=255)

                  cTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGATGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAAt                                                                                                                                    >  1:1565259/1‑150 (MQ=255)

                    ggTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCa                                                                                                                                  >  1:2102940/1‑150 (MQ=255)

                    ggTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCa                                                                                                                                  >  1:623420/1‑150 (MQ=255)

                       cGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACgg                                                                                                                               <  1:2361014/150‑1 (MQ=255)

                        gATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGc                                                                                                                              >  2:2109904/1‑150 (MQ=255)

                         ttttACAAAATGATTTCTGTGCTGGGGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGTCGAATCACGGCa                                                                                                                             <  1:1727854/149‑1 (MQ=255)

                          tttACAAAATGATTTCTGTGCTGGGGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAg                                                                                                                            <  2:1173446/150‑1 (MQ=255)

                            tACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGGCTGGCACCCGGCGAATCACGGCAGtt                                                                                                                          >  2:1321858/1‑150 (MQ=255)

                            tACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGtt                                                                                                                          >  2:1626005/1‑150 (MQ=255)

                              cAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGtttt                                                                                                                        >  2:2264355/1‑150 (MQ=255)

                               aaaaTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTg                                                                                                                       <  2:1735745/150‑1 (MQ=255)

                                aaaTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGc                                                                                                                      >  1:1547827/1‑150 (MQ=255)

                                aaaTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGc                                                                                                                      >  1:818484/1‑150 (MQ=255)

                                aaaTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAATCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGc                                                                                                                      >  2:773561/1‑150 (MQ=255)

                                aaaTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAATCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGc                                                                                                                      >  2:765178/1‑150 (MQ=255)

                                 aaTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGcc                                                                                                                     <  1:1093477/150‑1 (MQ=255)

                                  aTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCa                                                                                                                    >  2:623373/1‑150 (MQ=255)

                                   tGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAg                                                                                                                   <  1:19487/150‑1 (MQ=255)

                                     aTTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTc                                                                                                                 >  2:2426293/1‑150 (MQ=255)

                                      tttCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCa                                                                                                                >  2:1279857/1‑150 (MQ=255)

                                                ggtggCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAg                                                                                                               >  2:129797/1‑141 (MQ=255)

                                                ggtggCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAg                                                                                                               <  1:129797/141‑1 (MQ=255)

                                                 gtggCGCGCTCGCCGTGCTGTTATGTGAGAGTTCGGTGGATGTCGCTAACCTCCTGATTGTGTTGTTCCATTCGCGCGGTGAAGCGGTTATCGACAGTCAGGTCTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTag                                                                                                     <  2:842994/150‑1 (MQ=255)

                                                     cgcgcTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGcc                                                                                                 >  1:678378/1‑150 (MQ=255)

                                                              cGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGcc                                                                                        <  2:1649213/150‑1 (MQ=255)

                                                                tGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAg                                                                                      >  2:2144196/1‑150 (MQ=255)

                                                                 gCCGGAAGGTGCCAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAgg                                                                                     >  1:165499/1‑150 (MQ=255)

                                                                    ggAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCa                                                                                  >  1:2023273/1‑150 (MQ=255)

                                                                        ggTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTc                                                                              >  2:355909/1‑150 (MQ=255)

                                                                         gTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCg                                                                             >  1:441758/1‑150 (MQ=255)

                                                                         gTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCg                                                                             <  2:2264942/150‑1 (MQ=255)

                                                                           gACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGAc                                                                           >  1:567529/1‑150 (MQ=255)

                                                                               gTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGtt                                                                       >  2:321128/1‑150 (MQ=255)

                                                                               gTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGtt                                                                       <  1:2234115/150‑1 (MQ=255)

                                                                                   ggTGGCTGTCGCTAACCGCCTGATTGACTGGTGCCCGTCGCGCGGTGACTCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGGTTTGCCAGTCAGCACGGTGTAGAGCCTTATAGGCCAGGCCAACTCGACGGTTTGcc                                                                   <  2:2109616/150‑1 (MQ=255)

                                                                                      ggATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCcaca                                                                >  2:520229/1‑150 (MQ=255)

                                                                                        ttGTCGCTAACCGCCTGATTGACTGGTGCGAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGAGTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACaaa                                                              <  2:1042475/149‑1 (MQ=255)

                                                                                         tGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAAc                                                             <  2:386507/150‑1 (MQ=255)

                                                                                         tGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAAc                                                             >  1:2249610/1‑150 (MQ=255)

                                                                                          gTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAAcc                                                            >  2:1235961/1‑150 (MQ=255)

                                                                                           tCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCt                                                           >  1:1396352/1‑150 (MQ=255)

                                                                                               tAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTg                                                       <  2:332570/150‑1 (MQ=255)

                                                                                               tAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTg                                                       >  2:44850/1‑150 (MQ=255)

                                                                                                aaCCGGCTGCATGACCGGTGGCACGCGCGGGGGGAACGGGTTATCGCCATTCATGTACGGCACCCGGCGAAACACGGGAGTTTTTGCACTCATCAAGGTGTATATGCTTATACGCCACGCCAACTCGACGGTTTGCCACAAACCTTCTgg                                                      <  2:2070472/150‑1 (MQ=255)

                                                                                                aaCCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTgg                                                      >  1:1654093/1‑150 (MQ=255)

                                                                                                aaCCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTgg                                                      <  2:640016/150‑1 (MQ=255)

                                                                                                  ccGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGcc                                                    >  1:696347/1‑150 (MQ=255)

                                                                                                        gATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGTACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCa                                              >  2:431236/1‑150 (MQ=255)

                                                                                                            gACTGGTTCAAGTCGCGCGGTAAAGAGTTTATCGCCAGTCATCACTGTCACCCGGCGAATAAAGACAATTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtgt                                          >  2:621115/1‑150 (MQ=255)

                                                                                                            gACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtgt                                          >  2:1761125/1‑150 (MQ=255)

                                                                                                            gACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtgt                                          >  1:2413397/1‑150 (MQ=255)

                                                                                                            gACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtgt                                          >  1:1691679/1‑150 (MQ=255)

                                                                                                             aCTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtgtg                                         >  1:2124729/1‑150 (MQ=255)

                                                                                                             aCTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACCGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtgtg                                         <  1:2012557/150‑1 (MQ=255)

                                                                                                                ggTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGc                                      >  2:831511/1‑150 (MQ=255)

                                                                                                                ggTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGc                                      <  1:520229/150‑1 (MQ=255)

                                                                                                                ggTGCCAATCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGc                                      <  1:765178/150‑1 (MQ=255)

                                                                                                                ggTGCCAATCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGc                                      <  1:773561/150‑1 (MQ=255)

                                                                                                                  tGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAg                                    <  2:217514/150‑1 (MQ=255)

                                                                                                                     cAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAAc                                 >  1:1027710/1‑150 (MQ=255)

                                                                                                                      aGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACa                                <  1:102834/150‑1 (MQ=255)

                                                                                                                         cgcgcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTg                             >  1:98845/1‑150 (MQ=255)

                                                                                                                         cgcgcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTg                             >  1:2339877/1‑150 (MQ=255)

                                                                                                                            gcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAg                          >  2:72123/1‑150 (MQ=255)

                                                                                                                            gcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAg                          >  1:1891029/1‑150 (MQ=255)

                                                                                                                               gTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGgcg                       >  1:1123305/1‑150 (MQ=255)

                                                                                                                                   aGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGcaca                   >  2:210172/1‑150 (MQ=255)

                                                                                                                                   aGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGAGcaca                   >  2:1510312/1‑150 (MQ=255)

                                                                                                                                   aGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGAGcaca                   >  2:1519479/1‑150 (MQ=255)

                                                                                                                                      ggTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCACCTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAGCAGTGAAGGCGCACAAtt                >  2:4823/1‑150 (MQ=255)

                                                                                                                                      ggTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCTCAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAAtt                <  2:285296/150‑1 (MQ=255)

                                                                                                                                      ggTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAAtt                >  1:1793295/1‑150 (MQ=255)

                                                                                                                                       gTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTa               <  1:2109904/150‑1 (MQ=255)

                                                                                                                                       gTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTa               <  2:1977065/150‑1 (MQ=255)

                                                                                                                                       gTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTa               <  1:2135537/150‑1 (MQ=255)

                                                                                                                                       gTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTa               <  1:1411268/150‑1 (MQ=255)

                                                                                                                                        ttATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTAc              >  2:2335716/1‑150 (MQ=255)

                                                                                                                                         tATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACa             >  1:1368212/1‑150 (MQ=255)

                                                                                                                                         gtgtcccAGTCCTGACTGGCACCCGGCGAATCACGGCAGATTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACa             >  2:304148/6‑150 (MQ=255)

                                                                                                                                         gggtgCCAGTCCTGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACa             >  2:2470835/5‑150 (MQ=255)

                                                                                                                                           tCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATc           >  1:925538/1‑150 (MQ=255)

                                                                                                                                             gCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGCGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCg         >  2:1067085/1‑150 (MQ=255)

                                                                                                                                                   cAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTg   >  1:1456376/1‑150 (MQ=255)

                                                                                                                                                    aGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGc  <  2:1435053/150‑1 (MQ=255)

CCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGC  >  NC_000913/1851890‑1852187

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