Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,472 G→C 6.4% A20A (GCC→GCG)  iscS ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,661,4720GC 6.4% 317.2 / 8.1 94A20A (GCC→GCG) iscScysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base G (39/49);  new base C (3/3);  total (42/52)
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.89e-01

ACGATTTCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGG  >  NC_000913/2661323‑2661621
                                                                                                                                                     |                                                                                                                                                      
aCGATTTCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCg                                                                                                                                                        <  2:946002/150‑1 (MQ=255)
aCGATTTCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCg                                                                                                                                                        <  2:182292/150‑1 (MQ=255)
  gATTTCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGc                                                                                                                                                      <  1:272252/150‑1 (MQ=255)
  gATTTCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGc                                                                                                                                                      <  1:1431889/150‑1 (MQ=255)
  gATTTCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGc                                                                                                                                                      >  2:730014/1‑150 (MQ=255)
   aTTTCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCa                                                                                                                                                     <  1:796594/150‑1 (MQ=255)
    tttCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCaa                                                                                                                                                    >  1:1128675/1‑150 (MQ=255)
     ttCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAa‑c                                                                                                                                                  <  1:285546/150‑1 (MQ=255)
      tCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAa‑ca                                                                                                                                                 <  2:54787/150‑1 (MQ=255)
      tCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAa‑ca                                                                                                                                                 <  2:1551420/150‑1 (MQ=255)
          gcgGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CAcgcg                                                                                                                                             >  1:171521/1‑150 (MQ=255)
            ggATCAGCGCCGACCAGATCGTCAATCTGATTACGGGCTATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCggg                                                                                                                                           <  1:1354223/150‑1 (MQ=255)
             gATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGt                                                                                                                                          >  1:195817/1‑150 (MQ=255)
               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:212149/1‑150 (MQ=255)
                                       tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                >  2:1508031/1‑150 (MQ=255)
                                       tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                <  2:195817/150‑1 (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:298328/150‑1 (MQ=255)
                                                                cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                       <  2:1418135/150‑1 (MQ=255)
                                                                cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                       <  1:1006613/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                                                                       >  2:787839/1‑150 (MQ=255)
                                                                                gAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACAt                                                                       <  1:787839/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                      
ACGATTTCACGCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGG  >  NC_000913/2661323‑2661621

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: ‑