Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,311,784	C→A	A87S (GCC→TCC)	alsB ←	D‑allose ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,311,784	0	C	A	100.0%	233.3 / NA	66	A87S (GCC→TCC)	alsB	D‑allose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (40/26);  total (40/26)

AAGCTTCCACATTGCCGCCAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATA  >  NC_000913/4311649‑4311920                                                                                                                                        |                                                                                                                                         aaGCTTCCACATTGCCGCCAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATa                                                                                                                                      >  1:338985/1‑140 (MQ=255)  aGCTTCCACATTGCCGCCAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGgaga                                                                                                                                        >  1:515944/1‑137 (MQ=255)       ccACATTGCCGCCAGCCTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTg                                                                                                                                >  2:580331/1‑140 (MQ=255)         ccaTTGCCGCCAGTTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTa                                                                                                                              <  1:1195846/139‑1 (MQ=255)         acaTTGCCGCCAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTa                                                                                                                              <  1:16364/140‑1 (MQ=255)          aaTTGCCGCCAGCTTTTTTCAGATTATCCATGTCGATTTTTTCACCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGAT‑CCCTTGTAAt                                                                                                                            <  1:79103/139‑1 (MQ=255)           aTTGCCGCCAGCCTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAAt                                                                                                                            <  1:580331/140‑1 (MQ=255)             tGCCGCCAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTa                                                                                                                              <  1:943011/136‑1 (MQ=255)             tGCCGCCAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTa                                                                                                                              >  2:943011/1‑136 (MQ=255)              gccgccAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAAtttt                                                                                                                         >  2:1321608/1‑140 (MQ=255)               ccgccAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAAttttt                                                                                                                        >  2:1293299/1‑140 (MQ=255)                cgccAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTtta                                                                                                                       <  1:1321608/140‑1 (MQ=255)                  ccAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTttatt                                                                                                                     >  1:990660/1‑140 (MQ=255)                   cAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTttatta                                                                                                                    <  1:1354270/140‑1 (MQ=255)                                  tATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAAt                                                                                                     <  1:1293299/140‑1 (MQ=255)                                   aTCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGCTACCTTTGTAATTTTTATTACTGAGATCTTCAAATa                                                                                                    >  1:1202747/1‑140 (MQ=255)                                          tCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCaa                                                                                             <  2:585110/140‑1 (MQ=255)                                          tCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCaa                                                                                             <  2:282559/140‑1 (MQ=255)                                              ttttttCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTga                                                                                         >  1:1038681/1‑140 (MQ=255)                                              ttttttCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTga                                                                                         >  1:1092463/1‑140 (MQ=255)                                                          gaTTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCg                                                                             <  2:645789/140‑1 (MQ=255)                                                                    aTAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGTTAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAg                                                                   >  2:95931/1‑140 (MQ=255)                                                                    aTAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGc                                                                            <  2:1124003/131‑1 (MQ=255)                                                                    aTAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGc                                                                            >  1:1124003/1‑131 (MQ=255)                                                                    aTAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAg                                                                   >  2:59841/1‑140 (MQ=255)                                                                    aTAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAg                                                                   >  2:1110071/1‑140 (MQ=255)                                                                    aTAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAg                                                                   >  1:716000/1‑140 (MQ=255)                                                                      aaaTGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGga                                                                 <  2:564889/140‑1 (MQ=255)                                                                      aaaTGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGga                                                                 <  2:990660/140‑1 (MQ=255)                                                                         tGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCg                                                                             <  1:28484/125‑1 (MQ=255)                                                                         tGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCg                                                                             >  2:28484/1‑125 (MQ=255)                                                                          gCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAgg                                                             <  2:1092463/140‑1 (MQ=255)                                                                              tttttCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGcc                                                                           >  1:342158/1‑122 (MQ=255)                                                                              tttttCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGcc                                                                           <  2:342158/122‑1 (MQ=255)                                                                              tttttCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCaaa                                                         >  2:8743/1‑140 (MQ=255)                                                                              tttttCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCaaa                                                         >  2:416656/1‑140 (MQ=255)                                                                              tttttCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCaaa                                                         >  1:1208004/1‑140 (MQ=255)                                                                                          cgGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCt                                             >  1:1263776/1‑140 (MQ=255)                                                                                                aCAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGcc                                       <  1:95931/140‑1 (MQ=255)                                                                                                         aCCAGATTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGtttt                              >  1:371318/1‑140 (MQ=255)                                                                                                              aTTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGccc                                      >  1:164536/1‑127 (MQ=255)                                                                                                              aTTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGccc                                      <  2:164536/127‑1 (MQ=255)                                                                                                              aTTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCtt                         >  2:450884/1‑140 (MQ=255)                                                                                                              aTTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCtt                         >  2:420760/1‑140 (MQ=255)                                                                                                              aTTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCtt                         >  2:380381/1‑140 (MQ=255)                                                                                                              aTTCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCtt                         >  1:1152395/1‑140 (MQ=255)                                                                                                               ttCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTc                        <  1:59841/140‑1 (MQ=255)                                                                                                                tCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCa                       <  2:1242638/140‑1 (MQ=255)                                                                                                                tCACTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCa                       >  2:1029383/1‑140 (MQ=255)                                                                                                                   cTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATGTCAACGCTGACGCCCAGTGTTTTTGCTTCATCt                    >  1:179797/1‑140 (MQ=255)                                                                                                                   cTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCt                          <  2:320127/134‑1 (MQ=255)                                                                                                                   cTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCt                          >  1:320127/1‑134 (MQ=255)                                                                                                                   cTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCa                       <  2:420649/137‑1 (MQ=255)                                                                                                                   cTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCa                       >  1:420649/1‑137 (MQ=255)                                                                                                                   cTGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCt                    >  1:740941/1‑140 (MQ=255)                                                                                                                    tGAGGATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCtt                   >  1:836231/1‑140 (MQ=255)                                                                                                                        gATAATGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAAt               >  1:287245/1‑140 (MQ=255)                                                                                                                            aTGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCt           <  2:740941/140‑1 (MQ=255)                                                                                                                            aTGGAGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCt           <  1:450884/140‑1 (MQ=255)                                                                                                                                aGCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCttttt       >  1:400067/1‑140 (MQ=255)                                                                                                                                 gCGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCtttttt      <  2:1208004/140‑1 (MQ=255)                                                                                                                                  cGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTc     >  2:527916/1‑140 (MQ=255)                                                                                                                                  cGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTc     >  1:1377719/1‑140 (MQ=255)                                                                                                                                  cGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTc     >  2:91632/1‑140 (MQ=255)                                                                                                                                  cGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTc     >  1:534564/1‑140 (MQ=255)                                                                                                                                  cGAAGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTc     >  1:1045406/1‑140 (MQ=255)                                                                                                                                    aaGGAGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTG‑TTTTGCTTCATCTTCAATGCCTTTTTTCata  >  2:731664/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                         AAGCTTCCACATTGCCGCCAGCTTTTTTCAGATTATCCATGTCGATTTTTTCATCGAGATTAACCAGATAAATGCCTTTTTTCCATGCGCGGGCGACAGGCATGACCAGATTCACTGAGGATAATGGAGCGAAGGCGATACCTTTGTAATTTTTATTACTGAGATCTTCAAATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATA  >  NC_000913/4311649‑4311920

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