Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,870,532	A→G	T50A (ACA→GCA)	cdgI →	putative c‑di‑GMP binding protein CdgI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,870,532	0	A	G	100.0%	285.6 / NA	74	T50A (ACA→GCA)	cdgI	putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/37);  total (37/37)

ATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCT  >  NC_000913/1870403‑1870656                                                                                                                                  |                                                                                                                             aTTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAAtt                                                                                                                    >  2:376921/1‑140 (MQ=255)   ttCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTaa                                                                                                                  <  2:846122/140‑1 (MQ=255)      cATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAg                                                                                                               >  1:1501030/1‑140 (MQ=255)           ggTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCaa                                                                                                          >  2:42526/1‑140 (MQ=255)           ggTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCaa                                                                                                          >  1:1353433/1‑140 (MQ=255)             ttATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAg                                                                                                        <  1:1189715/140‑1 (MQ=255)              tATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGc                                                                                                       >  2:1231923/1‑140 (MQ=255)               aTTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCt                                                                                                      >  1:618198/1‑140 (MQ=255)               aTTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCt                                                                                                      >  2:1316384/1‑140 (MQ=255)                     tAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTa                                                                                                >  2:269153/1‑140 (MQ=255)                     tAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTa                                                                                                >  1:954650/1‑140 (MQ=255)                        atatTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTAttt                                                                                             <  1:1316384/140‑1 (MQ=255)                         tatTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTAtttt                                                                                            <  2:618198/140‑1 (MQ=255)                            tttttATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAg                                                                                         >  2:1293574/1‑140 (MQ=255)                             ttttATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGt                                                                                        >  1:972261/1‑140 (MQ=255)                              tttATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGtt                                                                                       <  1:992578/140‑1 (MQ=255)                                   cGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTg                                                                                  <  1:1202456/140‑1 (MQ=255)                                     ttAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCt                                                                                >  2:635862/1‑140 (MQ=255)                                       aaCGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGa                                                                              >  2:1209822/1‑140 (MQ=255)                                          gAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGt                                                                           <  1:1231923/140‑1 (MQ=255)                                                      tGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAattttatttt                                                               <  2:96194/140‑1 (MQ=255)                                                       ggTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTg                                                              <  2:1353433/140‑1 (MQ=255)                                                       ggTCAAAACTATATATCCCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTg                                                              <  1:42526/140‑1 (MQ=255)                                                        gTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGAc                                                                             <  1:1054312/124‑1 (MQ=255)                                                        gTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGAc                                                                             >  2:1054312/1‑124 (MQ=255)                                                        gTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAattttattt                                                                >  1:1138530/1‑137 (MQ=255)                                                        gTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAattttattt                                                                <  2:1138530/137‑1 (MQ=255)                                                         tCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGTTGACGTTAATTTTATTTTGCa                                                            >  1:975625/1‑140 (MQ=255)                                                          cAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCAt                                                           <  1:1293574/140‑1 (MQ=255)                                                            aaaCTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATta                                                         <  1:376921/140‑1 (MQ=255)                                                            aaaCTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATta                                                         <  1:387414/140‑1 (MQ=255)                                                                 atatatCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTAcc                                                    <  2:1452545/140‑1 (MQ=255)                                                                       ctctGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATtact                                                       >  1:672634/1‑131 (MQ=255)                                                                       ctctGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATtact                                                       <  2:672634/131‑1 (MQ=255)                                                                       ctctGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTtatat                                              >  1:897381/1‑140 (MQ=255)                                                                          tGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAGAAGATCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAAc                                           >  2:882776/1‑140 (MQ=255)                                                                          tGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAAc                                           >  1:235476/1‑140 (MQ=255)                                                                          tGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAAc                                           >  2:227291/1‑140 (MQ=255)                                                                          tGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAAc                                           >  2:895241/1‑140 (MQ=255)                                                                            ccATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATtact                                                       >  2:1099553/1‑126 (MQ=255)                                                                            ccATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATtact                                                       <  1:1099553/126‑1 (MQ=255)                                                                            ccATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGa                                         >  1:204458/1‑140 (MQ=255)                                                                            ccATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGa                                         >  1:777586/1‑140 (MQ=255)                                                                               tCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGaaaa                                      >  1:440164/1‑140 (MQ=255)                                                                                cAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGC‑TTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAaca                                    >  2:674871/1‑140 (MQ=255)                                                                                  aTAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGTTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAacac                                   <  2:975625/140‑1 (MQ=255)                                                                                  aTAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAacac                                   <  2:235476/140‑1 (MQ=255)                                                                                    aaaaTCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACaa                                 <  1:227291/140‑1 (MQ=255)                                                                                        tCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGt                             <  2:440164/140‑1 (MQ=255)                                                                                         ccAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAacac                                   >  2:1051525/1‑133 (MQ=255)                                                                                         ccAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAacac                                   <  1:1051525/133‑1 (MQ=255)                                                                                           aGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAacac                                   >  1:406056/1‑131 (MQ=255)                                                                                           aGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAacac                                   <  2:406056/131‑1 (MQ=255)                                                                                            gCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTaaaa                         >  1:1475222/1‑140 (MQ=255)                                                                                                tCACATTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCa                     <  1:197592/140‑1 (MQ=255)                                                                                                  acaTTGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATc                   >  2:304717/1‑140 (MQ=255)                                                                                                     ttGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCca                <  1:804087/140‑1 (MQ=255)                                                                                                     ttGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCca                <  1:1209822/140‑1 (MQ=255)                                                                                                     ttGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCca                <  2:535015/140‑1 (MQ=255)                                                                                                      tGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCcac               >  1:793405/1‑140 (MQ=255)                                                                                                      tGCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCcac               <  2:777586/140‑1 (MQ=255)                                                                                                       gCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCcaca              <  1:635862/140‑1 (MQ=255)                                                                                                       gCTCTCCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCcaca              <  2:204458/140‑1 (MQ=255)                                                                                                          ctcCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAg             >  2:724644/1‑138 (MQ=255)                                                                                                          ctcCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAg             <  1:724644/138‑1 (MQ=255)                                                                                                          ctcCTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGAt           >  1:858363/1‑140 (MQ=255)                                                                                                            ccTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTaaa                          >  2:609618/1‑123 (MQ=255)                                                                                                            ccTTTTTATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTaaa                          <  1:609618/123‑1 (MQ=255)                                                                                                                tttATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGa            <  1:1059281/133‑1 (MQ=255)                                                                                                                tttATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGa            >  2:1059281/1‑133 (MQ=255)                                                                                                                  tATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCcac               <  2:175402/128‑1 (MQ=255)                                                                                                                  tATGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCcac               >  1:175402/1‑128 (MQ=255)                                                                                                                   aTGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCt  <  2:793405/140‑1 (MQ=255)                                                                                                                   aTGGTTTCTATGGGTGCACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCt  <  2:858363/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                             ATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCT  >  NC_000913/1870403‑1870656

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