BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,720,797	A→C	5.1%	R40R (AGG→CGG) ‡	insJ →	IS150 transposase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720,797	0	A	C	5.1%	482.4 / 12.0	157	R40R (AGG→CGG) ‡	insJ	IS150 transposase A
Reads supporting (aligned to +/- strand): ref base A (64/85); new base C (2/6); total (66/91)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

Only 100 of 159 total aligned reads displayed.

TTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTT  >  NC_000913/3720648‑3720946

 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)

       aCGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATg                                                                                                                                                <  1:1482018/150‑1 (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:675605/150‑5 (MQ=255)

                 ggAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCctc                                                                                                                                      >  1:2033957/1‑150 (MQ=255)

                         cccTAATATGTAAAAGCAAAAATACCCTTTTGAAAATCGCCTTGAAGTCTTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaaa                                                                                                                              <  2:1239228/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)

                               tatGTCAAAGCCAAAATACCCTTTTGAAAAGAGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGATACAGTATCATCTCGGCACGTTTTGGTGTCCCCCGACCCCCGGTCAGGACATGGGGTGCACTCTATGAAAAACAggg                                                                                                                        >  1:131040/1‑150 (MQ=255)

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

                                    cAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGaaa                                                                                                                   >  2:804589/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)

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

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

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

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

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

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

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

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

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

                                                                   tGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAg                                                                                    <  2:1072910/150‑1 (MQ=255)

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

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

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

                                                                       gTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCt                                                                                <  2:1128647/150‑1 (MQ=255)

                                                                         cGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGa                                                                              >  1:249059/1‑150 (MQ=255)

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

                                                                            gAATCACTACTTCACAACTGATGATGGTTACAGGATCCTCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTACCAAACCTAAAGGCGTTAGTGCTGATcc                                                                           >  1:82917/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                                 actactTCACAACTGATGATGGTTACAGCATCATCTCGGCACGTTTTGGTGACCCCCGAACCCAGGTCAGGACATGGGTTGCCCACTATGAAAAACAAGGAGAAAAACGTTTAATTCACAAACCTAAAGGCGTTAGTGCTGATCCAGAGt                                                                      >  1:1292107/1‑150 (MQ=255)

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

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

                                                                                        cacaACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGT‑CCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTAtt                                                              >  2:1972090/1‑150 (MQ=255)

                                                                                            aCTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAg                                                           <  1:43684/150‑1 (MQ=255)

                                                                                             cTGATGATGGTTACAGGATCATCTCGGCACGATTTGGTGTCCCCCGAACCCAGTTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAGGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAgg                                                          <  2:101810/150‑1 (MQ=255)

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

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

                                                                                               tatgatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTc                                                        <  1:926419/149‑1 (MQ=255)

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

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

                                                                                                       ttACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGc                                                <  2:578873/150‑1 (MQ=255)

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

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

                                                                                                          cAGGATCATCTCGGCATGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCtgt                                             >  2:425712/1‑150 (MQ=255)

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

                                                                                                               tcatcTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGTAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCg                                        <  1:1656414/150‑1 (MQ=255)

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

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

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

                                                                                                                      ggCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAgcagca                                 <  1:1961681/150‑1 (MQ=255)

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

                                                                                                                        cACGTTTTGGTGTCCCCCGATCCCAGGTCAGGTCATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAGCTATAGGCGTTAGTGCTGATCCAGAGTGGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGcaca                               <  1:548242/150‑1 (MQ=255)

                                                                                                                        cACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGcaca                               <  1:319257/150‑1 (MQ=255)

                                                                                                                        cACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGcaca                               <  2:1941397/150‑1 (MQ=255)

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

                                                                                                                          cGTTTTGGTGTCCACCGAACCCATGTCATGACATGGGTTGCCCTCTATGAAAAACATCGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCAGATCCAGAGTGGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATg                             <  1:346043/150‑1 (MQ=255)

                                                                                                                           gTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGGTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGt                            >  2:1038291/1‑150 (MQ=255)

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

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

                                                                                                                                ggTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCtt                       <  1:1606329/150‑1 (MQ=255)

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

                                                                                                                                    tCCCCCGAACCCAGTTCAGGACATGGGTTGCACTCTATGAAAAACATGGATAATAAGGTTTAATTCCCAAACCTATAGGCTTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATc                   <  2:131040/150‑1 (MQ=255)

                                                                                                                                    tCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATc                   >  1:1409078/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)

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

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

                                                                                                                                         cGAACCCAGGTCAGGTCATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgct              <  2:619472/150‑1 (MQ=255)

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

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

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

                                                                                                                                             cccAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgctgctg          >  2:2022369/1‑150 (MQ=255)

                                                                                                                                                aGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTc       <  1:61371/150‑1 (MQ=255)

                                                                                                                                                  ggcAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGATTTGCGTATTAAGGGCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCAc     <  2:190585/148‑1 (MQ=255)

                                                                                                                                                  ggcAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCAc     <  1:1609592/148‑1 (MQ=255)

                                                                                                                                                  gTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCAc     >  1:1783879/1‑150 (MQ=255)

TTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTT  >  NC_000913/3720648‑3720946

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: `‑`