Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,763,451	C→T	10.3%	V187M (GTG→ATG) ‡	sufB ←	component of SufBCD Fe‑S cluster assembly scaffold

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,763,451	0	C	T	10.3%	205.7 / 13.3	68	V187M (GTG→ATG) ‡	sufB	component of SufBCD Fe‑S cluster assembly scaffold
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base T (3/4);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

AAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTT  >  NC_000913/1763305‑1763598                                                                                                                                                   |                                                                                                                                                    aaaCTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACgg                                                                                                                                                   >  4:565230/1‑149 (MQ=255)      gCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCg                                                                                                                                              >  3:144390/1‑149 (MQ=255)       cccGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCGTGGTGCCgg                                                                                                                                             <  3:1107917/149‑2 (MQ=255)          ggTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGat                                                                                                                                          <  4:144390/149‑1 (MQ=255)           gTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGata                                                                                                                                         >  2:1050232/1‑149 (MQ=255)                ttCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTAc                                                                                                                                    >  1:157158/1‑149 (MQ=255)                                    aGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCga                                                                                                                                             >  2:1380847/1‑120 (MQ=255)                                    aGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCga                                                                                                                                             <  1:1380847/120‑1 (MQ=255)                                      gTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTg                                                                                                              <  1:403838/149‑1 (MQ=255)                                             gTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGcc                                                                                                       <  4:845793/149‑1 (MQ=255)                                               tCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCGTGGTGCCgg                                                                                                                                             <  1:299921/109‑2 (MQ=255)                                               tCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCGTGGTGCCgg                                                                                                                                             >  2:299921/1‑108 (MQ=255)                                                  aTCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGc                                                                                                        <  2:1074568/143‑1 (MQ=255)                                                  aTCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGc                                                                                                        >  1:1074568/1‑143 (MQ=255)                                                     gggCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCaaa                                                                                               <  1:67830/149‑1 (MQ=255)                                                      ggCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                                                                              <  2:157158/149‑1 (MQ=255)                                                         aGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGaa                                                                                           <  3:565230/149‑1 (MQ=255)                                                          gcgcACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAAc                                                                                          <  1:1050232/149‑1 (MQ=255)                                                              aCGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTa                                                                                                                                     <  3:93351/102‑1 (MQ=255)                                                              aCGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTa                                                                                                                                     >  4:93351/1‑102 (MQ=255)                                                              aCGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATc                                                                                                          >  1:934058/1‑129 (MQ=255)                                                              aCGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATc                                                                                                          <  2:934058/129‑1 (MQ=255)                                                              aCGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGaa                                                                                      >  1:726795/1‑149 (MQ=255)                                                                               aataaaCGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCg                                                                     >  1:511714/1‑149 (MQ=255)                                                                                       tACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGtttt                                                             >  4:1239321/1‑149 (MQ=255)                                                                                                agGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATaa                                                    >  1:443022/1‑149 (MQ=255)                                                                                                agGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATaa                                                    <  2:443022/149‑1 (MQ=255)                                                                                                   cTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAagta                                                 <  3:1239321/149‑1 (MQ=255)                                                                                                      ccgccgCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGc                                                                                                                                  >  4:346833/1‑65 (MQ=255)                                                                                                      ccgccgCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGc                                                                                                                                  <  3:346833/65‑1 (MQ=255)                                                                                                      ccgccgCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCg                                                                                                                        <  4:91866/75‑1 (MQ=255)                                                                                                      ccgccgCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCg                                                                                                                        >  3:91866/1‑75 (MQ=255)                                                                                                           gcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTAc                                                                                                                                    >  4:576178/1‑58 (MQ=255)                                                                                                           gcATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTAc                                                                                                                                    <  3:576178/58‑1 (MQ=255)                                                                                                            cATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTcc                                                                                                                         >  1:1245405/1‑68 (MQ=255)                                                                                                            cATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTcc                                                                                                                         <  2:1245405/68‑1 (MQ=255)                                                                                                            cATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTc                                                           >  4:1350881/1‑130 (MQ=255)                                                                                                            cATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTc                                                           <  3:1350881/130‑1 (MQ=255)                                                                                                              ttAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCga                                                                                                                                             >  4:1026192/1‑46 (MQ=255)                                                                                                              ttAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCga                                                                                                                                             <  3:1026192/46‑1 (MQ=255)                                                                                                                  gcgcGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTg                                              <  1:1451472/137‑1 (MQ=255)                                                                                                                  gcgcGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTg                                              >  2:1451472/1‑137 (MQ=255)                                                                                                                         aaaGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATaa                                                    <  1:591501/124‑1 (MQ=255)                                                                                                                         aaaGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATaa                                                    >  2:591501/1‑124 (MQ=255)                                                                                                                          aagaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTcc                                                                             >  4:478954/1‑98 (MQ=255)                                                                                                                          aagaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTcc                                                                             <  3:478954/98‑1 (MQ=255)                                                                                                                            gaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCg                                                                     <  2:479787/104‑1 (MQ=255)                                                                                                                            gaagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCg                                                                     >  1:479787/1‑104 (MQ=255)                                                                                                                             aagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                                                                              >  3:1341541/1‑78 (MQ=255)                                                                                                                             aagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                                                                              >  1:373582/1‑78 (MQ=255)                                                                                                                             aagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                                                                              <  2:373582/78‑1 (MQ=255)                                                                                                                             aagTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAg                                                                                              <  4:1341541/78‑1 (MQ=255)                                                                                                                              agTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATg                      <  2:511714/149‑1 (MQ=255)                                                                                                                                ttGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGa                                <  3:837287/137‑1 (MQ=255)                                                                                                                                ttGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGa                                >  4:837287/1‑137 (MQ=255)                                                                                                                                 tGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCa                   <  2:726795/149‑1 (MQ=255)                                                                                                                                   tCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGa                                      >  2:1272310/1‑128 (MQ=255)                                                                                                                                   tCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGa                                      <  1:1272310/128‑1 (MQ=255)                                                                                                                                     aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAagtagt                                               <  2:452788/117‑1 (MQ=255)                                                                                                                                     aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAagtagt                                               >  1:452788/1‑117 (MQ=255)                                                                                                                                     aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGa                                <  2:980958/132‑1 (MQ=255)                                                                                                                                     aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGa                                >  1:980958/1‑132 (MQ=255)                                                                                                                                     aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCa               >  1:551940/1‑149 (MQ=255)                                                                                                                                     aTTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCa               >  3:790788/1‑149 (MQ=255)                                                                                                                                        tctcGGCACCGTGGTGCCGGGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTc                                                           <  1:384654/99‑1 (MQ=255)                                                                                                                                        tctcGGCACCGTGGTGCCGGGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTc                                                           >  2:384654/4‑102 (MQ=255)                                                                                                                                          tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGc                                                                         >  4:180210/2‑86 (MQ=255)                                                                                                                                          tcGGCACCGTGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGc                                                                         <  3:180210/85‑1 (MQ=255)                                                                                                                                                  acGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCttt  <  2:551940/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AAACTGCCCGGTTTTTTCTGCGTTAATGCGAAAATAGGTGGAAAGTTCCATCGGGCAGCGCACGCCTTTAGGCACATAAATAAACGTACCATCAGAGGCTACCGCCGCATTAAGCGCGGCAAAGAAGTTGTCATTCCCCGGCACCACGGTGCCGAGATATTTACGCACCAGTTCCGGGTGATCGTGGATCGCCTCACCAAAGGAACAGAAAATAATTCCCTGCTCCGCCAGTTTTTCGCGATAAGTAGTGGCAACCGAAACTGAGTCGAAAATGGCATCCACCGCCACCTCTTT  >  NC_000913/1763305‑1763598

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