Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,486 C→G 7.7% D16H (GAC→CAC)  iscS ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,661,4860CG 7.7% 305.0 / 13.3 91D16H (GAC→CAC) iscScysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base G (2/5);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAAT  >  NC_000913/2661338‑2661634
                                                                                                                                                     |                                                                                                                                                    
agattgcCGATCTGGTCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTcc                                                                                                                                                     >  1:989316/6‑150 (MQ=255)
 cAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGGGAGAACGGGAGGCCGGGTCACCAAAGGGTCCGTCCATCGTCCTAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCa                                                                                                                                                    <  2:677039/150‑1 (MQ=255)
     gcCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCgg                                                                                                                                                >  1:423900/1‑150 (MQ=255)
      ccGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGc                                                                                                                                               >  1:1590426/1‑150 (MQ=255)
         aCCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTg                                                                                                                                            <  2:924215/150‑1 (MQ=255)
           cAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGt                                                                                                                                          <  1:1190127/150‑1 (MQ=255)
             gATCTGCAATCTGATTACGCGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCCGGTGGACCCGCGTgttgc                                                                                                                                        <  1:731627/150‑5 (MQ=255)
               tCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCCGGTGGACCCGCGTgttgccg                                                                                                                                      <  1:104273/150‑7 (MQ=255)
               tCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGTgttgcc                                                                                                                                       >  2:204929/1‑145 (MQ=255)
                 ggCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGCAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGa                                                                                                                                    >  2:1552829/1‑150 (MQ=255)
                 ggCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGa                                                                                                                                    <  1:17871/150‑1 (MQ=255)
                 ggCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGTgttgccga                                                                                                                                     <  1:678368/150‑8 (MQ=38)
                  gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACTGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                                   <  1:760896/150‑1 (MQ=255)
                  gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                                   <  1:160087/150‑1 (MQ=255)
                   cAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGagt                                                                                                                                  <  1:773578/150‑1 (MQ=255)
                     aTCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGT‑GTTGCCGagaaa                                                                                                                                <  1:204929/150‑4 (MQ=37)
                      tCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGT‑GTTGCCGagaaaa                                                                                                                               <  1:730014/150‑5 (MQ=37)
                       cTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTc                                                                                                                              <  1:271923/150‑1 (MQ=255)
                        tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                             >  2:1508031/1‑150 (MQ=255)
                        tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                             <  2:195817/150‑1 (MQ=255)
                        tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                             >  2:212149/1‑150 (MQ=255)
                            tACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGAt                                                                                                                         >  1:1033729/1‑150 (MQ=255)
                             aCGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATa                                                                                                                        >  1:306774/1‑150 (MQ=255)
                                 gCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATc                                                                                                                    <  1:1430331/150‑1 (MQ=255)
                                   gATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCgg                                                                                                                  <  2:1271872/150‑1 (MQ=255)
                                     tatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTa                                                                                                                >  1:1418135/1‑150 (MQ=255)
                                      atCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACTGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTaa                                                                                                               <  1:1421252/150‑1 (MQ=255)
                                       tCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAt                                                                                                              >  1:1225428/1‑150 (MQ=255)
                                         tACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                                                                                                            >  1:226813/1‑150 (MQ=255)
                                          aCCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTc                                                                                                           <  2:306774/150‑1 (MQ=255)
                                          aCCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTc                                                                                                           >  2:1194213/1‑150 (MQ=255)
                                            cGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAt                                                                                                         <  2:1444600/150‑1 (MQ=255)
                                             gCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAtt                                                                                                        >  2:1006613/1‑150 (MQ=255)
                                               ttcttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGc                                                                                                      <  2:500935/150‑1 (MQ=255)
                                                 cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                                    <  1:1006613/150‑1 (MQ=255)
                                                 cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                                    <  2:1418135/150‑1 (MQ=255)
                                                 cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                                    <  1:298328/150‑1 (MQ=255)
                                                  ttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctct                                                                                                   <  2:1004501/150‑1 (MQ=255)
                                                   tcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCta                                                                                                  >  2:151244/1‑150 (MQ=255)
                                                      gCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATaa                                                                                               <  2:865838/150‑1 (MQ=255)
                                                         tGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                            <  2:556966/150‑1 (MQ=255)
                                                          gCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATaa                                                                                               <  1:1433426/146‑1 (MQ=255)
                                                          gCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATaa                                                                                               >  2:1433426/1‑146 (MQ=255)
                                                          gCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTc                                                                                           <  1:957305/150‑1 (MQ=255)
                                                          gCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTc                                                                                           >  2:1482125/1‑150 (MQ=255)
                                                            cAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTGATTTCATTGCTCTATAAACTCCg                                                                                         >  1:388041/1‑150 (MQ=255)
                                                            cAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCg                                                                                         >  2:968330/1‑150 (MQ=255)
                                                              gCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTa                                                                                       <  1:1601209/150‑1 (MQ=255)
                                                                 gAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACAt                                                                                    <  1:787839/150‑1 (MQ=255)
                                                                 gAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACAt                                                                                    >  2:787839/1‑150 (MQ=255)
                                                                    aCGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCAc                                                                                 >  1:1118553/1‑150 (MQ=255)
                                                                       gTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCa                                                                              >  2:184367/1‑150 (MQ=255)
                                                                         gagaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAAt                                                                            >  1:1080792/1‑150 (MQ=255)
                                                                          agaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATg                                                                           >  2:74031/1‑150 (MQ=255)
                                                                            aaCGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                         <  1:1552829/150‑1 (MQ=255)
                                                                              cgtgAGGCCGGGTTACCAAAGGTTCCTTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTTGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTATATAAACTCCGTACATCACTCAATGCAAg                                                                       <  1:1636164/147‑1 (MQ=255)
                                                                              cGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAg                                                                       <  2:719462/150‑1 (MQ=255)
                                                                                ggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGa                                                                     >  2:1222276/1‑150 (MQ=255)
                                                                                   ggCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATc                                                                  <  2:1406872/150‑1 (MQ=255)
                                                                                    gCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCa                                                                 <  2:585845/150‑1 (MQ=255)
                                                                                       gggTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGc                                                              <  1:1482125/150‑1 (MQ=255)
                                                                                       gggTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGc                                                              <  2:423900/150‑1 (MQ=255)
                                                                                           tACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTAcc                                                          >  1:1141904/1‑150 (MQ=255)
                                                                                            aCCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCg                                                         >  1:1393283/1‑150 (MQ=255)
                                                                                             ccAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCgg                                                        >  2:448872/1‑150 (MQ=255)
                                                                                                aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTg                                                     <  2:1225428/150‑1 (MQ=255)
                                                                                                  ggTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGa                                                   <  1:18169/150‑1 (MQ=255)
                                                                                                  ggTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGa                                                   <  1:1508031/150‑1 (MQ=255)
                                                                                                    ttCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCTGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATg                                                 <  1:1377931/150‑1 (MQ=255)
                                                                                                       cGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTAc                                              >  1:845655/1‑150 (MQ=255)
                                                                                                         tCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGa                                            >  2:1579761/1‑150 (MQ=255)
                                                                                                             tCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgt                                        >  2:1559374/1‑150 (MQ=255)
                                                                                                                tCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGtt                                     >  1:355870/1‑150 (MQ=255)
                                                                                                                  aTAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTa                                   <  2:989316/150‑1 (MQ=255)
                                                                                                                    aaaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACg                                 <  2:1219814/150‑1 (MQ=255)
                                                                                                                         gCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTa                            >  2:597454/1‑150 (MQ=255)
                                                                                                                           atcatTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAtt                          >  1:1237504/1‑150 (MQ=255)
                                                                                                                           atcatTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAtt                          <  1:968330/150‑1 (MQ=255)
                                                                                                                              atTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAg                       <  1:312763/150‑1 (MQ=255)
                                                                                                                                  tctcGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACt                   >  1:868292/1‑150 (MQ=255)
                                                                                                                                   ctcGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTc                  >  2:617687/1‑150 (MQ=255)
                                                                                                                                   ctcCGCAACCACGCCGGTGGACCCGCGT‑GTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTc                  >  1:542129/1‑150 (MQ=255)
                                                                                                                                    tcGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTcc                 >  2:339358/1‑150 (MQ=255)
                                                                                                                                      ggCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCgg               >  1:1460027/1‑150 (MQ=255)
                                                                                                                                       gCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGc              >  2:915660/1‑150 (MQ=255)
                                                                                                                                       gCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTA‑TTAGCACTCCGGcc             >  1:1174881/1‑150 (MQ=255)
                                                                                                                                        cAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGcc             >  2:137532/1‑150 (MQ=255)
                                                                                                                                              cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTACGGCCTGAtt        >  2:820581/1‑150 (MQ=255)
                                                                                                                                                   gTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCATGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGaa   <  2:1128675/150‑1 (MQ=255)
                                                                                                                                                   gTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGaa   >  1:647323/1‑150 (MQ=255)
                                                                                                                                                   gTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGaa   <  2:1080792/150‑1 (MQ=255)
                                                                                                                                                    tCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAAt  <  2:1393283/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAAT  >  NC_000913/2661338‑2661634

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