Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,694,840	G→A	H133Y (CAT→TAT)	bcsA ←	cellulose synthase, catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,694,840	0	G	A	100.0%	307.1 / NA	96	H133Y (CAT→TAT)	bcsA	cellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (46/50);  total (46/50)

AACTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATGCTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACC  >  NC_000913/3694707‑3694963                                                                                                                                      |                                                                                                                              aaCTGCGCCAGCGGGTTAAACGGCTGTGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACtgtg                                                                                                                          >  1:1342280/1‑139 (MQ=255)  aCTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTgtgt                                                                                                                         >  2:1657634/1‑139 (MQ=255)    tGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGttt                                                                                                                       <  2:1588726/139‑1 (MQ=255)      cgcCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTa                                                                                                                     >  1:1637290/1‑139 (MQ=255)        ccAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATg                                                                                                                   >  2:991085/1‑139 (MQ=255)        ccAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATg                                                                                                                   >  2:876044/1‑139 (MQ=255)        ccAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATg                                                                                                                   >  1:770918/1‑139 (MQ=255)             gggTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                              <  2:436264/139‑1 (MQ=255)             gggTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                              <  1:86343/139‑1 (MQ=255)             gggTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                              <  1:437004/139‑1 (MQ=255)                       gCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAGGCG‑‑CTCAgg                                                                                                  <  1:1111304/139‑1 (MQ=255)                        cTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGc                                                                                                 >  1:894233/1‑139 (MQ=255)                        cTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGc                                                                                                 >  1:818201/1‑139 (MQ=255)                          gAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAg                                                                                               <  2:894233/139‑1 (MQ=255)                          gAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAg                                                                                               <  1:1657634/139‑1 (MQ=255)                           aGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGc                                                                                              <  1:1455826/139‑1 (MQ=255)                             tAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTc                                                                                            >  1:348185/1‑139 (MQ=255)                                cGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATTGCTAACGCG‑‑CTCAGGCAGCTCGtt                                                                                         >  2:842978/1‑139 (MQ=255)                                cGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGtt                                                                                         <  1:876044/139‑1 (MQ=255)                                   aGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCAt                                                                                      <  2:687839/139‑1 (MQ=255)                                     aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAgg                                                                                                  >  1:832609/1‑125 (MQ=255)                                     aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAgg                                                                                                  <  2:832609/125‑1 (MQ=255)                                     aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATcc                                                                                    >  2:854675/1‑139 (MQ=255)                                     aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATcc                                                                                    <  2:1419943/139‑1 (MQ=255)                                     aTTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATcc                                                                                    >  2:823653/1‑139 (MQ=255)                                         aCGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTg                                                                                <  1:235716/139‑1 (MQ=255)                                         aCGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTg                                                                                <  1:892693/139‑1 (MQ=255)                                          cGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAg                                                                                               >  1:1505266/1‑123 (MQ=255)                                          cGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAg                                                                                               <  2:1505266/123‑1 (MQ=255)                                          cGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGa                                                                               >  1:1186211/1‑139 (MQ=255)                                               gAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCAt                                                                                      >  2:141360/1‑127 (MQ=255)                                               gAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCAt                                                                                      <  1:141360/127‑1 (MQ=255)                                               gAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTa                                                                          >  1:1146952/1‑139 (MQ=255)                                                          aaGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGaa                                                               >  1:798939/1‑139 (MQ=255)                                                                aCGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcgcg                                                                                                          <  1:511005/92‑1 (MQ=255)                                                                aCGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcgcg                                                                                                          >  2:511005/1‑92 (MQ=255)                                                                aCGATGATACCGAGGATCAACCGCCGCGCCCCGGCACGCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCTATTGATGGTAACGTCCACGAATAtttt                                                         >  2:1271958/1‑139 (MQ=255)                                                                 cgttgatACCGAGGATCAACCGCCGCGTCTCGGCACTCAAATGTCCTAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGATCGTTCATCCATTGATGGTAACGTCAACGTATATTTTg                                                        <  2:227630/136‑1 (MQ=255)                                                                   atgatACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCa                                                      <  2:1637290/139‑1 (MQ=255)                                                                   atgatACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCa                                                      <  1:524289/139‑1 (MQ=255)                                                                       tACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAg                                                                                               >  2:551997/1‑94 (MQ=255)                                                                       tACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAg                                                                                               <  1:551997/94‑1 (MQ=255)                                                                       tACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCt                                                                                             >  2:473473/1‑96 (MQ=255)                                                                       tACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCt                                                                                             <  1:473473/96‑1 (MQ=255)                                                                       tACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGtt                                                                                         <  1:1598158/100‑1 (MQ=255)                                                                       tACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGtt                                                                                         >  2:1598158/1‑100 (MQ=255)                                                                       tACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGAcc                                                  >  1:1517016/1‑139 (MQ=255)                                                                        aCCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcgcg                                                                                                          <  2:1351585/84‑1 (MQ=255)                                                                        aCCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAAcgcg                                                                                                          >  1:1351585/1‑84 (MQ=255)                                                                        aCCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCt                                                 >  1:1191074/1‑139 (MQ=255)                                                                          cGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                              <  1:708318/78‑1 (MQ=255)                                                                          cGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                              >  2:708318/1‑78 (MQ=255)                                                                          cGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                               <  2:1146952/139‑1 (MQ=255)                                                                          cGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                               <  1:1491345/139‑1 (MQ=255)                                                                          cGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                               >  1:1328924/1‑139 (MQ=255)                                                                                  aaCCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGccc                                       >  1:30858/1‑139 (MQ=255)                                                                                           gcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGccc                                       <  2:1613265/130‑1 (MQ=255)                                                                                           gcCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGccc                                       >  1:1613265/1‑130 (MQ=255)                                                                                               cGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                              <  1:1040734/57‑1 (MQ=255)                                                                                               cGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTaa                                                                                                              >  2:1040734/1‑57 (MQ=255)                                                                                                 gCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGtttt                                                                                                                      <  2:140077/47‑1 (MQ=255)                                                                                                 gCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGtttt                                                                                                                      >  1:140077/1‑47 (MQ=255)                                                                                                    cTCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCg                     <  2:709999/139‑1 (MQ=255)                                                                                                     tCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAg                                                     <  1:1149832/106‑1 (MQ=255)                                                                                                     tCAAATGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAg                                                     >  2:1149832/1‑106 (MQ=255)                                                                                                       aaaTGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGttt                  >  2:1353998/1‑139 (MQ=255)                                                                                                       aaaTGACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CGCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGttt                  >  2:311731/1‑139 (MQ=255)                                                                                                       aaaTGACCGAGCTCTTTCTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGttt                  >  1:1114438/1‑139 (MQ=255)                                                                                                           gACCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCttt              <  2:30858/139‑1 (MQ=255)                                                                                                            aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCtt               >  2:679515/1‑139 (MQ=255)                                                                                                            aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAgg                                                                                                  >  2:1385547/1‑54 (MQ=255)                                                                                                            aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAgg                                                                                                  <  1:1385547/54‑1 (MQ=255)                                                                                                            aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTg             >  1:50308/1‑139 (MQ=255)                                                                                                            aCCGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTg             >  1:1149232/1‑139 (MQ=255)                                                                                                              cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGAc                                                   <  1:789045/99‑1 (MQ=255)                                                                                                              cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGAc                                                   >  2:789045/1‑99 (MQ=255)                                                                                                              cGAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTgcg           <  2:475258/139‑1 (MQ=255)                                                                                                               gAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCAc                                                                  <  2:608591/83‑1 (MQ=255)                                                                                                               gAGCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCAc                                                                  >  1:608591/1‑83 (MQ=255)                                                                                                                 gCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCAt                                                                                      >  1:614170/1‑61 (MQ=255)                                                                                                                 gCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCAt                                                                                      <  2:614170/61‑1 (MQ=255)                                                                                                                 gCTCTTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGa        >  1:372618/1‑139 (MQ=255)                                                                                                                   tctTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATg      <  1:991085/139‑1 (MQ=255)                                                                                                                   tctTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATg      <  2:983103/139‑1 (MQ=255)                                                                                                                   tctTTTTTCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATg      <  2:818201/139‑1 (MQ=255)                                                                                                                     ttttttCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCa    <  2:1191074/139‑1 (MQ=255)                                                                                                                       ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                               >  2:57646/1‑94 (MQ=255)                                                                                                                       ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGa                                               <  1:57646/94‑1 (MQ=255)                                                                                                                       ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCAcc  <  1:1353998/139‑1 (MQ=255)                                                                                                                       ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCAcc  <  2:372618/139‑1 (MQ=255)                                                                                                                       ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCAcc  <  2:317677/139‑1 (MQ=255)                                                                                                                       ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCAcc  <  2:163709/139‑1 (MQ=255)                                                                                                                       ttttCTCATCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCAcc  <  2:1342280/139‑1 (MQ=255)                                                                                                                        tttCTCATCCAGATACTGTGTTTTATGGCTAACGCGCTCTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCAc   <  1:679515/139‑1 (MQ=255)                                                                                                                               tCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCa    <  1:1502649/129‑1 (MQ=255)                                                                                                                               tCCAGATACTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCa    >  2:1502649/1‑129 (MQ=255)                                                                                                                                      |                                                                                                                              AACTGCGCCAGCGGGTTAAACGGCTGAGTAACGCAGATTAACGCCAGAATCAGCGAGAAGGTGACGATGATACCGAGGATCAACCGCCGCGCCCCGGCACTCAAATGACCGAGCTCTTTTTTCTCATCCAGATGCTGTGTTTTATGGCTAACGCG‑‑CTCAGGCAGCTCGTTCATCCATTGATGGTAACGTCCACGAATATTTTGCAGACCTGAAAATGCCCGCCTGCGCGGTTTCGGCGTTTCTTTGCGCGATGCACC  >  NC_000913/3694707‑3694963

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