Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,425,640	A→C	6.6%	*88Y (TAA→TAC)  pseudogene (11/453 nt)	ynaK → ydaY →	Rac prophage; conserved protein pseudogene, Rac prophage;Phage or Prophage Related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,425,640	0	A	C	6.6%	256.0 / 9.9	91	*88Y (TAA→TAC)  pseudogene (11/453 nt)	ynaK ydaY	Rac prophage; conserved protein pseudogene, Rac prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (40/45);  new base C (2/4);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

AATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCG  >  NC_000913/1425492‑1425783                                                                                                                                                     |                                                                                                                                                aaTTCGGCTTGACTAATCCTGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGGCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATaa                                                                                                                                                 <  2:754372/149‑1 (MQ=255) aaTTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATaa                                                                                                                                                 <  2:885570/149‑1 (MQ=255)  aTTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGTGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAAc                                                                                                                                                >  2:1075010/1‑149 (MQ=255)     cGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATTGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGAt                                                                                                                                             <  1:533512/149‑1 (MQ=255)     cGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGAt                                                                                                                                             <  2:1253942/149‑1 (MQ=255)      ggCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGATCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATc                                                                                                                                            <  1:435213/149‑1 (MQ=255)       gCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTATGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCa                                                                                                                                           >  1:281513/1‑149 (MQ=255)        cTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAg                                                                                                                                          >  1:73136/1‑149 (MQ=255)           gACTAATCCGGTGCTGATTGACGATAAGGGCGAAATTATTGCGGGGCACGGTCGTGTTATGGCTTCTGTAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGTCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTAt                                                                                                                                       <  1:1202221/149‑1 (MQ=255)            aCTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATc                                                                                                                                      <  1:381521/149‑1 (MQ=255)             cTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCg                                                                                                                                     >  2:129514/1‑149 (MQ=255)             cTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCg                                                                                                                                     >  2:763326/1‑149 (MQ=255)              tAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGc                                                                                                                                    >  1:1384890/1‑149 (MQ=255)              tAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGc                                                                                                                                    >  1:1066581/1‑149 (MQ=255)              tAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGATGAAATGCTCAAAAAGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGCCGGATGAGCAGAAGCAGCGATAACGATCAGTATCGc                                                                                                                                    >  1:1061221/1‑149 (MQ=255)               aaTCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATACTGATCGTTATCGCt                                                                                                                                   <  1:1408801/149‑1 (MQ=255)                aTCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCta                                                                                                                                  >  1:1137707/1‑148 (MQ=255)                aTCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTc                                                                                                                                  >  2:1142429/1‑149 (MQ=255)                aTCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGAAGGATGAGCAGAAGCAGCGATAACGAACAGTATCGCTc                                                                                                                                  >  1:1324357/1‑149 (MQ=255)                 tCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATACTGATCGTTATCGCtgc                                                                                                                                 <  2:881529/149‑3 (MQ=255)                 tCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATACTGATCGTTATCGCtgc                                                                                                                                 <  1:26139/149‑3 (MQ=255)                    ggTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATACTGATCGTTATCGCtgcttc                                                                                                                              <  2:399945/149‑6 (MQ=255)                       gCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATTGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAAt                                                                                                                           >  2:971157/1‑149 (MQ=255)                       gCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCATAATGGATTCTGTTCCGGTCTTTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAAt                                                                                                                           <  1:364364/149‑1 (MQ=255)                         tGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGc                                                                                                                         >  2:34690/1‑149 (MQ=255)                         tGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGc                                                                                                                         <  1:1513551/149‑1 (MQ=255)                         tGATTGACGAAAAGCGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGc                                                                                                                         <  2:543495/149‑1 (MQ=255)                          gATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCa                                                                                                                        <  2:552343/149‑1 (MQ=255)                           aTTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCAt                                                                                                                       <  2:392929/149‑1 (MQ=255)                           aTTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCAt                                                                                                                       >  2:428985/1‑149 (MQ=255)                              gACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGTCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTATTGCATTaa                                                                                                                    <  2:1029398/149‑1 (MQ=255)                                           aaaTTATTGCGCGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGATCCGGGCATTGTTCTGTCTGGCCTGACGGCTGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACAt                                                                                                       <  1:234624/149‑1 (MQ=255)                                             attattGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTg                                                                                                     <  2:911953/149‑1 (MQ=255)                                              ttattGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTga                                                                                                    <  2:370038/149‑1 (MQ=255)                                               tattGCGGGGCACGGGCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGCTCCGGGCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTgag                                                                                                   <  2:214916/149‑1 (MQ=255)                                                    cGGGTCACGGTCGTGGTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTACGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAAt                                                                                              >  1:720009/1‑149 (MQ=255)                                                     gggTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAAAGGATTCTGTTCCGGTCATTGTTATGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGCTCAGTATAGCACCCGTAATGCATTAATCAGTCGCCACATTGAGAAAATg                                                                                             >  2:1217514/1‑149 (MQ=255)                                                       gTCACGGTCGTGTTATGGCGGCAGAAAAGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAACAGAAGCAGCGATAACGATCAGAATCGCTCCCGTAATGCATTAATCCGTAGCCACATTGATAAAATgct                                                                                           >  2:215325/1‑147 (MQ=255)                                                       gTCACGGACGTGGTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCACTGTTCTGTCTGGCCTTACGGATGAGCAGTAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGa                                                                                           <  1:296186/149‑1 (MQ=255)                                                        gcACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGAt                                                                                          <  2:241821/148‑1 (MQ=255)                                                          aCGGTCGTGTTATGGCGGCTGAAATGCTAAAAATGGATTCTGTTACGGTAATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGc                                                                                        >  1:1360871/1‑149 (MQ=255)                                                             gTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAg                                                                                     >  2:891932/1‑149 (MQ=255)                                                               cGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGtt                                                                                   <  1:800041/149‑1 (MQ=255)                                                                   ttATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCa                                                                               >  2:340509/1‑149 (MQ=255)                                                                   ttATGGCGGCTGAAATGCTCAAAATGGATTCAGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCa                                                                               >  1:855116/1‑149 (MQ=255)                                                                    tATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGGAc                                                                              <  1:34690/149‑1 (MQ=255)                                                                    tATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCAc                                                                              <  2:73136/149‑1 (MQ=255)                                                                        gcggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTTTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCg                                                                          <  1:842761/149‑1 (MQ=255)                                                                        gcggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCg                                                                          <  2:644693/149‑1 (MQ=255)                                                                          ggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGa                                                                        >  1:775020/1‑149 (MQ=255)                                                                          ggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGa                                                                        >  1:769753/1‑149 (MQ=255)                                                                           gcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGaa                                                                       <  1:1075010/149‑1 (MQ=255)                                                                               aaaTGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGaa                                                                   >  1:804004/1‑149 (MQ=255)                                                                               aaaTGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGaa                                                                   <  2:1154788/149‑1 (MQ=255)                                                                                aaTGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAg                                                                  <  2:1328628/149‑1 (MQ=255)                                                                                   gCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATTCAAGTTTGCACGTCGGAACGAAGGat                                                               >  2:469989/1‑148 (MQ=255)                                                                                     tCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCACGTAATGAATTAATCCGTCGCCACATTGAGAAAATCGATGCCAGTTTGCACGTCGGAACGAAGGAGtt                                                             >  1:566096/1‑149 (MQ=255)                                                                                        aaaTGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCTTTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGa                                                          >  1:912735/1‑149 (MQ=255)                                                                                        aaaTGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGACACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGa                                                          >  1:310127/1‑149 (MQ=255)                                                                                          aTGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGata                                                        <  2:1589812/149‑1 (MQ=255)                                                                                             gATTCTGTTCCGGTAATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTGGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATAttt                                                     <  1:215325/149‑1 (MQ=255)                                                                                               ttCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCa                                                   >  1:1574679/1‑149 (MQ=255)                                                                                               ttCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCa                                                   >  1:918330/1‑149 (MQ=255)                                                                                                 cTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCaaa                                                 <  2:12117/149‑1 (MQ=255)                                                                                                 cTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCaaa                                                 <  2:1393812/149‑1 (MQ=255)                                                                                                 cTGTTCCGGTCATTGTTATTTCTGGCCAGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGACCGAAGGAGTTTGATATTTCaaa                                                 <  1:779337/149‑1 (MQ=255)                                                                                                  tGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAg                                                <  1:971157/149‑1 (MQ=255)                                                                                                       cGGTCATTGTGCTGTCTGGCCTGACGGAGGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTc                                           <  1:249073/149‑1 (MQ=255)                                                                                                          tCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTgag                                                                                                   <  1:614801/90‑1 (MQ=255)                                                                                                          tCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTgag                                                                                                   >  2:614801/1‑90 (MQ=255)                                                                                                           cATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAg                                       >  1:1218740/1‑149 (MQ=255)                                                                                                                  ctgtctgGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGAtt                                <  1:1142429/149‑1 (MQ=255)                                                                                                                   tgtctgACCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGCACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTc                               >  1:1401782/1‑149 (MQ=255)                                                                                                                      ctgGCCTGACGGATGATCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAATGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGt                            <  1:1368665/149‑1 (MQ=255)                                                                                                                      ctgGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGt                            <  2:44144/149‑1 (MQ=255)                                                                                                                         gCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTtga                         >  2:897843/1‑149 (MQ=255)                                                                                                                         gCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATCGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTtga                         <  1:1271704/149‑1 (MQ=255)                                                                                                                             gACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTgatgat                     >  1:727071/1‑149 (MQ=255)                                                                                                                              aCGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGAtt                    >  1:63802/1‑149 (MQ=255)                                                                                                                              aCGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGAtt                    <  1:129514/149‑1 (MQ=255)                                                                                                                                  aTGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACt                >  2:39512/1‑149 (MQ=255)                                                                                                                                   ttatcAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTc               <  2:1066581/145‑1 (MQ=255)                                                                                                                                   tGAGCAGAAGCAGCGATAACGATCTGTATCGCTCGCGTTATGCATTAATCCGTCGCCACATTTAGAAAATGGATGACAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCATAGGTTTCCGAGTTGGTTTCTGTTGATGATTTACTc               <  2:1061221/149‑1 (MQ=255)                                                                                                                                   tGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGCTGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTc               <  2:809865/149‑1 (MQ=255)                                                                                                                                     aGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGATGTGGATTCTGTTGATGATTTACTCAt             <  2:281513/149‑1 (MQ=255)                                                                                                                                      gCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCAtt            >  2:524928/1‑149 (MQ=255)                                                                                                                                        aGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGa          <  1:1191354/149‑1 (MQ=255)                                                                                                                                         gAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGAt         >  2:12304/1‑149 (MQ=255)                                                                                                                                          cgggagcGATACTGATCGTTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATa        >  1:490721/5‑149 (MQ=255)                                                                                                                                            gcagcGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAAt      <  2:432669/149‑1 (MQ=255)                                                                                                                                             gagcGATACTGATCGTTATCGCTCCCGTACTGCATTAATCCGTCGCCAGATTGAGTAAATGGATGCCAGTTTGCAGGTCGGAACGATGGATTTTGATATTTCAAAGGTTTCCGAGTTGGATTCTGTTGATGATTTACTGATTGATAATg     >  1:248107/2‑149 (MQ=255)                                                                                                                                             caacGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAGCGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATg     <  1:129342/146‑1 (MQ=255)                                                                                                                                                cGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTGTCCGAGGTGGGTTCTGGGGATGTTTTACTCTTTGGTAAGGCCg  >  2:272231/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                AATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCG  >  NC_000913/1425492‑1425783

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