Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,996 T→A 8.0% L246Q (CTG→CAG)  wecC → UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,9960TA 8.0% 477.5 / 25.4 138L246Q (CTG→CAG) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (62/65);  new base A (3/8);  total (65/73)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

Only 100 of 142 total aligned reads displayed.

GAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACAAA  >  NC_000913/3971848‑3972144
                                                                                                                                                    |                                                                                                                                                    
gAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCt                                                                                                                                                      >  1:391840/1‑149 (MQ=255)
   ggTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGc                                                                                                                                                   >  3:1443378/1‑149 (MQ=255)
     tGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGa                                                                                                                                                 >  4:1358523/1‑149 (MQ=255)
        gtgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATc                                                                                                                                              <  3:499984/149‑1 (MQ=255)
         tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATca                                                                                                                                             <  1:1225799/149‑2 (MQ=255)
         tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATca                                                                                                                                             <  2:1137904/149‑2 (MQ=255)
           tgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           <  3:808099/149‑1 (MQ=255)
           tgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATcagt                                                                                                                                           <  3:331272/149‑4 (MQ=255)
            gtcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCa                                                                                                                                          <  2:391840/149‑1 (MQ=255)
             tcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAc                                                                                                                                         >  4:1091328/1‑149 (MQ=255)
               gtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAccc                                                                                                                                       <  1:457526/149‑1 (MQ=255)
               gtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAccc                                                                                                                                       <  3:368837/149‑1 (MQ=255)
                  aCTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAAt                                                                                                                                               >  4:188452/1‑138 (MQ=255)
                  aCTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAAt                                                                                                                                               <  3:188452/138‑1 (MQ=255)
                      aCTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGtt                                                                                                                                <  1:940175/149‑1 (MQ=255)
                           acggataagccgcgcctgctgggggttctgtgccacgataaaccacggatcaacagcaatgcagtgaccgcccacgccagggccaggctgaagaatATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAatat                                                                                                                           >  4:815547‑M1/97‑149 (MQ=255)
                               cggcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTGATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttctt                                                                                                                       <  3:374826/149‑1 (MQ=255)
                                    gAAATGTGTAAGCTCATCGAAAACAGCTTCCGCGATATGAATATAGCTTTTGCTAATGAATTGCCGTTGATTTGTGCCGATCAGGGGATTAACGTCTGTGAACTGATTGGCCTGGCCAATCGTCACCCTCGCGTTAATAttctagatcg                                                                                                                  >  3:1233122/1‑143 (MQ=255)
                                     aaaTGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcct                                                                                                                 >  1:279253/1‑149 (MQ=255)
                                          tgtAAGCTCACCGAAAACAGTTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccc                                                                                                            >  3:59646/1‑149 (MQ=255)
                                           gtAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccct                                                                                                           <  3:1358523/149‑1 (MQ=255)
                                            tAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccctg                                                                                                          <  4:1443378/149‑1 (MQ=255)
                                               gCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttct                                                                                                                        <  3:979897/132‑1 (MQ=255)
                                               gCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttct                                                                                                                        >  4:979897/1‑132 (MQ=255)
                                                cTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGt                                                                                                      <  2:279253/149‑1 (MQ=255)
                                                  cACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTgg                                                                                                    >  2:1056852/1‑149 (MQ=255)
                                                     cGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGt                                                                                                                                 <  3:1094668/117‑1 (MQ=255)
                                                     cGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCg                                                                                                 <  2:931792/149‑1 (MQ=255)
                                                       aaaaCAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGt                                                                                               >  2:549660/1‑149 (MQ=255)
                                                        aaaCAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTgg                                                                                                    >  4:105260/1‑143 (MQ=255)
                                                        aaaCAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTgg                                                                                                    <  3:105260/143‑1 (MQ=255)
                                                            aGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTg                                                                                          <  3:995468/149‑1 (MQ=255)
                                                               ttCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCAt                                                                                       >  3:404255/1‑149 (MQ=255)
                                                                        gtgAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGa                                                                              <  1:1056852/149‑1 (MQ=255)
                                                                          gatcaacagcaatgcagtgaccgcccacgccagggccaggctgaagaatATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttct                                                                                                                        <  1:1124048‑M1/56‑1 (MQ=255)
                                                                          gAATATCTCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATc                                                                            >  2:786560/1‑149 (MQ=255)
                                                                              atCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCt                                                                                   <  2:481869/138‑1 (MQ=255)
                                                                              atCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCt                                                                                   >  1:481869/1‑138 (MQ=255)
                                                                              atCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTg                                                                        <  1:786560/149‑1 (MQ=255)
                                                                                cGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGt                                                                      <  1:549660/149‑1 (MQ=255)
                                                                                cGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGt                                                                      <  4:59646/149‑1 (MQ=255)
                                                                                 gCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           >  1:1246517/1‑79 (MQ=255)
                                                                                 gCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           <  2:1246517/79‑1 (MQ=255)
                                                                                   ttttGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTgg                                                                                                    >  1:61447/1‑116 (MQ=255)
                                                                                   ttttGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTa                                                                   <  1:963001/149‑1 (MQ=255)
                                                                                         tAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGGGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGTCGTGGGCGGTCAc                                                                                            >  4:1151232/1‑118 (MQ=255)
                                                                                         tAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGttt                                                                    <  4:37391/142‑1 (MQ=255)
                                                                                         tAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTgg                                                             <  1:1402757/149‑1 (MQ=255)
                                                                                         tAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTgg                                                             <  2:1198803/149‑1 (MQ=255)
                                                                                         tAATGAATTGTCGCTGATTTGGGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCAc                                                                                            <  3:1151232/118‑1 (MQ=255)
                                                                                          aaTGAATCGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGc                                                            >  3:889887/1‑149 (MQ=255)
                                                                                           aTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCt                                                                                                                                                      <  4:836347/58‑1 (MQ=255)
                                                                                           aTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccctg                                                                                                          >  1:177127/1‑102 (MQ=255)
                                                                                           aTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccctg                                                                                                          <  2:177127/102‑1 (MQ=255)
                                                                                           aTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTa                                                                   <  4:351006/141‑1 (MQ=255)
                                                                                           aTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTAt                                                                  <  1:1183395/142‑1 (MQ=255)
                                                                                                 tGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAc                                                     <  2:1176409/149‑1 (MQ=255)
                                                                                                   tCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACt                                                                                           <  4:686452/109‑1 (MQ=255)
                                                                                                   tCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACt                                                                                           >  3:686452/1‑109 (MQ=255)
                                                                                                   tCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGttt                                                                    >  4:9903/1‑132 (MQ=255)
                                                                                                   tCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGttt                                                                    <  3:9903/132‑1 (MQ=255)
                                                                                                     gCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccct                                                                                                           >  2:841053/1‑91 (MQ=255)
                                                                                                     gCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccct                                                                                                           <  1:841053/91‑1 (MQ=255)
                                                                                                     gCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAccccc                                                 >  1:1126488/1‑149 (MQ=255)
                                                                                                     gCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAccccc                                                 >  3:243133/1‑149 (MQ=255)
                                                                                                     gCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAc                                                     <  3:1447823/145‑1 (MQ=255)
                                                                                                      cTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCca                                                >  2:1195652/1‑149 (MQ=255)
                                                                                                         aTTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTggg                                                                                                   >  3:575679/1‑95 (MQ=255)
                                                                                                         aTTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTggg                                                                                                   <  4:575679/95‑1 (MQ=255)
                                                                                                             gtgCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttc                                                                                                                         <  2:410793/69‑1 (MQ=255)
                                                                                                             gtgCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttc                                                                                                                         >  1:410793/1‑69 (MQ=255)
                                                                                                              tgCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGgcgc                                        <  3:844292/149‑1 (MQ=255)
                                                                                                                 cGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGc                                     >  3:1182706/1‑149 (MQ=255)
                                                                                                                 cGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGc                                     <  4:889887/149‑1 (MQ=255)
                                                                                                                      aGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCtt                                   <  1:720326/146‑1 (MQ=255)
                                                                                                                      aGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCtt                                   >  2:720326/1‑146 (MQ=255)
                                                                                                                          gATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATcc                               <  2:1263311/146‑1 (MQ=255)
                                                                                                                             TAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTa                            <  1:676643‑M1/146‑1 (MQ=255)
                                                                                                                             TAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTa                            >  2:676643‑M1/1‑146 (MQ=255)
                                                                                                                                 GTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTg                                                                        >  2:805080‑M1/1‑98 (MQ=255)
                                                                                                                                 GTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgcgcg                     <  4:88170‑M1/149‑1 (MQ=255)
                                                                                                                                     gggAACTGATTCGCCTGGCGAATCGTCACCCTCGTGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGaa                 <  1:1195652/149‑1 (MQ=255)
                                                                                                                                     gggAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGaa                 <  4:1182706/149‑1 (MQ=255)
                                                                                                                                       gAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGt               >  2:1273192/1‑149 (MQ=255)
                                                                                                                                          tctaTTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCGTTCCCTGGCGTGGGAGGTCCCGGGATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTAGCGCGTGCGTAGTGaa            <  4:351889/146‑1 (MQ=255)
                                                                                                                                           tGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAAc           <  4:1347847/149‑1 (MQ=255)
                                                                                                                                           tGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAACCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAAc           <  3:353102/149‑1 (MQ=255)
                                                                                                                                            gATTCGCCAGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTg                                                              >  1:833467/1‑97 (MQ=255)
                                                                                                                                            gATTCGCCAGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATc                                >  3:1004346/1‑127 (MQ=255)
                                                                                                                                            gATTCGCCAGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAAc           <  2:79187/148‑1 (MQ=255)
                                                                                                                                            gATTCGCCAGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAAc           >  1:79187/1‑148 (MQ=255)
                                                                                                                                             aTTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGaa                 >  4:165054/1‑141 (MQ=255)
                                                                                                                                              ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTggg                                                                                                   <  1:1233081/58‑1 (MQ=255)
                                                                                                                                              ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTgg                                                             >  4:342506/1‑96 (MQ=255)
                                                                                                                                              ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTgg                                                             <  3:342506/96‑1 (MQ=255)
                                                                                                                                              ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAAcaga        >  3:351889/1‑146 (MQ=255)
                                                                                                                                                cGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATca      <  4:404255/149‑1 (MQ=255)
                                                                                                                                                   cTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACaa   >  3:1267840/1‑149 (MQ=255)
                                                                                                                                                   cTGGCGAATCGGCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACaa   >  2:367045/1‑149 (MQ=255)
                                                                                                                                                    tGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACaaa  >  1:1276935/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACAAA  >  NC_000913/3971848‑3972144

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