Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,996 T→A 5.7% 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 5.7% 439.3 / 12.9 123L246Q (CTG→CAG) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (61/55);  new base A (3/4);  total (64/59)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

Only 100 of 124 total aligned reads displayed.

GTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACAAA  >  NC_000913/3971852‑3972144
                                                                                                                                                |                                                                                                                                                    
gTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCg                                                                                                                                                  >  4:657722/1‑149 (MQ=255)
    gtgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATc                                                                                                                                              <  1:445068/149‑1 (MQ=255)
     tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCg                                                                                                                                             >  3:575681/1‑149 (MQ=255)
     tgtgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATca                                                                                                                                             <  3:556273/149‑2 (MQ=255)
       tgtCTTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATcagt                                                                                                                                           <  4:391445/149‑4 (MQ=255)
       tgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           <  1:206415/149‑1 (MQ=255)
       tgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           <  3:627124/149‑1 (MQ=255)
       tgtCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           <  2:119181/149‑1 (MQ=255)
         tcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAc                                                                                                                                         >  4:888077/1‑149 (MQ=255)
         tcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAc                                                                                                                                         >  1:361114/1‑149 (MQ=255)
         tcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAc                                                                                                                                         <  1:644659/149‑1 (MQ=255)
         tcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAc                                                                                                                                         >  4:410748/1‑149 (MQ=255)
           gtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAccc                                                                                                                                       >  1:320032/1‑149 (MQ=255)
              aCTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTcg                                                                                                                                    >  3:755102/1‑149 (MQ=255)
               cTAACTCGCGTACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGCAACTGATTCGCCTGGCGAATCGTCACCCTcgc                                                                                                                                   <  4:590270/149‑1 (MQ=255)
               cTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTcgc                                                                                                                                   >  4:444178/1‑149 (MQ=255)
               cTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTcgc                                                                                                                                   <  4:6005/149‑1 (MQ=255)
                  aCTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGCTTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGtt                                                                                                                                <  4:534064/149‑1 (MQ=255)
                  aCTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGtt                                                                                                                                <  2:256217/149‑1 (MQ=255)
                  aCTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGtt                                                                                                                                >  2:547030/1‑149 (MQ=255)
                    tCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTaa                                                                                                                              >  2:839685/1‑149 (MQ=255)
                    tCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTaa                                                                                                                              <  1:36004/149‑1 (MQ=255)
                        ggACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAtt                                                                                                                          <  1:599/149‑1 (MQ=255)
                         gACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttc                                                                                                                         <  3:437438/149‑1 (MQ=255)
                          aCGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttct                                                                                                                        <  2:320032/149‑1 (MQ=255)
                            ggcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttcttc                                                                                                                      >  2:386568/1‑149 (MQ=255)
                            ggcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttcttc                                                                                                                      >  3:48019/1‑149 (MQ=255)
                             gcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCa                                                                                                                     <  1:902420/149‑1 (MQ=255)
                             gcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCa                                                                                                                     <  1:419197/149‑1 (MQ=255)
                              cggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAg                                                                                                                    >  1:892683/1‑149 (MQ=255)
                              cggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAg                                                                                                                    >  1:236265/1‑149 (MQ=255)
                               ggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGc                                                                                                                   <  1:839685/149‑1 (MQ=255)
                                gAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcc                                                                                                                  >  2:700141/1‑149 (MQ=255)
                                    tgtgTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggc                                                                                                              <  4:154227/149‑1 (MQ=255)
                                     gtgtAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggcc                                                                                                             >  1:647207/1‑149 (MQ=255)
                                      tgtAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATcagt                                                                                                                                           >  2:372402/1‑115 (MQ=255)
                                      tgtAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATcagt                                                                                                                                           <  1:372402/118‑4 (MQ=255)
                                       gtAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccct                                                                                                           <  4:48019/149‑1 (MQ=255)
                                       gtAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTCTGCCGAGCAGGGGATTAACGTCTGGGAACTGATTCGCCTAGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccct                                                                                                           >  4:600039/1‑149 (MQ=255)
                                           gCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATcagtt                                                                                                                                          >  1:129953/1‑110 (MQ=255)
                                           gCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAATcagtt                                                                                                                                          <  2:129953/114‑5 (MQ=255)
                                             tCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTg                                                                                                     >  3:617512/1‑149 (MQ=255)
                                             tCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTg                                                                                                     >  4:118329/1‑149 (MQ=255)
                                              cACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           >  1:126740/1‑110 (MQ=255)
                                              cACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTc                                                                                                                                           <  2:126740/110‑1 (MQ=255)
                                              cACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTgg                                                                                                    <  4:409668/149‑1 (MQ=255)
                                                ccGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGc                                                                                                  <  2:749543/149‑1 (MQ=255)
                                                 cGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAg                                                                                                                    <  3:219394/130‑1 (MQ=255)
                                                 cGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAg                                                                                                                    >  4:219394/1‑130 (MQ=255)
                                                  gAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCgg                                                                                                <  3:484577/149‑1 (MQ=255)
                                                   aaaaCAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGt                                                                                               <  3:657722/149‑1 (MQ=255)
                                                     aaCAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCa                                                                                             <  1:700141/149‑1 (MQ=255)
                                                           ttCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCAt                                                                                       <  2:606748/149‑1 (MQ=255)
                                                               gcgATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCt                                                                                   <  3:335638/149‑1 (MQ=255)
                                                                      gAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATc                                                                            <  1:386568/149‑1 (MQ=255)
                                                                       aaTATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATcc                                                                           >  1:417068/1‑149 (MQ=255)
                                                                         tatCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTAGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGt                                                                         >  4:556218/1‑149 (MQ=255)
                                                                          atCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGt                                                                                                      <  1:461245/119‑1 (MQ=255)
                                                                          atCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGt                                                                                                      >  2:461245/1‑119 (MQ=255)
                                                                            cGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGt                                                                      >  3:258892/1‑149 (MQ=255)
                                                                             gCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGtt                                                                     >  1:243712/1‑149 (MQ=255)
                                                                              cTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGttt                                                                    <  2:643585/149‑1 (MQ=255)
                                                                              cTTTTGCTAATAAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGttt                                                                    >  4:330235/1‑149 (MQ=255)
                                                                                   gCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGTCGAAGCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGt                                                               <  2:668735/149‑1 (MQ=255)
                                                                                   gCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGt                                                               <  3:888077/149‑1 (MQ=255)
                                                                                    cTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTg                                                              <  2:390898/149‑1 (MQ=255)
                                                                                     tAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTgg                                                             >  2:653597/1‑149 (MQ=255)
                                                                                     tAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTgg                                                             >  3:181923/1‑149 (MQ=255)
                                                                                     tAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTgg                                                             >  4:214966/1‑149 (MQ=255)
                                                                                      aaTGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGc                                                            >  1:340437/1‑149 (MQ=255)
                                                                                        tGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGtt                                                                                >  4:360706/1‑127 (MQ=255)
                                                                                        tGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGtt                                                                                <  3:360706/127‑1 (MQ=255)
                                                                                           aTTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGa                                                       >  2:518842/1‑149 (MQ=255)
                                                                                              gTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAcc                                                    <  2:773171/149‑1 (MQ=255)
                                                                                                 gCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAAccccc                                                 >  3:87977/1‑149 (MQ=255)
                                                                                                   tGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGTGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCcag                                               <  4:575681/149‑1 (MQ=255)
                                                                                                         gtgCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGgcg                                         >  2:157873/1‑149 (MQ=255)
                                                                                                              gATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCt                                    >  3:428900/1‑149 (MQ=255)
                                                                                                                  aGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATc                                >  3:655913/1‑149 (MQ=255)
                                                                                                                  aGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATc                                <  4:755102/149‑1 (MQ=255)
                                                                                                                      gATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTa                            >  4:133414/1‑149 (MQ=255)
                                                                                                                             gTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgcgcg                     >  3:626541/1‑149 (MQ=255)
                                                                                                                              tCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgcgcgc                    >  3:89129/1‑149 (MQ=255)
                                                                                                                              tCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgcgcgc                    >  4:257547/1‑149 (MQ=255)
                                                                                                                                 gggAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGaa                 >  1:179735/1‑149 (MQ=255)
                                                                                                                                    aaCTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGcaca                                                         >  3:864513/1‑106 (MQ=255)
                                                                                                                                    aaCTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGcaca                                                         <  4:864513/106‑1 (MQ=255)
                                                                                                                                    aaCTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTg              >  1:135178/1‑149 (MQ=255)
                                                                                                                                    aaCTGATTCGCCTGGCGAATCGTCACACTCGCGTTAATCTTCTTCAGCCTGGCCCTGGCGTGGGCGTTCACTGCATTGCTGTTGATCCGTGGTTTATCATGGCACAGATCCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTg              >  3:776825/1‑149 (MQ=255)
                                                                                                                                     gacgATTCGCCAGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGa             >  3:465641/4‑149 (MQ=255)
                                                                                                                                        gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCcagc                                              >  4:247945/1‑113 (MQ=255)
                                                                                                                                        gATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCcagc                                              <  3:247945/113‑1 (MQ=255)
                                                                                                                                         aTTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACt                                                                                           >  3:408753/1‑67 (MQ=255)
                                                                                                                                         aTTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGa         >  4:509704/1‑149 (MQ=255)
                                                                                                                                         aTTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGa         >  1:644347/1‑149 (MQ=255)
                                                                                                                                          ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGt                                                               <  4:864923/94‑1 (MQ=255)
                                                                                                                                          ttCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGAt        >  3:847235/1‑149 (MQ=255)
                                                                                                                                           tCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATc       >  3:45345/1‑149 (MQ=255)
                                                                                                                                            cGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATca      <  2:135178/149‑1 (MQ=255)
                                                                                                                                                tGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACaaa  <  3:459409/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTGAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGATCACAAA  >  NC_000913/3971852‑3972144

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: