BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,720,802	A→C	5.5%	T41T (ACA→ACC)	insJ →	IS150 transposase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720,802	0	A	C	5.5%	388.9 / 10.4	128	T41T (ACA→ACC)	insJ	IS150 transposase A
Reads supporting (aligned to +/- strand): ref base A (60/60); new base C (2/5); total (62/66)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

Only 100 of 133 total aligned reads displayed.

GGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTATGCT  >  NC_000913/3720653‑3720951

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

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

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

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

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

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

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

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

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

          aTGGAATACCCCCTAATATGTCAAAGCCAAAATCCCCTTTTGAAAAGCGCCTAGCAGTCGTGAATCACAACTTCACAACTCAAGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGccc                                                                                                                                             <  1:343068/150‑1 (MQ=255)

           tGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                                                                            <  2:1587644/150‑1 (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)

                   ccccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaa                                                                                                                                    >  2:1270194/1‑150 (MQ=255)

                   ccccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGGGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGGGGCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGaaa                                                                                                                                    <  2:441895/150‑1 (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)

                           atGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGga                                                                                                                            <  1:1451786/150‑1 (MQ=255)

                           atGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGga                                                                                                                            >  2:1329455/1‑150 (MQ=255)

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

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

                            tGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCGCTCTATGAAAAACATGgag                                                                                                                           <  2:1025832/150‑1 (MQ=255)

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

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

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

                                   gCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAgg                                                                                                                    <  2:1472813/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                                                                                                             >  2:622667/1‑150 (MQ=255)

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

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

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

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

                                                 ttGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGAGCATCTCGGCCCGTTTTGGTGTCCCCCGACCCCAGGTCAGGACATGGGTTGCCGTCTATGAAAACCATGGAGAAAAAGGTTTAATTCCCAAAc                                                                                                      >  1:1646074/1‑150 (MQ=255)

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

                                                  tGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAAcc                                                                                                     >  1:662706/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)

                                                                 aGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGc                                                                                      >  2:551547/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                                                                                >  2:1568685/1‑150 (MQ=255)

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

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

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

                                                                          tCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGCACGCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATACCCAAACCTAAAGGCGTTAGTGCTGATCCata                                                                             >  2:316464/1‑148 (MQ=255)

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

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

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

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

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

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

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

                                                                                         tgatgaTGGTTACAGGAGCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGt                                                              <  2:821353/150‑1 (MQ=255)

                                                                                           atgatgGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGGCg                                                            <  1:640062/150‑1 (MQ=255)

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

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

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

                                                                                                 gTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGGGAAAg                                                      <  1:1285077/150‑1 (MQ=255)

                                                                                                    aCAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCtg                                                   <  1:969584/150‑1 (MQ=255)

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

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

                                                                                                                cGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGagcagc                                       >  2:1215289/1‑150 (MQ=255)

                                                                                                                  gCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGcac                                     <  1:429713/150‑1 (MQ=255)

                                                                                                                      gTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGt                                 <  2:177450/150‑1 (MQ=255)

                                                                                                                       ttttGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTc                                >  2:883887/1‑150 (MQ=255)

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

                                                                                                                          tGGTGTCCGCCTAAACCAGGGCAGGACTTGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGATTAATTCCCAAACCTATAGGCGTTAGTGCTGATCCAGAGTTGCGTATTTAGGTCGTGAAAGCTGGGATCGAGCAGCACATGTCCCt                             <  1:826750/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                cccccGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCa                       >  1:541602/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                                             gTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCAc          <  1:1607580/150‑1 (MQ=255)

                                                                                                                                                 ggACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTa      <  2:662706/150‑1 (MQ=255)

                                                                                                                                                  gACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTAt     >  2:60457/1‑150 (MQ=255)

                                                                                                                                                  gACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTAt     >  1:935689/1‑150 (MQ=255)

                                                                                                                                                  gACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTAt     >  1:1505033/1‑150 (MQ=255)

                                                                                                                                                    cATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTAtgc   >  2:546923/1‑150 (MQ=255)

                                                                                                                                                     aTGGGTTGCCCTATATGAAAAACAGGGAGAAAAATGTTTAATTCCCAAACCTAAAGGCGTTAGTTCTGATCCAGAGTTGCTTATTAAGGGCGTGAATGCTGTTATCTAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTAtgct  <  2:1726474/150‑1 (MQ=255)

GGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTATGCT  >  NC_000913/3720653‑3720951

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