BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,720,793	G→T	6.6%	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	6.6%	369.6 / 12.4	130	Q38H (CAG→CAT)	insJ	IS150 transposase A
Reads supporting (aligned to +/- strand): ref base G (58/61); new base T (3/6); total (62/68)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 135 total aligned reads displayed.

CATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCA  >  NC_000913/3720644‑3720942

cATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAg                                                                                                                                                       >  1:821353/1‑150 (MQ=255)

cATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAg                                                                                                                                                       >  2:429713/1‑150 (MQ=255)

   tttGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTc                                                                                                                                                    <  2:38746/150‑1 (MQ=255)

   tttGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTc                                                                                                                                                    <  2:452246/150‑1 (MQ=255)

     tgttgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTcct                                                                                                                                                  >  1:903522/1‑148 (MQ=255)

       ttgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGa                                                                                                                                                <  1:895236/150‑1 (MQ=255)

        tgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGAc                                                                                                                                               >  2:1473811/1‑150 (MQ=255)

         ggACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACa                                                                                                                                              <  2:206163/150‑1 (MQ=255)

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

          gACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACAt                                                                                                                                             >  2:165238/1‑150 (MQ=255)

           aCGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATg                                                                                                                                            <  1:243169/150‑1 (MQ=255)

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

                gAAATGGAATATCCCCTAATATGTCAAAGCAAAAATACCCTTTTCAAAAGCGACTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttc                                                                                                                                       <  2:1584119/150‑2 (MQ=255)

                 aaaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCAAAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGc                                                                                                                                      >  1:1308863/1‑150 (MQ=255)

                  aaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcgg                                                                                                                                     <  1:1174336/150‑4 (MQ=255)

                  aaTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                                                                     >  2:720602/1‑150 (MQ=255)

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

                   aTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcggg                                                                                                                                    <  2:903522/150‑5 (MQ=255)

                    tGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCATGTCCTGACCTGGGttcgggg                                                                                                                                   <  2:403907/150‑6 (MQ=255)

                      gAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCctct                                                                                                                                 <  2:1328001/150‑1 (MQ=255)

                         aaGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATg                                                                                                                              <  1:1551325/149‑1 (MQ=255)

                          agacccTAAAATGTCAAAGCCAAAATACCCTTATGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTCATGATGGTTACAGCATCATCTAGGCACGTTTTGGTGTCCCCCGAACCCAGGACAGGACAAGTGTTGCCCTCTATGa                                                                                                                             >  2:826750/4‑150 (MQ=255)

                           gCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTAACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCAAGGTCAGGACATGGGTTGCCCTCTATGaa                                                                                                                            >  2:146470/1‑150 (MQ=255)

                            ccccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGGGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGGGGCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaa                                                                                                                           <  2:441895/150‑1 (MQ=255)

                             cccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaaa                                                                                                                          >  2:1374039/1‑150 (MQ=255)

                               cTAATATGTCAAAGCCAAAAAACCCTTTTGAAAAGCGCCTTGAAGCCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGTTTGCCCTCTATGAAAAAc                                                                                                                        <  1:295903/150‑1 (MQ=255)

                                   tatGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATgg                                                                                                                    <  1:1473811/150‑1 (MQ=255)

                                   tatGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATgg                                                                                                                    <  1:410665/150‑1 (MQ=255)

                                   tatGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGAGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGATTTGGTGTCCCCCGAACCCAGGTAAGGACATGGGTTGCCCTCTATGAAAAACATgg                                                                                                                    >  2:1313229/1‑150 (MQ=255)

                                   tatGACAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATgg                                                                                                                    <  1:356874/150‑1 (MQ=255)

                                     tGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACa                                                                                                                                              <  1:1099355/122‑1 (MQ=255)

                                     tGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACa                                                                                                                                              >  2:1099355/1‑122 (MQ=255)

                                     tGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGGCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGgag                                                                                                                  <  2:1289865/150‑1 (MQ=255)

                                      gTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAAGTCAGGACATGGGTTGCCCTCTATGAAAAACATGgaga                                                                                                                 >  2:307020/1‑150 (MQ=255)

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

                                          aaGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTACAAAACTGATCATGGTTTCAGGAACATCTCGGCACGTTTTGGTTACCGCCGAACCCAGGTCAGGACATGGTTTGCCCGCTATTAAAAACATGTAGaaaaa                                                                                                             >  2:1641939/1‑150 (MQ=255)

                                            gCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAgg                                                                                                           <  2:1472813/150‑1 (MQ=255)

                                            gCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAgg                                                                                                           <  1:1104560/150‑1 (MQ=255)

                                             ccAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGt                                                                                                          <  1:1348075/150‑1 (MQ=255)

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

                                             ccAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGATCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGt                                                                                                          <  1:1628271/150‑1 (MQ=255)

                                             ccAAAATAACCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTAACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGGCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGt                                                                                                          <  2:1568021/150‑1 (MQ=255)

                                                   tACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAAtt                                                                                                    <  1:1632025/150‑1 (MQ=255)

                                                   tACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAAtt                                                                                                    >  2:622667/1‑150 (MQ=255)

                                                    aCCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTc                                                                                                   <  1:722864/150‑1 (MQ=255)

                                                       cTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCa                                                                                                >  2:933420/1‑150 (MQ=255)

                                                          ttGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAc                                                                                             <  2:1375158/150‑1 (MQ=255)

                                                           tGAAAAGCGCCTTGAAGTCGTTAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAcc                                                                                            <  1:1674565/150‑1 (MQ=255)

                                                           tGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAcc                                                                                            <  1:149753/150‑1 (MQ=255)

                                                           tGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAcc                                                                                            >  1:662706/1‑150 (MQ=255)

                                                             aaaaGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTATTGGTGTCCCCCGAACCCCGGGCAGGACATGGGTTGCCCAGTATGAAAAACATGGAGAAAAAGGTTTAATTACAAAACCTa                                                                                          >  1:1285953/1‑150 (MQ=255)

                                                                  cgcCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGc                                                                                     >  2:205762/1‑150 (MQ=255)

                                                                     cTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGCCATGGGTTGCCCTCTATGAAAAACATGGAGAAAAACGTTTAATTCCCAAACCTAAAGGAGtt                                                                                  >  2:123116/1‑150 (MQ=255)

                                                                     cTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGtt                                                                                  <  2:58154/150‑1 (MQ=255)

                                                                        gAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAACTAAAGGCGTTAGt                                                                               >  2:1525367/1‑150 (MQ=255)

                                                                         aaGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTg                                                                              <  1:1656231/150‑1 (MQ=255)

                                                                          aGTCGTGAATCACTACTTCTCAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGc                                                                             >  1:6222/1‑150 (MQ=255)

                                                                          aGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGAGTTAGTGc                                                                             >  2:829546/1‑150 (MQ=255)

                                                                           gTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGAACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTACCAAACCTAAAGGCGTTAGTGCt                                                                            >  1:1306540/1‑150 (MQ=255)

                                                                           gTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATAATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTACCAAACCTAAAGGCGTTAGTGCt                                                                            >  2:591576/1‑150 (MQ=255)

                                                                             cGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTACCAAACCTAAAGGCGTTAGTGCTGa                                                                          >  1:438672/1‑150 (MQ=255)

                                                                              gTGAATCACTACTTCACAACTGATGATGGTTTCAGGATCATCTCGGCACGTTTTGGTGTCACCCGAACACAGGTCAGGACATGGGTTGCCATCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGAt                                                                         <  2:451039/150‑1 (MQ=255)

                                                                              gTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGAt                                                                         <  2:1671932/150‑1 (MQ=255)

                                                                                gAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATcc                                                                       >  1:995374/1‑150 (MQ=255)

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

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

                                                                                    cACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCagag                                                                   <  2:271424/150‑1 (MQ=255)

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

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

                                                                                         ctTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                                                              <  2:1020235/150‑1 (MQ=255)

                                                                                            cacaACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTAt                                                           >  2:243194/1‑150 (MQ=255)

                                                                                              aaaCTGATGATGGTTACAGTATTTTCTCGTCACGTTTTTGTGTGCCCAGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAACCATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTa                                                         <  1:260765/149‑1 (MQ=255)

                                                                                               aaCTGATGATGGTTACCGGATCATCTCGCCACGTTTTGGTGTCCCCAGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCATAGTTGCGTATTaa                                                        >  1:286790/1‑150 (MQ=255)

                                                                                                  tgatgaTGGTTAGAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGt                                                     <  1:970803/150‑1 (MQ=255)

                                                                                                     tgatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGt                                                  <  1:1374039/150‑1 (MQ=255)

                                                                                                      gatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTg                                                 <  1:165238/150‑1 (MQ=255)

                                                                                                      gatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTg                                                 >  2:1579749/1‑150 (MQ=255)

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

                                                                                                              cAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACCTGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCtgt                                         >  2:1243191/1‑150 (MQ=255)

                                                                                                                   tcatcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCg                                    >  2:561188/1‑150 (MQ=255)

                                                                                                                   tcatcTCGGCACGTTTTGGCGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTCTGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCg                                    <  2:1247144/150‑1 (MQ=255)

                                                                                                                        tcGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagcag                               <  1:243194/150‑1 (MQ=255)

                                                                                                                              cGTTTTGGTGTCCCCCGAACCCAAGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATg                         <  1:307020/150‑1 (MQ=255)

                                                                                                                                 tttGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTcc                      >  1:1464607/1‑150 (MQ=255)

                                                                                                                                   tGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAGGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCt                    <  1:449696/150‑1 (MQ=255)

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

                                                                                                                                    ggTGTCCCCCGAACCCAGGTCACGACATGGGATCCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCAAAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCtt                   >  2:1580054/1‑150 (MQ=255)

                                                                                                                                     gtgtCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTg                                                                              <  1:1563321/90‑1 (MQ=255)

                                                                                                                                     gtgtCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTg                                                                              >  2:1563321/1‑90 (MQ=255)

                                                                                                                                      tgtCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTaa                 >  2:660269/1‑150 (MQ=255)

                                                                                                                                         cccccGACCCCATGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTAAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCa              <  2:1046634/150‑1 (MQ=255)

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

                                                                                                                                             gcaaCCCATGTCCTGACCTGGGTTGCCCTCTATGAAAAACATGGAGAAAAATGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTACATTAATCAGgct          >  1:1722934/3‑150 (MQ=255)

                                                                                                                                             gcaaCCCATGTCCTGACCTGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgct          >  1:493448/3‑150 (MQ=255)

                                                                                                                                               aaCCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgc        >  1:1549337/1‑150 (MQ=255)

                                                                                                                                                aCCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgct       <  1:622667/150‑1 (MQ=255)

                                                                                                                                                aCCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgct       <  1:720602/150‑1 (MQ=255)

                                                                                                                                                   cAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgctgct    <  2:686710/150‑1 (MQ=255)

                                                                                                                                                   cAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgctgct    >  2:815506/1‑150 (MQ=255)

                                                                                                                                                     ggTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCa  >  2:206466/1‑150 (MQ=255)

CATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCA  >  NC_000913/3720644‑3720942

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