Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	571,444	G→A	*184* (TAG→TAA)	ybcL →	inactive polymorphonuclear leukocyte migration suppressor, DLP12 prophage, UPF0098 family secreted protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	571,444	0	G	A	100.0%	244.4 / NA	75	*184* (TAG→TAA)	ybcL	inactive polymorphonuclear leukocyte migration suppressor, DLP12 prophage, UPF0098 family secreted protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/39);  total (36/39)

TTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTAT  >  NC_000913/571314‑571578                                                                                                                                   |                                                                                                                                       ttAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGagag                                                                                                                                >  2:1703175/1‑139 (MQ=255)  tAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                               >  2:30430/1‑139 (MQ=255)         tGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTgg                                                                                                                        >  2:69230/1‑139 (MQ=255)          gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAagg                                                                                                                                      <  2:1369021/124‑1 (MQ=255)          gggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGAGACAGTAagg                                                                                                                                      >  1:1369021/1‑124 (MQ=255)           ggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                      <  1:30430/139‑1 (MQ=255)           ggCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                      <  1:34395/139‑1 (MQ=255)             ctctAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATACTGGCAag                                                                                                                    <  1:1607365/139‑1 (MQ=255)                 aaaaaCTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                               <  2:1470377/124‑1 (MQ=255)                 aaaaaCTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGt                                                                                                                               >  1:1470377/1‑124 (MQ=255)                     aCTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            <  2:595614/139‑1 (MQ=255)                          aaaGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTgaga                                                                                                                                 >  1:418791/1‑113 (MQ=255)                          aaaGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTgaga                                                                                                                                 <  2:418791/113‑1 (MQ=255)                                   tGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                          >  2:573899/1‑127 (MQ=255)                                   tGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                          <  1:573899/127‑1 (MQ=255)                                   tGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATc                                                                                              >  1:76948/1‑139 (MQ=255)                                     tAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                          <  1:680722/125‑1 (MQ=255)                                     tAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                          >  2:680722/1‑125 (MQ=255)                                     tAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCga                                                                                            >  1:1673749/1‑139 (MQ=255)                                            tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATa                                                                                     >  1:1003844/1‑139 (MQ=255)                                            tAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATa                                                                                     >  1:2726/1‑139 (MQ=255)                                              aCTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                   <  2:1673749/139‑1 (MQ=255)                                                tCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                                 <  1:1504568/139‑1 (MQ=255)                                                tCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                                 <  2:2726/139‑1 (MQ=255)                                                    gCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCa                                                                             <  2:76948/139‑1 (MQ=255)                                                        aGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                                 <  1:1755134/131‑1 (MQ=255)                                                        aGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                                 >  2:1755134/1‑131 (MQ=255)                                                          cGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCt                                                                                                                          <  2:1470960/88‑1 (MQ=255)                                                          cGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCt                                                                                                                          >  1:1470960/1‑88 (MQ=255)                                                             tAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatca                                                                    <  2:1003844/139‑1 (MQ=255)                                                              aGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatcat                                                                   <  1:1403175/139‑1 (MQ=255)                                                                  ggTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATgaga                                                               >  1:125103/1‑139 (MQ=255)                                                                         tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                   >  2:1707405/1‑112 (MQ=255)                                                                         tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACg                                                                                   <  1:1707405/112‑1 (MQ=255)                                                                         tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATa                                                                           >  1:669426/1‑120 (MQ=255)                                                                         tGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATa                                                                           <  2:669426/120‑1 (MQ=255)                                                                             tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGtt                                                    >  1:1348085/1‑139 (MQ=255)                                                                              aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                      <  2:1577095/72‑1 (MQ=255)                                                                              aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                      <  2:1214339/72‑1 (MQ=255)                                                                              aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                      >  1:1577095/1‑72 (MQ=255)                                                                              aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGt                                                                                                                >  1:1530562/1‑78 (MQ=255)                                                                              aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGt                                                                                                                <  2:1530562/78‑1 (MQ=255)                                                                              aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                                 >  1:910513/1‑109 (MQ=255)                                                                              aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGa                                                                                 <  2:910513/109‑1 (MQ=255)                                                                              aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGCGAGAGTATGCTGGCa                                                                                                                      >  1:1214339/1‑72 (MQ=255)                                                                                 gCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGc                                                                                                                       <  1:131182/68‑1 (MQ=255)                                                                                 gCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGc                                                                                                                       >  2:131182/1‑68 (MQ=255)                                                                                   taataaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCct                                              <  2:125103/139‑1 (MQ=255)                                                                                                aCCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATg                                                                                                                            <  2:992723/48‑1 (MQ=255)                                                                                                aCCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATg                                                                                                                            >  1:992723/1‑48 (MQ=255)                                                                                                aCCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGc                                                                                                                           >  1:1828028/1‑49 (MQ=255)                                                                                                aCCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGc                                                                                                                           <  2:1828028/49‑1 (MQ=255)                                                                                                aCCGCTGAGATAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCtt                                 >  2:350313/1‑139 (MQ=255)                                                                                                       agaTAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCg                          <  1:69230/139‑1 (MQ=255)                                                                                                         aTAACACCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTa                        <  1:1703175/139‑1 (MQ=255)                                                                                                             cacCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                      <  2:662267/41‑1 (MQ=255)                                                                                                             cacCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCa                                                                                                                      >  1:662267/1‑41 (MQ=255)                                                                                                              acCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGt                   <  2:1348085/139‑1 (MQ=255)                                                                                                              acCAGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGt                   <  1:350313/139‑1 (MQ=255)                                                                                                                 aGTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAg                >  1:487665/1‑139 (MQ=255)                                                                                                                  gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTtgtg                             <  2:1224669/125‑1 (MQ=255)                                                                                                                  gTTTATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAgg               >  1:111635/1‑139 (MQ=255)                                                                                                                  gTTTATGAGATAAAGGAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATTAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTtgtg                             >  1:1224669/1‑125 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            >  1:958580/1‑45 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGa                                                                                                            <  2:958580/45‑1 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                          <  1:1310492/47‑1 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                          >  2:1310492/1‑47 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAggg              >  1:615939/1‑139 (MQ=255)                                                                                                                   tttATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAggg              >  1:62775/1‑139 (MQ=255)                                                                                                                     tATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTa                                                                                                    <  1:1239062/51‑1 (MQ=255)                                                                                                                     tATGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTa                                                                                                    >  2:1239062/1‑51 (MQ=255)                                                                                                                       tGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                          >  1:298794/1‑43 (MQ=255)                                                                                                                       tGAGATAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                                                                                          <  2:298794/43‑1 (MQ=255)                                                                                                                        gagaTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGgaaaga         <  2:111635/139‑1 (MQ=255)                                                                                                                           aTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTa                        >  1:323656/1‑121 (MQ=255)                                                                                                                           aTAAAGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTa                        <  2:323656/121‑1 (MQ=255)                                                                                                                               aGTAAGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTAt  <  2:487665/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       TTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCATGAGATTCATAGGGTTAAGCCTCTTATGCCAGCTTTGTGTCGTATCCGTCAGGGAAAGAAAGTTAT  >  NC_000913/571314‑571578

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