Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,334,390 C→A 11.1% A179E (GCA→GAA)  cysB → N‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,334,3900CA11.1% 177.8 / 4.3 67A179E (GCA→GAA) cysBN‑acetylserine‑responsive cysteine regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base A (2/6);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATTGCTGATGCCGTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCCGCGTATCGTTTTCACG  >  NC_000913/1334252‑1334523
                                                                                                                                          |                                                                                                                                     
aaaTTGCTGATGCCGTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGGTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCcaaaaa                                                                                                                             >  2:715528/1‑149 (MQ=255)
aaaTTGCTGATGCCGTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCaaaaaa                                                                                                                             <  1:827900/149‑1 (MQ=255)
   ttGCTGATGCCGTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCCCCCGCTGGCAGGCAAAAAAGcc                                                                                                                          >  1:1499802/1‑149 (MQ=255)
   ttGCTGATGCCGTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGcc                                                                                                                          <  1:79907/149‑1 (MQ=255)
   ttGCTGATGCCGTCGCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGcc                                                                                                                          <  1:530732/149‑1 (MQ=255)
       tGATGCCGTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAACAAGCCAATa                                                                                                                      >  2:974377/1‑149 (MQ=255)
         aTGCCGTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATGTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCCCCCGCTGGAAGGCCAAAAAGCCAATAcc                                                                                                                    >  2:1460032/1‑149 (MQ=255)
          tGCCGTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCa                                                                                                                   <  1:1266308/149‑1 (MQ=255)
              gTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAGCCATTACAATTgaa                                                                                                              >  2:444150/1‑149 (MQ=255)
               gctctAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTgaa                                                                                                              <  1:1573882/148‑1 (MQ=255)
                ctctAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCCCTGGAAGGCAAAAGAGCCATTACCATTgaag                                                                                                             >  2:10715/1‑149 (MQ=255)
                     aaGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTCCCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGGAGAACTg                                                                                                        >  1:53066/1‑149 (MQ=255)
                       cgCAAGGCTTAGTTCGCTATCCCCACACAAGCGCTGCAGCTGTATCAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGgc                                                                                                      <  2:99255/148‑1 (MQ=255)
                        gCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGGTGCCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGgcg                                                                                                     >  1:442522/1‑149 (MQ=255)
                          aaTGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCAAGGTTACCCGATGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCa                                                                                                   >  1:414215/1‑149 (MQ=255)
                                 atatCGCTATCGCCACAGAGGCGCTGCATCTGTATGAAGATGTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTTTAGTCACCCCGGATCACCCGAAAAAAAAAAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCg                                                                                            <  2:1276014/146‑1 (MQ=255)
                                  tttCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGt                                                                                           >  2:696086/1‑149 (MQ=255)
                                   ttCGCTATCGCCCCAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTAACATTGAAAAACTGGCCCAACATCCGtt                                                                                          >  2:853610/1‑149 (MQ=255)
                                    tCGCTATCGCCACCGGAGCGCTGCATCTGTATGAAGATTTAGTGATGTTAACGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTg                                                                                         <  1:418016/149‑1 (MQ=255)
                                     cGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTgg                                                                                        <  1:1808913/149‑1 (MQ=255)
                                      gCTATCGCCACCGAAGCGCTGCCTCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCTGATCACCCGCTGGCAGGCAAAAGAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGt                                                                                       <  2:1108242/149‑1 (MQ=255)
                                      gCTAGCGCCACAGAAGCGCTGCATCGGTAGGAAGAGTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACTCGAAAAAAAAAAAAAAAGCCATTGCCATTGAAGAACTGGCGCAATATCCGTTggc                                                                                       <  2:1328303/149‑2 (MQ=255)
                                           cGCCACAAAAGCGCTGCAGCTGTATGAAGATTTAGTGATGGTACCGTGCTACGACTGGAATCGGGCTATTGTAGTCACTCCGGATCATCCAAAGAAAAAAAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACat                                                                                  <  2:857268/149‑1 (MQ=255)
                                               acaGAAGGGGTCCTTCTGCCTGAAGGGTTGGTGAGGCTCCCGTGCTCCCGCTGGAATGGGGCTATTGTAGTCCCGCCTTATCGCCGGCAAAAAAAAAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATAc                                                                              <  2:1154663/149‑1 (MQ=255)
                                               acaGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATAc                                                                              >  2:1188195/1‑149 (MQ=255)
                                               acaGAAGCGCTGCATCTGTAAGAAGAATTAGCGATGTTACCGTGCTCCCACTGGACTCGGGCTATTGTAGTCACTCCGGATCACACAAAAAAAAAAAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATAc                                                                              <  2:1777290/149‑1 (MQ=255)
                                                cctaaGCGCTGCCTCGGTATGGAGATGAAGTGACGTTACCGTGCTCCCCCTGTAATCGGGCTAGTGTATTCCCTCCGGATCTCCAGAAAAAAAAAAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATAcc                                                                             <  2:1577041/146‑1 (MQ=255)
                                                 aGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCt                                                                            >  2:1402944/1‑149 (MQ=255)
                                                 aGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCt                                                                            >  1:734888/1‑149 (MQ=255)
                                                      ggcTGTATCTGGATGAAGTTTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGc                                                                       <  2:583147/148‑1 (MQ=255)
                                                       gcTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCt                                                                      >  2:1451141/1‑149 (MQ=255)
                                                           cATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTAc                                                                  <  1:229888/149‑1 (MQ=255)
                                                           cATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTAc                                                                  <  1:733642/149‑1 (MQ=255)
                                                             tCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGAATCTTGCCGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCg                                                                <  1:1406751/149‑1 (MQ=255)
                                                                   tGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGtt                                                          <  2:1694007/149‑1 (MQ=255)
                                                                     aaGATTTAGTGATGTTACCGTGCTACCACCGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCa                                                        <  1:833598/149‑1 (MQ=255)
                                                                      aGATTTAGTGATGTTACCCTGCTACCACTGGACTCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAg                                                       <  2:1136285/149‑1 (MQ=255)
                                                                           tAGTGGTGTTACCGGGCTACCACGGGAATGGGGCTATTGTAGTCTCTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTg                                                  <  2:233300/149‑1 (MQ=255)
                                                                            aGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTgg                                                 >  1:1190413/1‑149 (MQ=255)
                                                                               gATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATa                                              <  2:599658/149‑1 (MQ=255)
                                                                                cggTTACCGGGATACCACTGGGAACGGGTAAATGCAGTCACGCCAAATCGCCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATAc                                             <  1:1229351/147‑1 (MQ=255)
                                                                                 ggTTACCGTGCTGCCACTGGGTGCGGTCTATGGTAGTCACTCCGGAGCACCCGGGGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACt                                            <  2:653573/148‑1 (MQ=255)
                                                                                        gTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCGGAACGGGATGCTGCCTTTa                                     >  2:460409/1‑149 (MQ=255)
                                                                                        gTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTa                                     >  2:1221838/1‑149 (MQ=255)
                                                                                        gTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTCCCTTTa                                     >  2:460587/1‑149 (MQ=255)
                                                                                        gTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCATTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTa                                     >  1:1337851/1‑149 (MQ=255)
                                                                                           ccaCCACTGGAACCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGACCTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATc                                  <  2:943134/147‑1 (MQ=255)
                                                                                            tACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATcg                                 >  2:1812247/1‑149 (MQ=255)
                                                                                             aCCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCACTATCCGTTGGTGACATATACCTTCGGCTTTACGGGACGTTCAGAACTGGTTACTGTCTTTAATcgc                                >  2:269649/1‑149 (MQ=255)
                                                                                              ccACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATcgcg                               >  2:103277/1‑149 (MQ=255)
                                                                                               cACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGGCGTTCAGAACTGGATACTGCCTTTAATCgcgc                              >  1:1767112/1‑149 (MQ=255)
                                                                                               cACTGGAATCGGGCTATTGTAGTCACTCCGGAACACCCGCTGGCAGGCCAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTCCCGGACGTTCAGAACTGGATCCTGCCTTTAATCgcgc                              >  2:113475/1‑149 (MQ=255)
                                                                                                  tGGAATCGGGCTTTTGTAGTCACTCCCGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAgg                           <  1:492106/149‑1 (MQ=255)
                                                                                                       tccggCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTaa                      <  1:1734939/146‑1 (MQ=255)
                                                                                                       tCGGGCTATTGCAGTCACGTCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTaa                      <  1:1653450/149‑1 (MQ=255)
                                                                                                         gggCTATTGTAGTCCCTCCGGAGCACCCGCTGGCAGGAAAAAAAGCCAGTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGATCAGAACTGGATACTGCCTTTAATCCCGCAGGGTTAAcg                    <  2:1105290/149‑1 (MQ=255)
                                                                                                           gCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAAcgcc                  <  1:425277/149‑1 (MQ=255)
                                                                                                               ttGTAGTCACTCCAACACACCCCCTGGCAGGCGAAAAAGCCACTACCTTTGAAGAACTTGCGCAGTATCCGTGGGTGGCATATACCTTCGGCTTTCCCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCcccgg              >  2:1062225/1‑148 (MQ=255)
                                                                                                                    gTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTACCGCCGCGTATCGt         >  1:161637/1‑149 (MQ=255)
                                                                                                                     tCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCCGCGTATCGtt        <  1:950595/149‑1 (MQ=255)
                                                                                                                      cACTCCGTGTCGCCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCCGCGTATCGttt       <  2:1482163/149‑1 (MQ=255)
                                                                                                                       aCTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCCGCGTATCGtttt      >  1:1216074/1‑149 (MQ=255)
                                                                                                                       aCTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCCGCGTATCGtttt      >  2:114564/1‑149 (MQ=255)
                                                                                                                         tCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAGTCGCGCAGGGTGAACGCCGGGTATCGTTTTCa    >  1:518568/1‑149 (MQ=255)
                                                                                                                         tCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCCGCGTATCGTTTTCa    <  1:577217/149‑1 (MQ=255)
                                                                                                                           cGGATGAACCGCTAGCAAACAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCCGCGTATCGTTTTCACg  <  2:1404967/149‑1 (MQ=255)
                                                                                                                           cGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCCGCGTATCGTTTTCACg  <  2:734888/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                     
AAATTGCTGATGCCGTCTCTAAAGGCAATGCTGATTTCGCTATCGCCACAGAAGCGCTGCATCTGTATGAAGATTTAGTGATGTTACCGTGCTACCACTGGAATCGGGCTATTGTAGTCACTCCGGATCACCCGCTGGCAGGCAAAAAAGCCATTACCATTGAAGAACTGGCGCAATATCCGTTGGTGACATATACCTTCGGCTTTACCGGACGTTCAGAACTGGATACTGCCTTTAATCGCGCAGGGTTAACGCCGCGTATCGTTTTCACG  >  NC_000913/1334252‑1334523

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