Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,322,700	A→G	5.4%	intergenic (‑195/+33)	secG ← / ← glmM	preprotein translocase membrane subunit/phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,322,700	0	A	G	5.4%	295.8 / 5.8	93	intergenic (‑195/+33)	secG/glmM	preprotein translocase membrane subunit/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (41/46);  new base G (2/3);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

AATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACA  >  NC_000913/3322552‑3322848                                                                                                                                                     |                                                                                                                                                     aaTAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACa                                                                                                                                                      >  1:679222/1‑149 (MQ=255)            aGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtt                                                                                                                                           <  2:93612/149‑1 (MQ=255)             gTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttt                                                                                                                                          <  3:648481/149‑1 (MQ=255)              tttttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCtttt                                                                                                                                         <  1:100014/149‑1 (MQ=255)              tttttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCtttt                                                                                                                                         <  3:125475/149‑1 (MQ=255)                 ttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCGATTGTCGCCttttttc                                                                                                                                      <  2:987185/149‑3 (MQ=255)                 ttAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAg                                                                                                                                      <  2:1072030/149‑1 (MQ=255)                    aTCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCa                                                                                                                                   <  2:1066919/149‑1 (MQ=255)                     tCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAg                                                                                                                                  <  3:899333/149‑1 (MQ=255)                     tCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAg                                                                                                                                  >  2:247266/1‑149 (MQ=255)                       tttCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGtt                                                                                                                                >  3:779818/1‑149 (MQ=255)                          cccACTGAAGCGTGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCACTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATc                                                                                                                             >  3:307692/1‑149 (MQ=255)                          cccACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGGTGCGTTCGCTGAGAAACACGGCGACAATCGCCGCCTTTTTAGCCAGTTATc                                                                                                                             >  3:106388/1‑149 (MQ=255)                            cACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTa                                                                                                                           >  4:352504/1‑149 (MQ=255)                             aCTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTaa                                                                                                                          >  1:33509/1‑149 (MQ=255)                             aCTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTaa                                                                                                                          >  1:156730/1‑149 (MQ=255)                                       gtgtGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAAc                                                                                                                >  3:269035/1‑149 (MQ=255)                                       gtgtGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAAc                                                                                                                >  2:1089688/1‑149 (MQ=255)                                       gtgtGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAAc                                                                                                                <  1:1040286/149‑1 (MQ=255)                                       gtgtGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAAc                                                                                                                <  1:887093/149‑1 (MQ=255)                                          tgAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                                                                             <  2:339357/149‑1 (MQ=255)                                               aCTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTa                                                                                                        <  2:167288/149‑1 (MQ=255)                                               aCTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTa                                                                                                        <  1:928165/149‑1 (MQ=255)                                                cTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTa                                                                                                                           >  4:965962/1‑129 (MQ=255)                                                cTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTa                                                                                                                           <  3:965962/129‑1 (MQ=255)                                                cTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                                                                       <  4:985845/149‑1 (MQ=255)                                                cTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTAc                                                                                                       <  4:1023850/149‑1 (MQ=255)                                                    ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                                                                   >  4:646983/1‑149 (MQ=255)                                                    ccAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCa                                                                                                   <  1:24778/149‑1 (MQ=255)                                                     cAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCAt                                                                                                  >  4:107390/1‑149 (MQ=255)                                                        aGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                                                                               <  2:755823/149‑1 (MQ=255)                                                        aGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCgg                                                                                               >  4:47919/1‑149 (MQ=255)                                                         gCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGc                                                                                              >  3:294454/1‑149 (MQ=255)                                                           gAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGa                                                                                            <  3:95245/149‑1 (MQ=255)                                                             aTGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATg                                                                                          <  4:1059732/149‑1 (MQ=255)                                                              tGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGc                                                                                         <  3:895698/149‑1 (MQ=255)                                                                aCCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCgg                                                                                       <  3:47919/149‑1 (MQ=255)                                                                 ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAATCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGGCAATCGCCGCCTTTTTAGCCAGtt                                                                                                                                <  1:428665/107‑1 (MQ=255)                                                                 ccTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCAGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGtt                                                                                                                                >  2:428665/1‑107 (MQ=255)                                                                   ttCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGCCGATGCGgtgt                                                                                    <  3:698820/149‑1 (MQ=255)                                                                   ttCGCAAGGGCAATTTTAACGCAATGTATCAACTGCCGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGgtgt                                                                                    <  2:156230/149‑1 (MQ=255)                                                                     cGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGc                                                                                  <  1:795049/149‑1 (MQ=255)                                                                          gggCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAAtt                                                                             >  4:764691/1‑149 (MQ=255)                                                                            gCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCa                                                                           >  2:519232/1‑149 (MQ=255)                                                                             cAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAg                                                                          >  2:202326/1‑149 (MQ=255)                                                                                ttttATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCa                                                                       >  3:13972/1‑149 (MQ=255)                                                                                   tagcGCATTGTATCAACTGCGGAAAGAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAAGGGCTTTTACTGCATCGTCGATGCGGTGTGCAAATTCAGTCACCt                                                                    >  2:967628/4‑149 (MQ=255)                                                                                         aTTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCtt                                                              <  4:767990/149‑1 (MQ=255)                                                                                            gTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGt                                                           >  4:69278/1‑149 (MQ=255)                                                                                              aTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCt                                                         <  2:791956/149‑1 (MQ=255)                                                                                              aTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCt                                                         >  3:879238/1‑149 (MQ=255)                                                                                              aTCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCGTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCt                                                         >  2:881705/1‑149 (MQ=255)                                                                                               tCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCtt                                                        >  1:955692/1‑149 (MQ=255)                                                                                                cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc                                                       <  2:518758/149‑1 (MQ=255)                                                                                                cAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTc                                                       <  4:13972/149‑1 (MQ=255)                                                                                                  aCTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGc                                                     <  2:636345/149‑1 (MQ=255)                                                                                                   cTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGcc                                                    <  1:627199/149‑1 (MQ=255)                                                                                                    tGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                                                                             >  4:869615/1‑91 (MQ=255)                                                                                                    tGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                                                                             <  3:869615/91‑1 (MQ=255)                                                                                                       ggAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTcc                                                <  4:303065/149‑1 (MQ=255)                                                                                                         aaaaaaaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGGCAGGTATCTAACGCATTAAACGGCTTTTACTGCATTGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTccac                                              >  2:305854/1‑149 (MQ=255)                                                                                                            aaaaCAGCAAAACCGGTGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCAt                                           >  2:1067059/1‑149 (MQ=255)                                                                                                              aaCAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCa                                         >  3:546219/1‑149 (MQ=255)                                                                                                               aCAGCAAAACCGATGTGTTTGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCAc                                        <  4:150552/149‑1 (MQ=255)                                                                                                                 aGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGACAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCAcgc                                      <  4:101554/149‑1 (MQ=255)                                                                                                                   cAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGa                                    >  4:156799/1‑149 (MQ=255)                                                                                                                    aaaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGaa                                   >  2:948131/1‑149 (MQ=255)                                                                                                                     aaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAAt                                  >  4:131437/1‑149 (MQ=255)                                                                                                                     aaaCCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAAt                                  <  1:202326/149‑1 (MQ=255)                                                                                                                        ccGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTaa                               <  1:519232/149‑1 (MQ=255)                                                                                                                        ccGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTaa                               <  3:391876/149‑1 (MQ=255)                                                                                                                         cGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAAc                              <  4:188604/149‑1 (MQ=255)                                                                                                                           aTGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACgg                            >  1:213326/1‑149 (MQ=255)                                                                                                                            tgtgTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGt                           <  4:546219/149‑1 (MQ=255)                                                                                                                              tgtTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTc                         <  4:269035/149‑1 (MQ=255)                                                                                                                                 tCGCTGAGAAAAAAGGCGACTATCGCCGCCCTTTTAGGCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCGCCTGCGCTTCCTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGt                      <  2:470479/149‑1 (MQ=255)                                                                                                                                 tCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGt                      <  3:352504/149‑1 (MQ=255)                                                                                                                                    cTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGcc                   <  2:696402/149‑1 (MQ=255)                                                                                                                                    cTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGcc                   <  3:116031/149‑1 (MQ=255)                                                                                                                                     tGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCg                  >  4:283083/1‑149 (MQ=255)                                                                                                                                      gagaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCgg                 >  3:193899/1‑149 (MQ=255)                                                                                                                                       agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                                                                             >  3:501622/1‑56 (MQ=255)                                                                                                                                       agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGc                                                                                                             <  4:501622/56‑1 (MQ=255)                                                                                                                                       agaAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGa                <  2:200071/149‑1 (MQ=255)                                                                                                                                        ctaaaaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGAt               >  1:764414/3‑149 (MQ=255)                                                                                                                                         aaaaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGAtt              >  3:825169/1‑149 (MQ=255)                                                                                                                                           aaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCTGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTa            >  2:345260/1‑149 (MQ=255)                                                                                                                                           aaaaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTa            >  2:713012/1‑149 (MQ=255)                                                                                                                                            aaaGGCGGCGATTGTCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAATCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTAc           >  2:566376/1‑149 (MQ=255)                                                                                                                                             aaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACg          <  3:69278/149‑1 (MQ=255)                                                                                                                                             aaGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACg          >  2:42173/1‑149 (MQ=255)                                                                                                                                              aGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgc         >  2:949245/1‑149 (MQ=255)                                                                                                                                               ggCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgca        <  1:323290/149‑1 (MQ=255)                                                                                                                                                  gACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACgcagca     >  2:58552/1‑149 (MQ=255)                                                                                                                                                     aaTCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAacc  >  4:778487/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAATGACCTTCGCAAGGGCAATTTTATCGCATTGTATCAACTGCGGAAAAAAACAGCAAAACCGATGTGTTCGCTGAGAAAAAAGGCGACAATCGCCGCCTTTTTAGCCAGTTATCTAACGCTTTAAACGGCTTTTACTGCATCGGCGATGCGGTGTGCAAATTCAGTCACCTGCGCTTCGTCTTCGCCTTCCACCATCACGCGAATTAACGGTTCGGTGCCGGATTTACGCAGCAACA  >  NC_000913/3322552‑3322848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑