BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	161,947	A→C	5.1%	F29V (TTC→GTC)	ligT ←	2'‑5' RNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	161,947	0	A	C	5.1%	583.1 / 8.9	170	F29V (TTC→GTC)	ligT	2'‑5' RNA ligase
Reads supporting (aligned to +/- strand): ref base A (54/105); new base C (7/4); total (61/109)
Fisher's exact test for biased strand distribution p-value = 5.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

CATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATGTTAACCCTCT  >  NC_000913/161800‑162094

catcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgt                                                                                                                                                   <  1:698735/150‑1 (MQ=255)

 atcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGACgtc                                                                                                                                                  <  2:528694/150‑2 (MQ=255)

  tcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGACgtgt                                                                                                                                                 >  1:2391381/1‑150 (MQ=255)

  tcaAGCGAGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATAATCGGCGGCGACCGGACGTCACGCCTCAGGTGGGAAgtgt                                                                                                                                                 >  1:2088374/1‑150 (MQ=255)

    aaGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgtgt                                                                                                                                               <  2:2526722/150‑1 (MQ=255)

    aaGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgtgt                                                                                                                                               <  2:1386635/150‑1 (MQ=255)

    aaGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgtgt                                                                                                                                               <  1:1836631/150‑1 (MQ=255)

     aGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGtgtgtg                                                                                                                                              <  2:706251/150‑1 (MQ=255)

       cGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGgc                                                                                                                                            >  2:570564/1‑150 (MQ=255)

            gtgtGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGcc                                                                                                                                       <  2:1288151/150‑1 (MQ=255)

            gtgtGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGcc                                                                                                                                       <  2:1339470/150‑1 (MQ=255)

             tgtgAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGACGGGTGTGGCGCGCCa                                                                                                                                      >  2:2164010/1‑150 (MQ=255)

             tgtgAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCa                                                                                                                                      >  2:86255/1‑150 (MQ=255)

               tgAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAAt                                                                                                                                    >  2:1511439/1‑150 (MQ=255)

                gAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATg                                                                                                                                   <  1:1359635/150‑1 (MQ=255)

                gAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATg                                                                                                                                   <  1:712103/150‑1 (MQ=255)

                gAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATg                                                                                                                                   <  1:2247899/150‑1 (MQ=255)

                gAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATg                                                                                                                                   <  1:997375/150‑1 (MQ=255)

                         ttGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTg                                                                                                                          <  1:393925/150‑1 (MQ=255)

                         ttGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTg                                                                                                                          <  1:1902427/150‑1 (MQ=255)

                           gACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGtt                                                                                                                        <  1:1895198/150‑1 (MQ=255)

                             cGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTcg                                                                                                                      >  2:1512993/1‑150 (MQ=255)

                             cGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTcg                                                                                                                      <  2:1863958/150‑1 (MQ=255)

                             cGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTcg                                                                                                                      <  2:1864031/150‑1 (MQ=255)

                             cGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTcg                                                                                                                      >  1:561408/1‑150 (MQ=255)

                              gAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTcgc                                                                                                                     <  1:1054015/150‑1 (MQ=255)

                               aaTCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTcgcg                                                                                                                    <  1:895535/150‑1 (MQ=255)

                                   cgtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAtt                                                                                                                <  2:652458/150‑1 (MQ=255)

                                   cgtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAtt                                                                                                                <  2:356872/150‑1 (MQ=255)

                                    gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGACGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAttt                                                                                                               >  1:1729179/1‑150 (MQ=255)

                                          cTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAg                                                                                                         <  1:93403/150‑1 (MQ=255)

                                          cTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAg                                                                                                         <  1:93366/150‑1 (MQ=255)

                                           tAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAgg                                                                                                        >  2:320034/1‑150 (MQ=255)

                                            aaaaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGt                                                                                                       <  2:1968812/150‑1 (MQ=255)

                                              aaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTaa                                                                                                     <  1:2548843/150‑1 (MQ=255)

                                              aaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTaa                                                                                                     <  2:2517304/150‑1 (MQ=255)

                                              aaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTaa                                                                                                     <  2:2559118/150‑1 (MQ=255)

                                              aaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTaa                                                                                                     <  2:1446562/150‑1 (MQ=255)

                                              aaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTaa                                                                                                     <  2:329157/150‑1 (MQ=255)

                                              aaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTaa                                                                                                     <  2:1579887/150‑1 (MQ=255)

                                              aaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTaa                                                                                                     >  1:1093443/1‑150 (MQ=255)

                                                    aaGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGAGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGAt                                                                                               >  1:662519/1‑150 (MQ=255)

                                                     aGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATa                                                                                              <  1:1639189/150‑1 (MQ=255)

                                                     aGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATa                                                                                              <  1:1882662/150‑1 (MQ=255)

                                                          cTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCaaa                                                                                         >  1:1677659/1‑150 (MQ=255)

                                                             ctctTTCTCGGCGCTCACTTCGCCTCAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGaa                                                                                      <  1:1940330/150‑1 (MQ=255)

                                                                tttCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAg                                                                                   >  1:2328354/1‑150 (MQ=255)

                                                                tttCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAg                                                                                   >  2:781970/1‑150 (MQ=255)

                                                                 ttCTCGGCGCTCACTTCGCCTAAAAATGCCTGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGa                                                                                  <  1:2029232/150‑1 (MQ=255)

                                                                 ttCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGa                                                                                  <  1:2024923/150‑1 (MQ=255)

                                                                    tcGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGAAGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGt                                                                               >  1:1555154/1‑150 (MQ=255)

                                                                     cGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGtt                                                                              >  2:1194806/1‑150 (MQ=255)

                                                                      ggCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAATTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTg                                                                             <  2:364958/150‑1 (MQ=255)

                                                                       gcgcTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGc                                                                            <  1:2457635/150‑1 (MQ=255)

                                                                       gcgcTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGc                                                                            <  1:2222869/150‑1 (MQ=255)

                                                                        cgcTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCg                                                                           <  2:945297/150‑1 (MQ=255)

                                                                        cgcTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCg                                                                           <  1:2300227/150‑1 (MQ=255)

                                                                           tCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGtt                                                                        >  1:1163732/1‑150 (MQ=255)

                                                                            cACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGACGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTc                                                                       >  1:1629940/1‑150 (MQ=255)

                                                                               ttCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGa                                                                    >  2:2249932/1‑150 (MQ=255)

                                                                                  gCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAt                                                                 <  2:1265695/150‑1 (MQ=255)

                                                                                  gCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAt                                                                 <  2:1962859/150‑1 (MQ=255)

                                                                                  gCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAt                                                                 <  2:1784459/150‑1 (MQ=255)

                                                                                  gCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAt                                                                 <  2:616420/150‑1 (MQ=255)

                                                                                      aaaaaTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGa                                                             <  1:1257831/150‑1 (MQ=255)

                                                                                      aaaaaTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGa                                                             <  1:1500397/150‑1 (MQ=255)

                                                                                        aaaTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAgg                                                           <  2:2342384/150‑1 (MQ=255)

                                                                                        aaaTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAgg                                                           >  1:316817/1‑150 (MQ=255)

                                                                                        aaaTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAgg                                                           >  1:316937/1‑150 (MQ=255)

                                                                                          aTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCa                                                         <  1:375164/150‑1 (MQ=255)

                                                                                               agagTCAGATGCATATTATCGGCGGCGACCGGAAGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGTTAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCgg                                                    <  2:1171233/150‑1 (MQ=255)

                                                                                               agagTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCgg                                                    >  1:758305/1‑150 (MQ=255)

                                                                                                   tCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTAt                                                >  1:281316/1‑150 (MQ=255)

                                                                                                    cAGATGCAAAGTATCGGCGGCGACCGGACGTCCCGCCTCAAGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATg                                               <  2:1933254/150‑1 (MQ=255)

                                                                                                     aGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGa                                              <  1:1512993/150‑1 (MQ=255)

                                                                                                     aGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAAGTGAGGCACTCGGTTATGa                                              >  1:2423188/1‑150 (MQ=255)

                                                                                                        tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAtt                                           >  1:1454519/1‑150 (MQ=255)

                                                                                                           aaaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAcc                                        >  1:566191/1‑150 (MQ=255)

                                                                                                             aTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACcgc                                      <  1:476293/150‑1 (MQ=255)

                                                                                                             aTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACcgc                                      <  1:674978/150‑1 (MQ=255)

                                                                                                                aTCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGAt                                   <  2:612846/150‑1 (MQ=255)

                                                                                                                    gcggcgACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTa                               <  2:1788180/150‑1 (MQ=255)

                                                                                                                    gcggcgACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTa                               <  2:1700909/150‑1 (MQ=255)

                                                                                                                    gcggcgACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTa                               >  1:135170/1‑150 (MQ=255)

                                                                                                                    gcggcgACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTa                               >  1:1008865/1‑150 (MQ=255)

                                                                                                                       gcgACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACaa                            <  1:1211962/150‑1 (MQ=255)

                                                                                                                         gACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAAtg                          >  2:408310/1‑150 (MQ=255)

                                                                                                                           ccGGACGTCCCGCCTCAGGTGGGACGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAAtgtg                        >  2:102220/1‑150 (MQ=255)

                                                                                                                           ccGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAAtgtg                        <  2:1248641/150‑1 (MQ=255)

                                                                                                                           ccGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAAtgtg                        <  1:936422/150‑1 (MQ=255)

                                                                                                                             ggACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGgc                      <  2:2143149/150‑1 (MQ=255)

                                                                                                                             ggACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGgc                      <  1:658614/150‑1 (MQ=255)

                                                                                                                              gACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGgcg                     >  1:416592/1‑150 (MQ=255)

                                                                                                                                  tCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaa                 >  1:2558301/1‑150 (MQ=255)

                                                                                                                                     cGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa              <  1:1460698/150‑1 (MQ=255)

                                                                                                                                      gCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAAt             >  1:1487219/1‑150 (MQ=255)

                                                                                                                                      gCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAAt             >  1:897359/1‑150 (MQ=255)

                                                                                                                                       ccTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATg            >  2:1710679/1‑150 (MQ=255)

                                                                                                                                       ccTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATg            <  1:447883/150‑1 (MQ=255)

                                                                                                                                                 gAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATGTTAACCctct  <  2:2577829/150‑1 (MQ=255)

CATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATGTTAACCCTCT  >  NC_000913/161800‑162094

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