Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,996 T→A 7.8% 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 7.8% 298.3 / 13.1 90L246Q (CTG→CAG) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (38/44);  new base A (2/5);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACAAA  >  NC_000913/3971848‑3972144
                                                                                                                                                    |                                                                                                                                                    
gAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCt                                                                                                                                                      <  4:193141/149‑1 (MQ=255)
 aaGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCt                                                                                                                                                      <  2:181618/148‑1 (MQ=255)
  aGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCtag                                                                                                                                                    <  3:908393/149‑3 (MQ=255)
   ggTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCt                                                                                                                                                      <  3:908138/146‑1 (MQ=255)
   ggTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCt                                                                                                                                                      >  4:908138/1‑146 (MQ=255)
   ggTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGc                                                                                                                                                   <  3:436235/149‑1 (MQ=255)
       aGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAAt                                                                                                                                               <  4:342465/149‑1 (MQ=255)
         tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATca                                                                                                                                             <  3:398776/149‑2 (MQ=255)
         tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATca                                                                                                                                             <  3:646100/149‑2 (MQ=255)
         tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATc                                                                                                                                              <  1:126105/148‑1 (MQ=255)
         tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATc                                                                                                                                              <  1:422761/148‑1 (MQ=255)
          gtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGt                                                                                                                                            >  3:936441/1‑149 (MQ=255)
          gtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGt                                                                                                                                            <  3:745169/149‑1 (MQ=255)
           tgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           <  3:201751/149‑1 (MQ=255)
           tgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           <  3:796989/149‑1 (MQ=255)
           tgtCGTCACCAACTCGCGGACGGTGGAAATGTGTAAGCTCACCGAAAACGGCTTCCGCGATGTGAATATCGCTTTCGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           <  3:806137/149‑1 (MQ=255)
            gtcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATcagt                                                                                                                                           <  1:64482/148‑4 (MQ=255)
             tcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCa                                                                                                                                          <  2:636572/148‑1 (MQ=255)
              cgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAc                                                                                                                                         >  1:689342/1‑148 (MQ=255)
                tcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGAGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCGccc                                                                                                                                       >  2:359582/1‑148 (MQ=255)
                   cTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTcgc                                                                                                                                   <  4:39744/149‑1 (MQ=255)
                   cTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTcg                                                                                                                                    >  2:675804/1‑148 (MQ=255)
                    tAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTcgcg                                                                                                                                  >  4:119363/1‑149 (MQ=255)
                        tCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTa                                                                                                                               <  1:835446/148‑1 (MQ=255)
                             gACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAtt                                                                                                                          >  1:43266/1‑148 (MQ=255)
                              aCGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttct                                                                                                                        <  3:341565/149‑1 (MQ=255)
                                     aaaTGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcc                                                                                                                  <  1:40001/148‑1 (MQ=255)
                                       aTGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctg                                                                                                                >  1:682579/1‑148 (MQ=255)
                                         gtgtAAGCTCACCGAAAACAGCTTCGGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCGGATTTGTGCCGATCAGGGGATTAACGTCTGGAAACTGATTCGCCGGGCGAACCGTCACCCTCGCGTTAATATTCTTCAGcctggc                                                                                                              <  2:128498/148‑1 (MQ=255)
                                            tAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccctg                                                                                                          >  4:927262/1‑149 (MQ=255)
                                             aaGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGtt                                                                                                                                >  4:78233/1‑126 (MQ=255)
                                             aaGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGtt                                                                                                                                <  3:78233/126‑1 (MQ=255)
                                               gCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCTTCACCCTCGCGTTAATATTCTTCAGcctggccctggc                                                                                                        >  2:868701/1‑148 (MQ=255)
                                                 tCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGt                                                                                                      >  2:579380/1‑148 (MQ=255)
                                                  cACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTgg                                                                                                    >  4:209676/1‑149 (MQ=255)
                                                    ccGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTggg                                                                                                   >  1:740477/1‑148 (MQ=255)
                                                     cGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGc                                                                                                  >  1:906794/1‑148 (MQ=255)
                                                      gAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCgg                                                                                                <  3:885712/149‑1 (MQ=255)
                                                       aaaaCAGCTTCCGCTATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGt                                                                                               <  4:894292/149‑1 (MQ=255)
                                                          aCAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCa                                                                                             <  1:152873/148‑1 (MQ=255)
                                                           cAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCAc                                                                                            <  2:1460/148‑1 (MQ=255)
                                                             gCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGc                                                                                         <  4:783404/149‑1 (MQ=255)
                                                              cTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAat                                                                                                                             >  1:725443/1‑112 (MQ=255)
                                                              cTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAat                                                                                                                             <  2:725443/112‑1 (MQ=255)
                                                                tCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCAtt                                                                                      >  3:845971/1‑149 (MQ=255)
                                                                      aTGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGtt                                                                                <  4:292537/149‑1 (MQ=255)
                                                                        gtgAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGa                                                                              <  4:213259/149‑1 (MQ=255)
                                                                         tgAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGa                                                                              <  1:415566/148‑1 (MQ=255)
                                                                               tCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGTCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTg                                                                        <  2:689342/148‑1 (MQ=255)
                                                                                cGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTgg                                                                       <  1:599470/148‑1 (MQ=255)
                                                                                cGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTgg                                                                       <  1:633194/148‑1 (MQ=255)
                                                                                   ttttGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTa                                                                   >  3:909969/1‑149 (MQ=255)
                                                                                       gCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCg                                                                <  1:230180/148‑1 (MQ=255)
                                                                                        cTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGt                                                               <  2:523872/148‑1 (MQ=255)
                                                                                           ttGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGACCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGc                                                            >  2:593001/2‑148 (MQ=255)
                                                                                            tGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGca                                                           <  1:556132/148‑1 (MQ=255)
                                                                                              aaTTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGcaca                                                         >  2:660716/1‑148 (MQ=255)
                                                                                                ttGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGaa                                                      >  4:267080/1‑149 (MQ=255)
                                                                                                   tCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTTTCGTg                                                              <  4:609098/138‑1 (MQ=255)
                                                                                                   tCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTg                                                              >  3:609098/1‑138 (MQ=255)
                                                                                                   tCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAccc                                                   >  3:335779/1‑149 (MQ=255)
                                                                                                   tCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAccc                                                   >  4:499512/1‑149 (MQ=255)
                                                                                                      cTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAccccc                                                 <  1:810215/148‑1 (MQ=255)
                                                                                                      cTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAccccc                                                 <  1:654758/148‑1 (MQ=255)
                                                                                                      cTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCca                                                <  4:406329/149‑1 (MQ=255)
                                                                                                         aTTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCcagca                                             <  4:845971/149‑1 (MQ=255)
                                                                                                         aTTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCcagc                                              <  2:267548/148‑1 (MQ=255)
                                                                                                                ccGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGcgcg                                       >  1:304276/1‑148 (MQ=255)
                                                                                                                ccGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCgg                                      <  3:927262/149‑1 (MQ=255)
                                                                                                                 cGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGc                                     <  3:25982/149‑1 (MQ=255)
                                                                                                                 cGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTACAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGc                                     <  4:600349/149‑1 (MQ=255)
                                                                                                                   aTCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCtt                                   >  3:915010/1‑149 (MQ=255)
                                                                                                                    tCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCtt                                   >  2:707119/1‑148 (MQ=255)
                                                                                                                      aGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTAt                                 <  1:67639/148‑1 (MQ=255)
                                                                                                                       ggggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATcc                               <  3:119363/149‑1 (MQ=255)
                                                                                                                       ggggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATcc                               >  3:310248/1‑149 (MQ=255)
                                                                                                                        gggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATcc                               >  1:895698/1‑148 (MQ=255)
                                                                                                                         ggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGt                                                                                 <  1:325469/97‑1 (MQ=255)
                                                                                                                         ggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGt                                                                                 >  2:325469/1‑97 (MQ=255)
                                                                                                                         ggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCg                              >  1:744237/1‑148 (MQ=255)
                                                                                                                           aTTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTa                            <  2:67586/148‑1 (MQ=255)
                                                                                                                            ttAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTAc                           >  2:605960/1‑148 (MQ=255)
                                                                                                                                cGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgcgc                      <  4:776581/149‑1 (MQ=255)
                                                                                                                                 gTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgcgcg                     >  3:835956/1‑149 (MQ=255)
                                                                                                                                  tCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgcgcg                     >  1:471929/1‑148 (MQ=255)
                                                                                                                                   cTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgcgcgc                    >  1:382180/1‑148 (MQ=255)
                                                                                                                                    tGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGa                  >  4:514013/1‑149 (MQ=255)
                                                                                                                                      ggAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTCGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAg                <  3:587946/149‑1 (MQ=255)
                                                                                                                                       gtgacgATTCGCCAGGCGAATCGGCACCCTCGCGGTAATATTCGTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCGGTTGATCCGGGGTTTATCGTGGCCCAGAACCCCCAGCAGGCGCGGCTTATCCGTGCCGCGCGCGAAGt               >  4:723455/6‑149 (MQ=255)
                                                                                                                                       gAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAg                >  2:914430/1‑148 (MQ=255)
                                                                                                                                           tGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGTGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAAc           <  3:353164/149‑1 (MQ=255)
                                                                                                                                            gATTCGCCAGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACg          >  4:720127/1‑149 (MQ=255)
                                                                                                                                              ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGa         <  1:359582/148‑1 (MQ=255)
                                                                                                                                               tCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGAt        >  1:393705/1‑148 (MQ=255)
                                                                                                                                                  ccTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATcac     >  2:371650/1‑148 (MQ=255)
                                                                                                                                                    tGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACaaa  <  4:97018/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAGGTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACAAA  >  NC_000913/3971848‑3972144

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