Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,661,091	G→A	L147L (CTC→CTT)	iscS ←	cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,661,091	0	G	A	98.9%	328.7 / ‑6.0	91	L147L (CTC→CTT)	iscS	cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (42/48);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCTTCGCGC  >  NC_000913/2660995‑2661187                                                                                                 |                                                                                                 aaCGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAgg                                                                                               >  1:1178889/1‑100 (MQ=255) aaCGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAgg                                                                                               <  2:1966353/100‑1 (MQ=255)     tGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATgg                                                                                           >  2:2289426/1‑100 (MQ=255)      gatagataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGgt                                                                                          >  1:2542161/1‑100 (MQ=255)      gatagataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGgt                                                                                          >  1:5540/1‑100 (MQ=255)      gatagataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGgt                                                                                          >  1:2166883/1‑100 (MQ=255)        tagataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGgtgt                                                                                        >  1:2547558/1‑100 (MQ=255)        tagataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGgtgt                                                                                        <  1:2657453/100‑1 (MQ=255)          gataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTgtcg                                                                                      <  2:663911/100‑1 (MQ=255)          gataATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCAGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTgtcg                                                                                      >  1:1808377/1‑100 (MQ=255)               tGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACg                                                                                 >  2:2795289/1‑100 (MQ=255)                gCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGATATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGc                                                                                >  1:1087286/1‑100 (MQ=255)                 ccACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCa                                                                               >  2:469653/1‑100 (MQ=255)                  cACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCAt                                                                              >  2:1942529/1‑100 (MQ=255)                  cACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCAt                                                                              >  2:2047216/1‑100 (MQ=255)                     gAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCgc                                                                           <  1:300629/100‑1 (MQ=255)                      aGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCgct                                                                          <  1:883853/100‑1 (MQ=255)                        cACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCgctgc                                                                        <  2:1928076/100‑1 (MQ=255)                        cACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCgctgc                                                                        >  2:531589/1‑100 (MQ=255)                         aCGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCgctgct                                                                       <  1:2102747/100‑1 (MQ=255)                         aCGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCgctgct                                                                       <  2:1178897/100‑1 (MQ=255)                         aCGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCgctgct                                                                       <  1:531584/100‑1 (MQ=255)                           ggCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTc                                                                     <  2:2074856/100‑1 (MQ=255)                           ggCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTc                                                                     >  2:2257018/1‑100 (MQ=255)                            gCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCa                                                                    >  1:2803093/1‑100 (MQ=255)                            gCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCa                                                                    <  1:2287057/100‑1 (MQ=255)                            gCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCa                                                                    <  1:206009/100‑1 (MQ=255)                             cacaTTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCGCGCATCGCTGCTTCaa                                                                   <  2:1552305/100‑1 (MQ=255)                              acaTTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACAAAAGGATGGTGTCGTCACGCATCGCTGCTTCaa                                                                   <  1:565637/99‑1 (MQ=255)                               caTTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGt                                                                 >  1:2178470/1‑100 (MQ=255)                                 tttCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTc                                                               <  2:799637/100‑1 (MQ=255)                                  ttCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCt                                                              >  1:1818597/1‑100 (MQ=255)                                   tCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCtt                                                             <  1:2120012/100‑1 (MQ=255)                                     gCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTc                                                           <  2:681629/100‑1 (MQ=255)                                      ccGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCa                                                          >  1:1964881/1‑100 (MQ=255)                                      ccGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCa                                                          >  2:1689839/1‑100 (MQ=255)                                          tAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTc                                                      <  2:1796934/100‑1 (MQ=255)                                            gCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGa                                                    >  1:961778/1‑100 (MQ=255)                                             ccGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGAt                                                   <  1:557306/100‑1 (MQ=255)                                             ccGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGAt                                                   <  2:1132798/100‑1 (MQ=255)                                                  atatCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGc                                              <  2:470142/100‑1 (MQ=255)                                                  atatCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGc                                              <  2:2166898/100‑1 (MQ=255)                                                   tatCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGcc                                             <  1:2069498/100‑1 (MQ=255)                                                    atCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCg                                            <  2:1746049/100‑1 (MQ=255)                                                     tCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGt                                           >  1:2324169/1‑100 (MQ=255)                                                      ccTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGtt                                          <  1:2746204/100‑1 (MQ=255)                                                      ccTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGtt                                          >  1:796678/1‑100 (MQ=255)                                                      ccTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGtt                                          >  1:2812500/1‑100 (MQ=255)                                                       cTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTa                                         >  1:1499206/1‑100 (MQ=255)                                                        tGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTAc                                        >  1:1548411/1‑100 (MQ=255)                                                         gCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACg                                       >  2:1617497/1‑100 (MQ=255)                                                          caccacGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGc                                      <  2:152498/100‑1 (MQ=255)                                                          caccacGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGc                                      <  1:530785/100‑1 (MQ=255)                                                            ccacGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTg                                    >  2:1183558/1‑100 (MQ=255)                                                               cGCCGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCgg                                 >  2:283153/1‑100 (MQ=255)                                                                 ccGATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTg                               <  2:2626333/100‑1 (MQ=255)                                                                   gATTTCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGcc                             >  2:1824786/1‑100 (MQ=255)                                                                     tttCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATGATGCCGTTACGCTGCGGTGCCag                           <  2:1974933/100‑1 (MQ=255)                                                                     tttCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCag                           >  1:1017301/1‑100 (MQ=255)                                                                      ttCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCagg                          <  2:1619139/100‑1 (MQ=255)                                                                       tCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggt                         <  2:1699096/100‑1 (MQ=255)                                                                       tCGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggt                         <  2:2324188/100‑1 (MQ=255)                                                                       tCGTTATTTACGTGCATGATGGACACAAGGATGATGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGAGGTGCCaggt                         >  1:1850919/1‑100 (MQ=255)                                                                        cGTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggta                        <  2:1046940/100‑1 (MQ=255)                                                                         gTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggtag                       <  1:2674660/100‑1 (MQ=255)                                                                         gTTATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggtag                       >  2:1460348/1‑100 (MQ=255)                                                                          ttATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggtagg                      >  1:2629852/1‑100 (MQ=255)                                                                           tATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggtaggt                     <  2:1087294/100‑1 (MQ=255)                                                                           tATTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCaggtaggt                     >  1:137299/1‑100 (MQ=255)                                                                            aTTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTg                    <  2:2547578/100‑1 (MQ=255)                                                                            aTTTACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTg                    <  2:134623/100‑1 (MQ=255)                                                                             tttACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGa                   >  2:1884482/1‑100 (MQ=255)                                                                              ttACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGAc                  <  2:2098000/100‑1 (MQ=255)                                                                              ttACGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGAc                  >  2:222240/1‑100 (MQ=255)                                                                                 cGTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTc               >  1:2628999/1‑100 (MQ=255)                                                                                  gTGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCa              <  1:203565/100‑1 (MQ=255)                                                                                   tGCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCaa             >  2:255559/1‑100 (MQ=255)                                                                                    gCATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCaaa            >  1:2653755/1‑100 (MQ=255)                                                                                     cATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCaaaa           <  2:1609877/100‑1 (MQ=255)                                                                                     cATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCaaaa           >  1:2393743/1‑100 (MQ=255)                                                                                     cATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCaaaa           <  2:796684/100‑1 (MQ=255)                                                                                     cATGATGGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCaaaa           <  1:1222192/100‑1 (MQ=255)                                                                                      atgatgGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAAc          >  1:2734016/1‑100 (MQ=255)                                                                                       tgatgGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAAcc         <  2:1945239/100‑1 (MQ=255)                                                                                         atgGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCtt       <  1:469648/100‑1 (MQ=255)                                                                                         atgGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCtt       >  1:1862721/1‑100 (MQ=255)                                                                                          tgGACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCTTc      >  2:360370/1‑100 (MQ=255)                                                                                           ggACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCTTcg     <  2:571886/100‑1 (MQ=255)                                                                                           ggACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCTTcg     <  1:1377594/100‑1 (MQ=255)                                                                                           ggACACAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCTTcg     <  2:1548423/100‑1 (MQ=255)                                                                                             acacAAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCTTcgcg   <  1:1933946/100‑1 (MQ=255)                                                                                              cacaAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCTTCgcgc  >  2:2036246/1‑100 (MQ=255)                                                                                                 |                                                                                                 AACGTGATAGATAATGCCACGAGCACGGCACATTTCGCCGATAGCCGCGATATCCTGCACCACGCCGATTTCGTTATTTACGTGCATGATGGACACGAGGATGGTGTCGTCACGCATCGCTGCTTCAAGTTCTTTCAGGTCGATAATGCCGTTACGCTGCGGTGCCAGGTAGGTGACTTCAAAACCTTCGCGC  >  NC_000913/2660995‑2661187

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