Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,491 G→C 6.2% P14R (CCG→CGG)  iscS ← cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,661,4910GC 6.2% 409.3 / 14.4 113P14R (CCG→CGG) iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (53/53);  new base C (4/3);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 115 total aligned reads displayed.

ACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTC  >  NC_000913/2661347‑2661629
                                                                                                                                                 |                                                                                                                                          
aCCAGATCGGCAATCTGATTCCGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCg                                                                                                                                         <  3:941188/148‑1 (MQ=255)
     aTCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGtt                                                                                                                                    >  2:49292/1‑148 (MQ=255)
      tCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGt                                                                                                                                        <  4:571893/144‑1 (MQ=255)
        ggCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCg                                                                                                                                 <  2:823866/148‑1 (MQ=255)
         gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCgg                                                                                                                                           >  2:129940/1‑137 (MQ=255)
         gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCgg                                                                                                                                           <  1:129940/137‑1 (MQ=255)
               tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGagtagt                                                                                                                          <  1:49292/148‑1 (MQ=255)
                gATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTc                                                                                                                         >  2:721327/1‑148 (MQ=255)
                   tACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCCGGTGGACCCGCGTg                                                                                                                                       <  4:928504/131‑1 (MQ=255)
                   tACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCCGGTGGACCCGCGTg                                                                                                                                       >  3:928504/1‑131 (MQ=255)
                           atatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                      >  4:869428/1‑140 (MQ=255)
                           atatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                      <  3:869428/140‑1 (MQ=255)
                             atCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCCGGTGGACCCGCGt                                                                                                                                        >  1:1037843/1‑120 (MQ=255)
                                  ccGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGt                                                                                                                                        >  1:1108321/1‑115 (MQ=255)
                                  ccGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGt                                                                                                                                        <  2:1108321/115‑1 (MQ=255)
                                    gCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                                                                                                     >  2:327243/1‑148 (MQ=255)
                                      ttcttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                      >  2:23597/1‑129 (MQ=255)
                                      ttcttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                      <  1:23597/129‑1 (MQ=255)
                                      ttcttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTCTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAtt                                                                                                   >  1:1030165/1‑148 (MQ=255)
                                        cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGc                                                                                                 <  1:721327/148‑1 (MQ=255)
                                         ttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTgt                                                                                                                                      <  1:548915/111‑2 (MQ=255)
                                         ttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTgt                                                                                                                                      >  2:548915/1‑110 (MQ=255)
                                          tcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAt                                                                                                    >  2:1013267/1‑143 (MQ=255)
                                          tcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAt                                                                                                    <  1:1013267/143‑1 (MQ=255)
                                          tcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                               >  2:293845/1‑148 (MQ=255)
                                          tcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGACCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                               >  2:895823/1‑148 (MQ=255)
                                            aGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                >  3:540894/1‑113 (MQ=255)
                                             gCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                >  1:164465/1‑112 (MQ=255)
                                             gCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                <  2:164465/112‑1 (MQ=255)
                                             gCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTAGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGCAATTTCATTGCTCtat                                                                                            <  4:956094/148‑1 (MQ=255)
                                               cTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATaa                                                                                          <  1:293845/148‑1 (MQ=255)
                                                     gCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTg                                                                                                  >  3:1111808/1‑134 (MQ=255)
                                                     gCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTg                                                                                                  <  4:1111808/134‑1 (MQ=255)
                                                         aaaCGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATa                                                                                                                   <  1:398749/113‑1 (MQ=255)
                                                         aaaCGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATa                                                                                                                   >  2:398749/1‑113 (MQ=255)
                                                          aaCGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACAt                                                                               >  2:845997/1‑148 (MQ=255)
                                                                gagaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGc                                                                                                                                  <  4:662777/91‑1 (MQ=255)
                                                                   aaCGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGCACATCACTCAATg                                                                      <  1:845997/148‑1 (MQ=255)
                                                                    aCGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGc                                                                                                                                  >  3:298064/1‑87 (MQ=255)
                                                                    aCGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGc                                                                                                                                  <  4:298064/87‑1 (MQ=255)
                                                                    aCGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAt                                                                                                    >  3:619282/1‑117 (MQ=255)
                                                                    aCGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAt                                                                                                    <  4:619282/117‑1 (MQ=255)
                                                                     cGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                    <  4:887868/148‑1 (MQ=255)
                                                                     cGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                    >  3:887868/1‑148 (MQ=255)
                                                                       ggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCgg                                                                                                             <  1:620538/105‑1 (MQ=255)
                                                                       ggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCgg                                                                                                             >  2:620538/1‑105 (MQ=255)
                                                                          ggCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTcc                                                                                     <  1:774996/126‑1 (MQ=255)
                                                                           gCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAAt                                                              <  1:895823/148‑1 (MQ=255)
                                                                            ccGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATc                                                             <  4:232033/148‑1 (MQ=255)
                                                                                 ttACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                      <  1:1082809/86‑1 (MQ=255)
                                                                                 ttACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                      >  2:1082809/1‑86 (MQ=255)
                                                                                  tACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAt                                                                                                         <  2:156821/98‑1 (MQ=255)
                                                                                  tACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAt                                                                                                         >  1:156821/1‑98 (MQ=255)
                                                                                     cAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                      >  3:974488/1‑82 (MQ=255)
                                                                                     cAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                    <  3:525953/132‑1 (MQ=255)
                                                                                     cAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                    >  4:525953/1‑132 (MQ=255)
                                                                                      aaaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTgtt                                                                                                                                     >  4:902201/1‑65 (MQ=18)
                                                                                      aaaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                              <  3:460897/73‑1 (MQ=255)
                                                                                      aaaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                              >  4:460897/1‑73 (MQ=255)
                                                                                      aaaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                      >  1:208067/1‑81 (MQ=255)
                                                                                      aaaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCgg                                                   >  4:844654/1‑148 (MQ=255)
                                                                                        aGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTg                                                                                                                                   >  3:528975/1‑66 (MQ=255)
                                                                                        aGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTg                                                                                                                                   <  4:528975/66‑1 (MQ=255)
                                                                                        aGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCa                                                            <  4:817403/137‑1 (MQ=255)
                                                                                        aGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCa                                                            >  3:817403/1‑137 (MQ=255)
                                                                                        aGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGc                                                  >  2:592290/1‑147 (MQ=255)
                                                                                             ccGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTc                                                                                      >  3:440098/1‑106 (MQ=255)
                                                                                             ccGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTc                                                                                      <  4:440098/106‑1 (MQ=255)
                                                                                                tCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAgg                                                          <  4:180468/131‑1 (MQ=255)
                                                                                                tCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAgg                                                          >  3:180468/1‑131 (MQ=255)
                                                                                                 ccATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGCTAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                    <  2:333945/120‑1 (MQ=255)
                                                                                                 ccATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTc                                                                                                      >  4:172217/1‑86 (MQ=255)
                                                                                                 ccATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTc                                                                                                      <  3:172217/86‑1 (MQ=255)
                                                                                                 ccATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                    >  1:333945/1‑120 (MQ=255)
                                                                                                   aTCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCAc                                                                            <  4:555547/110‑1 (MQ=255)
                                                                                                   aTCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCAc                                                                            >  3:555547/1‑110 (MQ=255)
                                                                                                     cGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                                                                                                     <  2:314762/83‑1 (MQ=255)
                                                                                                     cGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                                                                                                     >  1:314762/1‑83 (MQ=255)
                                                                                                        cATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                                                                                                       <  4:525616/78‑1 (MQ=255)
                                                                                                        cATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgtgt                                 >  2:375072/1‑148 (MQ=255)
                                                                                                           aaaCTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGc                                                  <  1:349473/128‑1 (MQ=255)
                                                                                                                gCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGaa                                                               <  4:58501/110‑1 (MQ=255)
                                                                                                                gCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGaa                                                               >  3:58501/1‑110 (MQ=255)
                                                                                                                   tcatTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGttt                               >  1:931646/1‑139 (MQ=255)
                                                                                                                   tcatTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGttt                               <  2:931646/139‑1 (MQ=255)
                                                                                                                   tcatTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAt                      >  3:544663/1‑148 (MQ=255)
                                                                                                                    catTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAtt                     >  3:23663/1‑148 (MQ=255)
                                                                                                                     atTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                    >  3:866272/1‑100 (MQ=255)
                                                                                                                        ttCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGc                 <  4:23663/148‑1 (MQ=255)
                                                                                                                              gCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCa                                                                             >  3:289459/1‑82 (MQ=255)
                                                                                                                                aaCACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAg                                                           >  2:147908/1‑98 (MQ=255)
                                                                                                                                 acacGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGa                                                                >  1:316920/1‑92 (MQ=255)
                                                                                                                                 acacGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGa                                                                <  2:316920/92‑1 (MQ=255)
                                                                                                                                 acacGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                         <  3:297886/131‑1 (MQ=255)
                                                                                                                                 acacGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                         >  4:297886/1‑131 (MQ=255)
                                                                                                                                     cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAtt                     >  2:1112419/1‑132 (MQ=255)
                                                                                                                                     cgcgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAtt                     <  1:1112419/132‑1 (MQ=255)
                                                                                                                                       cgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCa                                                            <  4:889719/91‑1 (MQ=255)
                                                                                                                                       cgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCa                                                            >  3:889719/1‑91 (MQ=255)
                                                                                                                                        gggTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTc  <  1:375072/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
ACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTC  >  NC_000913/2661347‑2661629

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