BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,624,896	A→T	7.3%	T42S (ACC→TCC)	ydeJ →	conserved protein YdeJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,624,896	0	A	T	7.3%	362.8 / 13.3	122	T42S (ACC→TCC)	ydeJ	conserved protein YdeJ
Reads supporting (aligned to +/- strand): ref base A (50/63); new base T (2/7); total (52/70)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 124 total aligned reads displayed.

ACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGA  >  NC_000913/1624747‑1625044

aCATATAAAAGCGTGGACAGAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCa                                                                                                                                                      >  2:473472/1‑150 (MQ=255)

    ataAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACcgg                                                                                                                                                  >  2:279607/1‑150 (MQ=255)

     taAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACcggc                                                                                                                                                 <  2:940814/150‑1 (MQ=255)

      aaaaGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACcgtcg                                                                                                                                                >  2:1204122/1‑150 (MQ=255)

      aaaaGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACcggcg                                                                                                                                                <  2:485190/150‑1 (MQ=255)

         aGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGGTGCATGATTCGGCGGTa                                                                                                                                             <  2:685005/150‑1 (MQ=255)

           cGTGTACAAAAAACAGTGAACATTAATAGATATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTTAGTCAATGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAg                                                                                                                                           <  2:1341365/150‑1 (MQ=255)

           cGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTTTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGGTGCATGATTCGGCGGTagt                                                                                                                                           <  2:200534/150‑3 (MQ=255)

           cGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCTATTCATGCACCGGCGGTAAg                                                                                                                                           <  2:716128/150‑1 (MQ=255)

           cGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGGTGCATGATTCGGCGGTagt                                                                                                                                           <  1:909445/150‑3 (MQ=255)

              ggACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTg                                                                                                                                        >  1:1125320/1‑150 (MQ=255)

               gACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGGTGCATGATTCGGCGGTagttaag                                                                                                                                       <  2:594239/150‑7 (MQ=255)

               gACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGGTGCATGATTCGGCGGTagttaag                                                                                                                                       <  2:737309/150‑7 (MQ=255)

               gACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTgg                                                                                                                                       <  2:1491155/150‑1 (MQ=255)

               gACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCCTAACTACCGCCGAATCATGCACCGGCGGTAAGTTgg                                                                                                                                       >  2:56949/1‑150 (MQ=255)

                aCAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGc                                                                                                                                      <  2:381976/150‑1 (MQ=255)

                aCAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGc                                                                                                                                      >  2:1027939/1‑150 (MQ=255)

                 cAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCt                                                                                                                                     >  2:874948/1‑150 (MQ=255)

                 cAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCt                                                                                                                                     >  1:87810/1‑150 (MQ=255)

                 cAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCt                                                                                                                                     >  1:259074/1‑150 (MQ=255)

                 cAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGAATTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCt                                                                                                                                     >  1:1159548/1‑150 (MQ=255)

                        cAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGccc                                                                                                                              <  1:456446/150‑1 (MQ=255)

                          gTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtg                                                                                                                            <  1:25725/150‑1 (MQ=255)

                          gTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATAACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtg                                                                                                                            >  2:823748/1‑150 (MQ=255)

                           tGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtgt                                                                                                                           >  1:993920/1‑150 (MQ=255)

                            gAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtgtg                                                                                                                          >  2:1269056/1‑150 (MQ=255)

                             aaCATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtgtgt                                                                                                                         >  1:716538/1‑150 (MQ=255)

                              aCATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCtgtgtg                                                                                                                        <  2:1354158/150‑1 (MQ=255)

                                 ttAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAg                                                                                                                     <  2:1507574/150‑1 (MQ=255)

                                  tAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGc                                                                                                                    <  2:399776/150‑1 (MQ=255)

                                   aaTAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCt                                                                                                                   >  1:572690/1‑150 (MQ=255)

                                    aTAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTg                                                                                                                  >  1:817570/1‑150 (MQ=255)

                                     tAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGa                                                                                                                 <  2:297984/150‑1 (MQ=255)

                                      agagATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGaa                                                                                                                <  2:469419/150‑1 (MQ=255)

                                         taTAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGAt                                                                                                             <  1:260128/149‑1 (MQ=255)

                                         gaTAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGAt                                                                                                             <  2:1280821/150‑1 (MQ=255)

                                             aaaTTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGGTGCATGATTCGGCGGTagtt                                                                                                                                          <  2:144133/117‑4 (MQ=38)

                                             aaaTTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGGTGCATGATTCGGCGGTagtt                                                                                                                                          >  1:144133/1‑114 (MQ=38)

                                              aaTTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACAcc                                                                                                        <  2:259074/150‑1 (MQ=255)

                                                   tACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACAGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAAt                                                                                                   >  2:878804/1‑150 (MQ=255)

                                                    aCAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAAtt                                                                                                  >  2:451641/1‑150 (MQ=255)

                                                     cAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAAttt                                                                                                 >  1:1497195/1‑150 (MQ=255)

                                                        tttGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTAc                                                                                              <  2:537025/149‑1 (MQ=255)

                                                         ttGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACg                                                                                             <  2:1412952/150‑1 (MQ=255)

                                                         ttGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACg                                                                                             <  2:732312/150‑1 (MQ=255)

                                                         ttGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACg                                                                                             <  1:390830/150‑1 (MQ=255)

                                                           gTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGt                                                                                           >  1:616008/1‑150 (MQ=255)

                                                            tAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTg                                                                                          >  2:969978/1‑150 (MQ=255)

                                                                aaCCGATACTATTCAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATTCACCGGCGGTAAGTTGGCTAGCGCCCTTTGTGAAGCGGAAGATACTCCCAAATTTTACGGTGCAgg                                                                                      <  2:1305108/149‑1 (MQ=255)

                                                                    gATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCttt                                                                                  >  1:1308356/1‑150 (MQ=255)

                                                                        cTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTa                                                                              >  1:1419071/1‑150 (MQ=255)

                                                                        cTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTa                                                                              >  1:361429/1‑150 (MQ=255)

                                                                          aTTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACt                                                                            <  2:1308356/150‑1 (MQ=255)

                                                                           ttGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACtt                                                                           <  2:616008/150‑1 (MQ=255)

                                                                           ttGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACtt                                                                           >  2:227556/1‑150 (MQ=255)

                                                                             gAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTc                                                                         >  1:974822/1‑150 (MQ=255)

                                                                             gAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAAACATGCACCGGCGGTAACTTGGCTAGCGAACTGTGTGCAGCTGAAGATACACCCAAATTATACGGTGAAGGCTATGTTACTTTc                                                                         >  1:1026543/1‑150 (MQ=255)

                                                                              aaaaCCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCa                                                                        <  1:36412/150‑1 (MQ=255)

                                                                               aaaCCTGACATCCGCGTTGAGTCAAAGACTTATCTCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCAc                                                                       <  1:1204122/150‑1 (MQ=255)

                                                                                   cTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGAt                                                                   >  2:1045441/1‑150 (MQ=255)

                                                                                    tGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATc                                                                  <  1:992077/150‑1 (MQ=255)

                                                                                    tGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATc                                                                  <  2:442653/150‑1 (MQ=255)

                                                                                    ggACATCCGCGGTGAGGCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGCGCAGGCTTTGTTACTTTCACCGATc                                                                  <  2:862577/149‑1 (MQ=255)

                                                                                     ttgaTCCGCGATAAGTCTTTGACTCAACGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCa                                                                 >  2:503387/4‑150 (MQ=255)

                                                                                       cATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAgg                                                               <  1:341031/150‑1 (MQ=255)

                                                                                       cATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAgg                                                               <  1:969978/150‑1 (MQ=255)

                                                                                        aTCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGc                                                              <  1:1499404/150‑1 (MQ=255)

                                                                                           cgcgTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCaaa                                                           >  2:1369680/1‑150 (MQ=255)

                                                                                                tGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGa                                                      <  2:193274/150‑1 (MQ=255)

                                                                                                tGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGa                                                      >  1:536346/1‑150 (MQ=255)

                                                                                                  aGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGaaa                                                    >  1:1083865/1‑150 (MQ=255)

                                                                                                   gTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCTCCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGaaaa                                                   >  1:633496/1‑150 (MQ=255)

                                                                                                   gTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGaaaa                                                   >  1:1066930/1‑150 (MQ=255)

                                                                                                     cAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATc                                                 >  1:918503/1‑150 (MQ=255)

                                                                                                     cAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATc                                                 >  2:282592/1‑150 (MQ=255)

                                                                                                      aaaGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATcc                                                <  1:451641/150‑1 (MQ=255)

                                                                                                        aGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTc                                              <  2:974822/150‑1 (MQ=255)

                                                                                                        aGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTc                                              <  2:87810/150‑1 (MQ=255)

                                                                                                          aCTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAg                                            <  2:1155566/150‑1 (MQ=255)

                                                                                                           cTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGc                                           <  1:1352142/150‑1 (MQ=255)

                                                                                                            ttATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCg                                          <  1:878804/150‑1 (MQ=255)

                                                                                                             tATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGt                                         >  1:1441632/1‑150 (MQ=255)

                                                                                                                 gcgGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGc                                     <  1:1141768/150‑1 (MQ=255)

                                                                                                                 gcgGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGc                                     <  1:832661/150‑1 (MQ=255)

                                                                                                                 gcgGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGc                                     >  1:1405653/1‑150 (MQ=255)

                                                                                                                    gATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAg                                  >  1:786177/1‑150 (MQ=255)

                                                                                                                     aTCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGc                                 <  1:473472/150‑1 (MQ=255)

                                                                                                                      tCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCa                                <  2:1125320/150‑1 (MQ=255)

                                                                                                                       cAATTACGCTTAACTACCGCCGAATCATGCACCGGCTGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCaa                               >  1:552810/1‑150 (MQ=255)

                                                                                                                             cGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCtt                         <  1:1027939/150‑1 (MQ=255)

                                                                                                                              gCTTAACTACCGCCGAATCTTGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGGGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTg                        <  2:884269/150‑1 (MQ=255)

                                                                                                                                 tAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAAc                     >  2:388992/1‑150 (MQ=255)

                                                                                                                                    cTACCGCCGAATCATGCACAGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTCCGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCACTCTGTTGAACGat                  >  2:772320/1‑150 (MQ=255)

                                                                                                                                     tACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGata                 <  2:918503/150‑1 (MQ=255)

                                                                                                                                      aCCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGatat                >  1:263762/1‑150 (MQ=255)

                                                                                                                                         gccgAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACTGTGCATGCTTTGTTAGTTTCACCTATCATTCAAAGATGAAAATCCTTAGCGTAATCCAGGAATCTCTTTTACGTTATTCt             <  2:1237046/150‑1 (MQ=255)

                                                                                                                                                 aTGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGgtgag     <  1:1045441/150‑1 (MQ=255)

                                                                                                                                                  tGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGgtgagt    <  1:446652/150‑1 (MQ=255)

                                                                                                                                                  tGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGgtgagt    >  2:1426504/1‑150 (MQ=255)

                                                                                                                                                    cACCGGCTGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGgtgagtga  <  2:824674/150‑1 (MQ=255)

ACATATAAAAGCGTGGACAAAAAACAGTGAACATTAATAGAGATAAAATTGTACAACTTGTAGATACCGATACTATTGAAAACCTGACATCCGCGTTGAGTCAAAGACTTATCGCGGATCAATTACGCTTAACTACCGCCGAATCATGCACCGGCGGTAAGTTGGCTAGCGCCCTGTGTGCAGCTGAAGATACACCCAAATTTTACGGTGCAGGCTTTGTTACTTTCACCGATCAGGCAAAGATGAAAATCCTCAGCGTAAGCCAGCAATCTCTTGAACGATATTCTGCGGTGAGTGA  >  NC_000913/1624747‑1625044

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