BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,720,793	G→T	5.8%	Q38H (CAG→CAT)	insJ →	IS150 transposase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720,793	0	G	T	5.8%	427.2 / 13.5	155	Q38H (CAG→CAT)	insJ	IS150 transposase A
Reads supporting (aligned to +/- strand): ref base G (64/82); new base T (2/7); total (66/89)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-02

Only 100 of 155 total aligned reads displayed.

TTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTC  >  NC_000913/3720648‑3720941

ttgttgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACACTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCa                                                                                                                                                  >  1:1453925/1‑150 (MQ=255)

 tgttgGACGATGAAATGTAATAGCCCCTAATATGTCAAAGCAAAAATACCCTTTTGAAAAGCTCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAg                                                                                                                                                 <  2:1622252/150‑1 (MQ=255)

 tgttgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTTATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAg                                                                                                                                                 <  2:1089542/150‑1 (MQ=255)

  gttgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAgg                                                                                                                                                >  1:1450366/1‑150 (MQ=255)

    tgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATGACGGCACGTTTTGGTGTCCCCCGAACCCTGGTCAGGAc                                                                                                                                              >  2:165802/1‑150 (MQ=255)

    tgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGAc                                                                                                                                              <  2:952587/150‑1 (MQ=255)

     ggACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACa                                                                                                                                             >  1:565213/1‑150 (MQ=255)

     ggACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATAATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACa                                                                                                                                             >  1:823267/1‑150 (MQ=255)

         gATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATggg                                                                                                                                         <  1:1495316/150‑1 (MQ=255)

          aTGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATCGTTACAGGATCATCTCGGCACGTTTTGTTGTCCCCCGAACCCAGGTCAGGACATGGGt                                                                                                                                        >  2:43684/1‑150 (MQ=255)

           tGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGtt                                                                                                                                       >  2:1080046/1‑150 (MQ=255)

            gAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttc                                                                                                                                      <  2:1649377/150‑2 (MQ=255)

            gAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTg                                                                                                                                      >  2:1207644/1‑150 (MQ=255)

            gAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTg                                                                                                                                      <  2:573105/150‑1 (MQ=255)

             aaaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCTGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGTGTTGc                                                                                                                                     >  2:1722613/1‑150 (MQ=255)

             aaaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcg                                                                                                                                     <  1:1317543/150‑3 (MQ=255)

             aaaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGc                                                                                                                                     >  1:654323/1‑150 (MQ=255)

             aaaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTAACAGGATCAACTCTGAACTTTTTGGTGTCACACGACCCCAGGTCAGCACATTGGTTGc                                                                                                                                     >  2:1357897/1‑150 (MQ=255)

             aaaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACACTTATGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGCTGAAGGTTCCAGTATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCCTGTCCTGACCTGGGttcg                                                                                                                                     <  1:820212/150‑3 (MQ=255)

             aaaTGGAAAAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGc                                                                                                                                     <  2:1981231/150‑1 (MQ=255)

               aTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcggg                                                                                                                                   <  1:906695/150‑5 (MQ=255)

                tGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                                                                  <  1:1288728/150‑1 (MQ=255)

                    aTAGCCCCTAAAATGTCAAAGCCAAAATACCATTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTAt                                                                                                                              <  2:1449110/150‑1 (MQ=255)

                         cccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGGCGTGAATCACTACTTCACAACTGATGATGGGTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaaa                                                                                                                         <  1:803012/150‑1 (MQ=255)

                          ccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaaaa                                                                                                                        >  2:1020792/1‑150 (MQ=255)

                            tAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACa                                                                                                                      >  2:319257/1‑150 (MQ=255)

                             aaTATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACAt                                                                                                                     >  2:61371/1‑150 (MQ=255)

                               tatGTCAAAGCCAAAATACCCTTTTTAAACGCGCCTTGAAGTCGTGAATAACTACATCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCATGACATGGGTTGCCCTCTATGAAAAACATgg                                                                                                                   <  2:1034153/150‑1 (MQ=255)

                                 tGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGgag                                                                                                                 <  2:1148661/150‑1 (MQ=255)

                                  gTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGgaga                                                                                                                >  2:1961680/1‑150 (MQ=255)

                                    cAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGaaa                                                                                                              >  2:804589/1‑150 (MQ=255)

                                    cAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGaaa                                                                                                              >  1:748792/1‑150 (MQ=255)

                                         ccAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGt                                                                                                         >  1:1173201/1‑150 (MQ=255)

                                          cAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGCTCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGtt                                                                                                        >  2:1591237/1‑150 (MQ=255)

                                          cAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGtt                                                                                                        <  2:1935296/150‑1 (MQ=255)

                                          cAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGtt                                                                                                        >  1:745959/1‑150 (MQ=255)

                                          cAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACAACGTCACAACGGATGATGGTTACAGGATCATCTCGGCCCGTTTTGGTGGCCCCCGAACCCAGGTCAGGACAGGGGTTGCCCTCGATGAAAAACATGGAGAAAAAGGtt                                                                                                        >  2:1063959/1‑150 (MQ=255)

                                            aaaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTa                                                                                                      >  1:1474251/1‑150 (MQ=255)

                                                aCCCTTTTGAAAAGCGCCCTGCAGTCGGGAATCACTACTACACAACTGATGATGGTTACAGCATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACCTGGAGAAAAAGGTTTAATTc                                                                                                  <  1:1361976/150‑1 (MQ=255)

                                                 cccTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTcc                                                                                                 >  1:1941398/1‑150 (MQ=255)

                                                  ccTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACACGGAGAAAAAGGTTTAATTccc                                                                                                >  2:285276/1‑150 (MQ=255)

                                                     tttGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCaat                                                                                             <  1:235723/150‑2 (MQ=255)

                                                     tttGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCaaa                                                                                             >  1:1917449/1‑150 (MQ=255)

                                                      ttGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGTCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAc                                                                                            <  1:48586/150‑1 (MQ=255)

                                                       tGAAAAGCGCCTTGCAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTTGTGTCCCCCGAACCGAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAcc                                                                                           <  1:321299/150‑1 (MQ=255)

                                                              cgcCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGGTTTGGGGTCCCCCGCACCCAGGTCAGGACATGGGTTGCGCTCTAGGAAAAACATGGAGAAAAAGGTTTATTTACCAAAGCGAAAggg                                                                                    >  1:325028/1‑149 (MQ=255)

                                                                ccTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGt                                                                                  <  2:1258391/150‑1 (MQ=255)

                                                                 cTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCTTCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGtt                                                                                 <  1:477827/150‑1 (MQ=255)

                                                                    gAAGTCGTGAATCACTACTTACCAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATTGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGt                                                                              >  2:161114/1‑150 (MQ=255)

                                                                     aaGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTg                                                                             >  2:2050453/1‑150 (MQ=255)

                                                                     aaGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTg                                                                             >  1:932450/1‑150 (MQ=255)

                                                                      aGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGc                                                                            <  1:1080047/150‑1 (MQ=255)

                                                                      aGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGc                                                                            <  1:465910/150‑1 (MQ=255)

                                                                      aGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTATGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCGCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGc                                                                            <  1:392229/150‑1 (MQ=255)

                                                                       gTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCt                                                                           >  2:1538873/1‑150 (MQ=255)

                                                                         cGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCACCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGtgatga                                                                         >  1:101751/1‑150 (MQ=255)

                                                                           tGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATc                                                                       >  1:1976982/1‑150 (MQ=255)

                                                                            gAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATcc                                                                      <  1:1468270/150‑1 (MQ=255)

                                                                              aTCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGTGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCag                                                                    >  2:2036765/1‑150 (MQ=255)

                                                                              aTCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCag                                                                    >  1:946205/1‑150 (MQ=255)

                                                                              aTCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCag                                                                    <  2:500021/150‑1 (MQ=255)

                                                                                 actactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGt                                                                 <  1:1670396/150‑1 (MQ=255)

                                                                                  ctactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGtt                                                                >  2:309750/1‑150 (MQ=255)

                                                                                     ctTCACAACTGATGATGGTTACAGGAACATCTCGGCAAGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                             <  1:785740/150‑1 (MQ=255)

                                                                                     attCACAACTGATGATGGTTACATGATCATCTCGGAACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                             <  1:1305676/149‑1 (MQ=255)

                                                                                     attCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                             <  2:280290/149‑1 (MQ=255)

                                                                                          caACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTa                                                        <  2:1981586/150‑1 (MQ=255)

                                                                                              tgatgaTGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGt                                                    <  1:1993619/150‑1 (MQ=255)

                                                                                              tgatgaTGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCAAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGCAAAAGGTTTAATTACCAAACCTACAGGCGTTAGAGCTGATCCAGAGTTGCGAATTAAggg                                                    >  2:867501/1‑149 (MQ=255)

                                                                                                  gatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCATTGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTg                                                >  2:1217314/1‑150 (MQ=255)

                                                                                                  gatgGGTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGGTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTg                                                <  2:1349356/150‑1 (MQ=255)

                                                                                                     ggTTACAGGATCATCTCGGCACG‑TTTGGTTTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGAGTATTAAGGTCGTGAAAg                                            <  1:1612468/150‑1 (MQ=255)

                                                                                                      gTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAg                                            <  1:1020793/150‑1 (MQ=255)

                                                                                                        tACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCt                                          >  1:256424/1‑150 (MQ=255)

                                                                                                          cAGGCTCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCGCTCTATGAAAAACATGGAGAAAATGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGCTCCAGAGTTGCGTATTAAGGTCGTGAAAGCtgt                                        >  1:318625/1‑150 (MQ=255)

                                                                                                          cAGGATCATCTCTGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGTGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCtgt                                        <  1:1722614/150‑1 (MQ=255)

                                                                                                          cAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCtgt                                        <  1:1354182/150‑1 (MQ=255)

                                                                                                           aGGATCATCTCGGCACGTTTTGGTGTCCCCCTAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCtgtg                                       <  1:1609700/150‑1 (MQ=255)

                                                                                                                 atcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCAAAACCTAAAGTCGTTATAGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGag                                 >  2:1055976/1‑150 (MQ=255)

                                                                                                                   ctcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagca                               <  2:1533606/150‑1 (MQ=255)

                                                                                                                     cGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagcagc                             <  2:745958/150‑1 (MQ=255)

                                                                                                                      ggCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGATAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAgcagca                            >  1:57566/1‑150 (MQ=255)

                                                                                                                        cACGTTTTTGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGcaca                          <  2:1736825/150‑1 (MQ=255)

                                                                                                                        cACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGTGTTGCCCTCTATGAAAAACATGGAGAAAAAGGGTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGcaca                          <  1:622064/150‑1 (MQ=255)

                                                                                                                          cGTTTTGGTGTCCCCCGACCCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAATAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATg                        >  1:1212337/1‑150 (MQ=255)

                                                                                                                              ttGGTGTCCCCCGAACCCAGGGCAGGACATGGGGTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTccc                    <  1:1482959/150‑1 (MQ=255)

                                                                                                                               tGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCt                   <  2:654323/150‑1 (MQ=255)

                                                                                                                                   gtCCCCCGAACCCGGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAACCCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAAt               >  2:859412/1‑150 (MQ=255)

                                                                                                                                    tCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATc              <  2:1450365/150‑1 (MQ=255)

                                                                                                                                    tCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCAATGAAAAACATGGAGAAAATGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATc              <  1:411939/150‑1 (MQ=255)

                                                                                                                                      ccccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAg            >  1:1099452/1‑150 (MQ=255)

                                                                                                                                       gggcaaCCCATGTCCTGACCTGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAgg           >  2:612341/5‑150 (MQ=255)

                                                                                                                                       gggcaaCCCATGTCCTGACCTGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTAAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTAGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAgg           >  2:6352/5‑150 (MQ=255)

                                                                                                                                       ccgcaaCCCTGGTCAGTACATGGGTTGCCCTTTATTAAAAAATTGTATAAAAAGGTTTAA‑TACCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgc          <  1:1170918/146‑1 (MQ=255)

                                                                                                                                       cccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAgg           >  1:1488587/1‑150 (MQ=255)

                                                                                                                                        ccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgc          <  1:1837220/150‑1 (MQ=255)

                                                                                                                                         cGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTg                                                                             <  1:1300352/82‑1 (MQ=255)

                                                                                                                                         cGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTg                                                                             >  2:1300351/1‑82 (MQ=255)

                                                                                                                                            aCCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGAAAGTGCTGATCCAGAGTTTCGTATTAAGGTCGTAAAAGCTGTGATCGAGCAGCACATGTCCCTAAATCAGgctgct      >  1:1779115/1‑150 (MQ=255)

                                                                                                                                              ccAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgctgc    <  2:1917448/150‑1 (MQ=255)

TTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTC  >  NC_000913/3720648‑3720941

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