Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,109,766	A→G	K311K (AAA→AAG)	sbp →	sulfate/thiosulfate ABC transporter periplasmic binding protein Sbp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,109,766	0	A	G	100.0%	235.4 / NA	66	K311K (AAA→AAG)	sbp	sulfate/thiosulfate ABC transporter periplasmic binding protein Sbp
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/32);  total (34/32)

ACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAAGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCT  >  NC_000913/4109630‑4109899                                                                                                                                         |                                                                                                                                      aCTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGgcg                                                                                                                                    >  1:399959/1‑140 (MQ=255)     gCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCaaa                                                                                                                                >  2:1138700/1‑140 (MQ=255)      ccAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCaaaa                                                                                                                               <  2:197346/140‑1 (MQ=255)          aaGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAagag                                                                                                                           >  1:1036461/1‑140 (MQ=255)            ggTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGgcgc                                                                                                                                   >  1:1037577/1‑130 (MQ=255)            ggTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGgcgc                                                                                                                                   <  2:1037577/130‑1 (MQ=255)                   aaaTTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCt                                                                                                                  <  2:81747/140‑1 (MQ=255)                      ttGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAAc                                                                                                               <  2:130554/140‑1 (MQ=255)                      ttGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAAc                                                                                                               >  1:454556/1‑140 (MQ=255)                           gcgAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAAcggcgg                                                                                                          <  1:1167399/140‑1 (MQ=255)                           gcgAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAAcggcgg                                                                                                          <  1:924951/140‑1 (MQ=255)                           gcgAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAAcggcgg                                                                                                          <  2:341007/140‑1 (MQ=255)                           gagAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAAcggcgg                                                                                                          <  1:165811/138‑1 (MQ=255)                                     actacCGTCCGCGCGACGCTGAGGTGCCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCCAAAAGAGCATTTTGCTAACGGCGGTACGTTCgat                                                                                                <  2:1036461/140‑1 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGt                                                                                                     <  2:422568/128‑1 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGt                                                                                                     >  1:315917/1‑128 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGt                                                                                                     <  2:315917/128‑1 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGt                                                                                                     >  1:422568/1‑128 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCgatcagatca                                                                                         >  2:819450/1‑140 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCgatcagatca                                                                                         >  2:557277/1‑140 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCgat                                                                                                <  1:364833/133‑1 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCgat                                                                                                <  2:1093799/133‑1 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCgat                                                                                                >  2:364833/1‑133 (MQ=255)                                            tCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCgat                                                                                                >  1:1093799/1‑133 (MQ=255)                                              cgcgcgACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGGTCGATCAGATCAGc                                                                                       >  2:1112578/1‑140 (MQ=255)                                                   gACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACg                                                                                  <  1:1138700/140‑1 (MQ=255)                                                        tGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGAt                                                                             <  1:562767/140‑1 (MQ=255)                                                             tGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTccc                                                                        <  2:562019/140‑1 (MQ=255)                                                              ggCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTcccc                                                                       <  2:6592/140‑1 (MQ=255)                                                               gCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTccc                                                                        >  1:834871/1‑138 (MQ=255)                                                               gCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTccc                                                                        <  2:834871/138‑1 (MQ=255)                                                               gCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCa                                                                      <  2:828075/140‑1 (MQ=255)                                                               gCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCa                                                                      >  1:21138/1‑140 (MQ=255)                                                                cGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAg                                                                     <  1:960939/140‑1 (MQ=255)                                                                   aaaaGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGAt                                                                  <  2:21138/140‑1 (MQ=255)                                                                    aaaGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATa                                                                 >  2:632231/1‑140 (MQ=255)                                                                          cGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATc                                                           <  1:557277/140‑1 (MQ=255)                                                                                   gTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCggt                                                  <  1:1123581/140‑1 (MQ=255)                                                                                             cTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTcgc                                        >  1:114706/1‑140 (MQ=255)                                                                                               gAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGa                                      <  2:114706/140‑1 (MQ=255)                                                                                                aaGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGAc                                     >  2:431362/1‑140 (MQ=255)                                                                                                  gTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCggtggt                                               >  1:1050252/1‑128 (MQ=255)                                                                                                  gTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCggtggt                                               <  2:1050252/128‑1 (MQ=255)                                                                                                      ttCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGGTTCTCGCGACCGGGtt                               >  1:457788/1‑140 (MQ=255)                                                                                                       tCACCATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGttt                              >  2:717785/1‑140 (MQ=255)                                                                                                          ccATTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGtttttt                           >  1:405722/1‑140 (MQ=255)                                                                                                            aTTGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTAt                         >  1:966805/1‑140 (MQ=255)                                                                                                             ttGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTTTCGCGACCGGGTTTTTTAtt                        >  2:250119/1‑140 (MQ=255)                                                                                                             ttGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGAcc                                    >  1:169755/1‑128 (MQ=255)                                                                                                             ttGATGAAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGAcc                                    <  2:169755/128‑1 (MQ=255)                                                                                                                atgaAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGAGCAACAAACGCTGATTTTCCCAGGGTAATTATTAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGt                     >  2:285409/1‑140 (MQ=255)                                                                                                                  gaAGAGTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCa                   <  1:250119/140‑1 (MQ=255)                                                                                                                     gagTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCCGGATAATTATCAAAACCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCATGa                >  2:1117872/1‑140 (MQ=255)                                                                                                                     gagTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGa                >  1:418870/1‑140 (MQ=255)                                                                                                                     gagTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAAAGCTGATTTTCCCAGGATAATTATCAAAACCGGGGGGTTTTCGCGACCGGGGTTTTTTTTTTTCACGa                >  2:283172/1‑140 (MQ=255)                                                                                                                       gTTCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGAtt              >  1:51422/1‑140 (MQ=255)                                                                                                                        ttCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGAttt             >  1:1107676/1‑140 (MQ=255)                                                                                                                        ttCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGAttt             >  1:432444/1‑140 (MQ=255)                                                                                                                        ttCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGAttt             >  2:147011/1‑140 (MQ=255)                                                                                                                        ttCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGAtt              <  1:632011/139‑1 (MQ=255)                                                                                                                        ttCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGAtt              >  2:632011/1‑139 (MQ=255)                                                                                                                        ttCGGCGGCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGGGGTTTCCCGCGACCGGGGTTTTTATTTGTCACGAttt             >  2:346620/1‑140 (MQ=255)                                                                                                                          cggcggCTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTg           <  1:819450/140‑1 (MQ=255)                                                                                                                            gcggcTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCg         <  2:480901/140‑1 (MQ=255)                                                                                                                             cggcTGGACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGt        <  2:966805/140‑1 (MQ=255)                                                                                                                                   gACGAAGGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCt  <  1:632231/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      ACTCGCCAGAAGGTCAGGAAATTGCCGCGAAAAACTACTACCGTCCGCGCGACGCTGAGGTGGCGAAAAAGTACGAAAATGCGTTTCCAAAGCTGAAGTTATTCACCATTGATGAAGAGTTCGGCGGCTGGACGAAAGCGCAAAAAGAGCATTTTGCTAACGGCGGTACGTTCGATCAGATCAGCAAACGCTGATTTCCCCAGGATAATTATCAAACCCGGTGGTTTCTCGCGACCGGGTTTTTTATTTGTCACGATTTTGCGTTACCCT  >  NC_000913/4109630‑4109899

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